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FOREWORD 


METALLURGY is an ancient profession and one in which science has 
played an increasingly important part. Metallurgists had acquired a vast 
body of empirical knowledge and practical skill long before the develop¬ 
ment of the formal sciences requisite to its intellectual understanding. For 
centuries the chemist learned from the metallurgist; but in the nineteenth 
century the tables were turned: the acceleration of metallurgical progress 
came from the application of chemical science. The methods of analysis 
and the chemical preparation of pure materials made possible, for the first 
time, realistic studies of the effect of composition and impurities on tech¬ 
nical behavior and properties. The twentieth century has seen the growth 
of X-ray diffraction and other methods of the physicist and their applica¬ 
tion to metals. It seems likely that the important developments of the 
next generation will come about by the application of the principles of 
physics, chemistry, and mathematics to the study of structure, properties, 
and behavior of metals. 

Late in 1945 the University of Chicago established three institutes for 
research in fields involving more than one established discipline. One of 
these was the Institute for the Study of Metals. Just as in the new fields 
of nucleonics and biophysics, so also in the old field of metallurgy there 
seemed to be need for an organization in which, under a common roof and 
exposed to mutual stimulation, physicists, chemists, and more specialized 
scientists could attack a broad problem. It is mutually advantageous for 
the physicists to learn and utilize some of the large body of rather em¬ 
pirical metallurgical information and for the metallurgist to employ the 
rigor of the methods of the physicist and physical chemist. 

The advance of knowledge in any field is closely related to its literature, 
for personal contact is too limited in both time and space. Soon after the 
Institute was established, it became obvious that its aims could be fur¬ 
thered by monographs on various aspects of the science of metals, and a 
series was established under the joint auspices of the Institute and the 
University of Chicago Press. Works by any author are eligible, provided 
that they deal in a competent manner with some aspect of metal science or 
closely related fields. Textbooks and discussions of industrial “know-how” 
are excluded because excellent outlets for such material already exist. 
Three kinds of monographs are to be encouraged: critical reviews of par- 
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ticuJar topics that seem ripe for research—the kind of work that will gen¬ 
erally be written by a scientist as he makes a detailed survey of the status 
of a field before undertaking intensive original work on it; extended sum¬ 
maries written after prolonged work in a specialized field at the time when 
a scientist wishes to record his work in perspective as a basis for future 
investigations by himself and others; and, lastly, books discussing the ap¬ 
plication of particular scientific principles or techniques to metallurgy and 
related subjects. 

The present book, Elasticity and Anelasticity of Metals, by Clarence 
Zener, deals with a subject on the border line of metaUurgy, physics, and 
engineering. Very largely as a result of Dr. Zener’s work, internal friction 
has become an important tool in the study of the structure and behavior 
of metals. There are already many different types of anelastic deformation 
kno^vn; and it seems likely that the acoustic spectrum will come to occupy 
a position in solid-state studies comparable to that possessed by the opti¬ 
cal spectrum in the theory of atomic and molecular structure. Though 
some of the data in this book have appeared before, this is the first time 
that a connected and thorough treatment of the whole subject has been 
written, and it should prove of considerable value to physicists, metal¬ 
lurgists, and engineers alike. The section on microelasticity, hitherto un¬ 
published, will be found to have considerable relevance to the stability and 
transformation of metallic phases. This is a work that the committee is 
particularly happy to sponsor as the first member of the series. 

Cyril Stanley Smith, Chairman 
Charles S. Barrett 
Earl A. Long 

Committee on Monographs 
Institute for the Study of Metals 
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TPL\T metals manifest deviations from perfect elastic behavior, even at 
small stress levels, has been known for over a centur>'. Only during the 
past decade, however, has marked progress been made in the interpreta¬ 
tion of these deviations in terms of the microstructure of the metal. Al¬ 
though not many years ago these deviations were investigated only by 
erudite students of metals and were regarded by metallurgists as interest¬ 
ing but quite useless, today they are being studied to answer very definite 
questions regarding the kinetics of changes in the microstructure. In pre¬ 
paring this monograph the author has attempted to present in a logical 
manner the science of the nonelastic behavior of metals at low stress 
levels, to review all the pertinent work in this field, and to point out the 
role which studies of this behavior will play in the future development of 
metallurgical science. Anticipating that all who may want to profit from 
this monograph will not wish to follow the formal development of the 
theory, the author has endeavored to write in such a manner that the for¬ 
mal developments may be omitted by an intelligent reader without loss of 
a physical understanding of the subject. The elastic and nonelastic prop¬ 
erties of metals are so interrelated that it has seemed inappropriate to dis¬ 
cuss the latter at length without at least a brief introduction to the former. 

Nonelastic behavior carries with it a connotation of a permanent set 
after removal of all stress. Although metals behave nonelastically even at 
low stress levels in the sense that stress and strain are not single-valued 
functions of each other, e.g., strain may lag behind stress in periodic vibra¬ 
tions, no permanent set remains after removal of all stresses. It has there¬ 
fore seemed appropriate to introduce a new term that is free of this con¬ 
notation of permanent set. Anelasticity has been chosen to denote that 
property of solids in virtue of which stress and strain arc not single-valued 
functions of one another in that low stress range in which no permanent 
set occurs and in which the relation of stress to strain is still linear. 

The author is indebted to the Office of Naval Research for its con¬ 
tinued support of research upon the anelasticity of metals, a support 
which has aided materially in the preparation of this monograph. 


vn 





TABLE OF CONTENTS 


Introduction 


1 


PART ONE. ELASTICITY OF METALS 

I. Formal Relations between Stress and Str.\in ... 7 

A. Strain. 7 

B. Stress. 9 

C. Relation between Stress and Strain . . 11 

D. Effect of Crystal Symmetry . . .. .14 

II. Low-Temperature Elastic Constants IN Cubic Metals ... 16 

A. Review of Data.16 

B. Relation of Elastic Constants to Interatomic Forces and to Structure 19 

III. Temperature Dependence of Elastic Constants.24 

A. Single-Particle Model , . 26 

B. Real Crystals. 29 

IV. Microelasticity. 32 

PART TWO. ANELASTICITY OF METALS 

V. Formal Theory OF Anelasticity.41 

A. Generalizations of Elastic Equations.41 

B. Boltzmann’s Superposition Principle.48 

C. Relaxation Spectrum. . . .55 

VI. Measurement of Relaxation Spectrum ... .60 

A. Relation between Measures of Internal Friction . 60 

B. Measurement of Internal Friction.63 

C. Variation of Frequency.66 

VII. Physical Interpretation of Anelasticity.69 

A. Homogeneous Relaxation.69 

1. Thermodynamical Variables.72 

II. Orthogonal Thermodynamic Potentials.76 

III. General Theory of Relaxation by Diffusion.79 

IV. Relaxation by Thermal Diffusion.89 

V. Relaxation by Atomic Diffusion.96 

VI. Relaxation by Magnetic Diffusion.101 

VII. Relaxation of Ordered Distributions . . ... 106 

VIII. Relaxation of Preferential Distributions.Ill 


IX 






























X 


CONTENTS 


B 


Index 


. Inhomogeneous Relaxation. 

I. Concept of Two-Component Systems 

II. Stress Relaxation along Previously Formed Slip Bands 

III. Stress Relaxation across Grain Boundaries 
IV • Stress Relaxation across Twin Interfaces 

INDEX 


126 

126 

132 

147 

159 


167 



INTRODUCTION 


THE modern theory of elasticity may be traced to a law formulated by 
Hooke’ in the seventeenth centur>': “Ut tensio sic vis.” According to this 
law, if a “force” is applied to a body, the resulting deformation is propor¬ 
tional thereto. 

The “force” in Hooke’s law may be an actual force or any related quan¬ 
tity, such as torque, couple, etc. The deformation is usually so defined 
that the work done upon the body when the deformation changes by a 
small amount is equal to the force times the increment of the deformation. 
Thus, if F is the force, the deformation, Z), is so defined that 

5 work = F* hD . 

Thus, for a force, torque, or couple, D refers to displacement, angle of 
twist, and angle of bend, respectively. 

The mathematical formulation of Hooke’s law is 

MD = F (i) 

where the coefficient A/, known as the “elastic modulus,” is independent 
of the force F. The modulus does, however, depend upon the type of the 
applied force. Thus M may refer to the tensile modulus, the torsion modu¬ 
lus, or the bending modulus. 

As originally formulated, Hooke’s law is applicable only to quasi-static 
forces, i.e., to forces very slowly applied. This limitation is evident when 
it is realized that the deformation, D, in Hooke’s law is influenced by 
distortions throughout the entire body. The deformation cannot, there¬ 
fore, attain its value associated with the applied force until, in effect, 
every part of the body has been notified of the force and the resulting dis¬ 
tortions have been signaled back to the region of application of the force. 
It is evident that the equilibrium deformation cannot be established in a 
time less than the time required for an elastic wave to travel from the re¬ 
gion of application of the force to the most remote part of the body and 
back again. As an example, an impulsive transverse force applied to a 
taut string produces a rectangular wave which extends rapidly toward the 
supports. In place of Hooke’s law, the relation between force and displace¬ 
ment in this case is 

^-= Constant XF. 
at 

1. Robert Hooke, De potenlia restitutiva (London, 1678). 

1 
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This equation is valid until waves have been reflected back from the end 
supports. Hooke’s law may not be applicable even after elastic waves have 
been reflected from the most remote parts of the specimen. The elastic 
waves may continue to travel to and fro. In other words, the system may 
be bCt into vibration by the force. A force may be regarded as applied in a 
quasi-stadc manner, with no attendant vibrations, only if the time taken 
to apply it is long compared with the period of the slowest mode of vibra¬ 
tion. In order that the original formulation of Hooke’s law may be ap¬ 
plicable, the force must therefore be applied the more slowly, the larger 
the body. 

During the early part of the nineteenth century’ Hooke’s law was gen¬ 
eralized in such a manner as to remove the restriction that forces must be 
applied in a quasi-static manner. In this generalization the body is con¬ 
sidered as partitioned into a great number of elementary regions of very 
small extent. The distortion of each region is then determined solely by 
the forces acting upon it by the neighboring regions, irrespective of how 
rapidly these forces change. In effect, one regards the elementary regions 
as so small that all forces change only a negligible amount during a time 
interval equal to the period of their slowest mode of vibration. The change 
in shape of the elementary regions may be completely described in terms 
of certain quantities called “strain components.” The generalized Hooke’s 
law then relates the strain components to the force components acting on 
the elementary regions or, more properly, to the surface density of the 
force components, known as “stress components.” The variation of the 
stress throughout the body is determined by this generalized Hooke’s law 
together with the equations of motion of each elementary region. 

As an example of the above generalization, reference will be made to 
the simple case of a force applied to one end of a wire, the other end of 
which is rigidly fixed. When a force F is applied to the end of the wire, the 
change in length of the wire will be proportional to F only if the force is 
applied quasi-statically. According to the above-outlined method, one 
focuses attention on an element of the wire of original length bX. Since 

is arbitrarily small, the acceleration of the elementary region will be 

finite only if the forces acting upon its two ends are equal in magnitude 

and opposite in direction. These forces will be denoted by/. The strain e of 

the element is defined in terms of the original and final coordinates X 
and U, 

dU , 


Then 


f = Constant X e 


(iii) 
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is the generalized Hooke’s law appropriate for this special case. The varia¬ 
tion of the force / along the specimen is given by 


mdW _ df 

dr- dX ’ 



where m is the mass per unit original length. Equations (ii)-(iv), together 
with the appropriate boundary conditions, completely specify the re¬ 
sponse of the wire to a force F applied at its end in an arbitrary manner. 

The classical theory of elasticity is essentially the development of all 
the consequences arising from the application of Hooke’s law to elemen¬ 
tary regions. Many excellent treatises- have been written on this classical 
theory, as well as an extensive history* of its development. In the discus¬ 
sion of elasticity presented in this monograph we are concerned only with 
the coefficients which relate stress to strain and with the dependence of 
these coefficients on the various physical parameters. 

The authors of the classical theory- of elasticity were under no illusion 
as to the lack of strict applicability of their fundamental assumptions to 
real solids. The value of the classical theory' lies not in its precise descrip¬ 
tion of the behavior of solids under applied forces but in its description of 
this behavior with sufficient accuracy for most practical purposes. 

According to the classical theory of elasticity, stress and strain are 
uniquely related. Also applied force and deformation are uniquely re¬ 
lated, provided that the force is applied in a quasi-static manner, i.e., 
provided that it is varied so slowly that no vibrations result. Weber/ as 
early as 1825 , found, in his investigations upon galvanometer suspensions, 
small deviations from perfect elasticity. Upon release of the couple, the 
suspension did not at once return to its zero-point but approached this 
zero-point only asymptotically. Such behavior he called an “elastic after¬ 
effect” (elasUsche Nachwirkung), an early summaiy- of which was pre¬ 
sented by Auerbach.® Other effects which are contrary to the theory of 
elasticity but which are shown by real solids are internal friction, stress re¬ 
laxation, and variation of modulus with frequency of measurements. All 

2. E.g., A. E. H. Love, Malhemctical Theory o/Elasticity (4th ed.; Cambridge: Cambridge 
University Press, 1927J; R. V. Southwell, Theory of Elasticity (London: Oxford University 
Press, 1936). 

3. I. Todhunter and K. Pearson, History of the Theory of Elasticity (Cambridge: Cambridge 
University Press, 1886). 

4. W. Weber, “Uber die Elastizitat der Seidenfaden,” Poggendorf's Ann., XXXV (1834), 
247; and “Uber die Elastizitat fester Kdrper,” ibid., XXIV (2d ser., 1841), 1. 

5. F. Auerbach, “Elastische Nachwirkung,” H'inkehnann's Handbiich der Physik, I (Bres¬ 
lau, 1891), 769-831. 
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such effects are different manifestations of the lack of uniqueness of the 
relation between stress and strain. The property of a solid in virtue of 
which stress and strain are not uniquely related in the pre-plastic range 
is called “anelasticity.” In the discussion of anelasticity presented in this 
monograph we are concerned primarily with the physical origin of the 
various anelastic effects. 
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ELASTICITY OF METALS 




I 

FORMAL RELATIONS BETWEEN STRESS AND STR.\IN 


NO ESSENTLAL additions have been made to the formal theory of the 
elasticity of crj'stals since the publication of Voigt’s^ excellent treatise. 
This formal theoiy^ will therefore be only briefly reviewed, with no attempt 
at completeness. 

A. STRAIN 

The concept of strain components has been introduced in order that 
small deformations of a body may be uniquely specified. In preparation 
for a definition of the strain components, we shall first define the dis¬ 
placement vector. As a result of deformation, a particle which was origi¬ 
nally at X, y, z will now find itself at a: + I/, y + F, 2 + IF. The vector 
{U, F, IF) is called the “displacement vector,” and its components are 
functions of x, y, and 2. We now confine our attention to a small element 
of volume, as indicated in Figure 1 , and regard the origin of coordinates 
as lying within the element before deformation. 

The displacement of every part of this small element may then be ex¬ 
pressed in terms of twelve constants, namely, as 


,, ,, , ac/ , dU _,dU 

+ + etc., 

where the derivatives are taken at the origin. The vector {Uq, Fq, IFo) 
gives the displacement of the origin of coordinates. The displacement of 
every part of the body with respect to the origin, (U\ V', W'), is then 
given in terms of nine constants, namely, their nine partial derivatives. 

Under certain important conditions a more direct physical interpreta¬ 
tion may be ascribed to nine other constants which may be formed by a 
linear combination of the above nine. These are as follows: 



du 


dF 


dW 


dx ’ 

Cffl, — 

dy ■ 


dz 


dW , 

dV 


dU ^ 

aiF 

Cyi = 

dy 

dz ’ 


= -Z— 
dz 

dx ' 


dV 

dU 


_dW 

dV 



dy' 


dy 

dz ■ 


dU 

dW 


_dV 

dU 

Wy = 

dz 

dx ' 


dx 

dy- 



1 . W. Voigt, Lehrbitth der Kristallphysik (Berlin; Teubner, 1910). 
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These nine constants may be given direct physical interpretations when 
they are small-order quantities, i.e., very small compared with unity. 
Under this condition the deformation relative to the origin may be con¬ 
sidered as taking place in two steps. The first deformation is represented 
by the displacement vector, 


Us = + § e^y + | , ’ 

Ts ‘2 “f“ ^VVy H” 2 ^ » I 

116'= 2 5 -\- e„z , . 




Fig. 1.—Illustration of displacement vector in two dimensions 

in which the coefficients form a symmetrical matrix. The second deforma¬ 
tion is represented by the displacement vector, 

UA ~ ) 1 

V A ~ tWjA' ' f ^ ( 3 ) 

11 .4 ” -}- 2^Jcy f ^ 

in which the coefficients form an antisymmetrical matrix. 

It may readily be shown that the displacement vector {Us, Vs, HV) 
represents a distortion for which a particular set of coordinate a.xes at¬ 
tached to the body do not change direction during the deformation. These 
axes are known as the “principal axes’’ of the distortion. The six quanti¬ 
ties, e«, Cyy, Cti, eyt, Cix, aud Ctv, are known as the “strain components.” The 
first strain component may be regarded as the relative change in length 
of a line originally parallel to the a;-axis; and an analogous interpretation 
may be given the second and third strain components. The strain com¬ 
ponent Cyi may be regarded as the change in the angle between the y- and 
the 2-axis during deformation if these axes are regarded as fixed to the 
body during the distortion; and a similar interpretation may be given 
to 

The numerical values of the strain components depend upon the orien- 
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tation of the coordinate axes to which they refer. It is often desired to 
express the strain components referred to one set of axes in terms of the 
strain components referred to a second set of axes. Let the transforma¬ 
tion between the two sets of axes be specified by the following scheme: 


Thus 

and, conversely, 



X 

y 

2 

x' 

h 

mi 


y' 

h 

nti 

«2 

z' 

h 


«3 


x' = hx-j- miy-j-HiZ 
X = lix'l^y'liZ' . 



The transformation equations between the two sets of strain components 
are then obtained by transforming the set of equations (2) into one con¬ 
taining Usj Vs, and x\ y\ and z'. Thus we multiply the three equa¬ 
tions by /i, 4 , and /a, respectively, and then form their sum, thereby ob¬ 
taining Us in terms of x, y, and 2. We now transform x, y, and z into 
x\ y', and t! . The result is 


^zz' = eyym\->r e„n5-f ey,wiWi-f 


^ ; I 


f 





= 2 eMz-\- 2 evyW2W3-|- 2 e„«2W3+ fiyr (W2«3+ W3«2) 

+ (H2/3’4‘«3f2) 4* I etc. 

It may readily be shown that the displacement vector {Ua, Wa) 
represents a deformation in which the body suffers a rotation as a rigid 
body about its origin. The quantities Wx, ojy, and Wx are the three com¬ 
ponents of the rotation vector and therefore transform with a rotation of 
coordinate axes according to scheme ( 4 ). 



B. STRESS 

The concept of stress has been introduced in order that we may specify 
completely the traction which acts across all planes passing through a 
given point within a solid. In preparation for the definition of stress com¬ 
ponents we shall consider the forces acting across an element of surface 
within the body. In particular, we shall consider the forces with which the 
material on one side of an element of internal area, dA, acts upon the 
material on the other side. This force will be proportional to dA, and may 
therefore be written as TdA, as illustrated in Figure 2 . Since the sign of T 
will depend upon which side we are considering as exerting the force, am¬ 
biguity will be removed by assigning to T the subscript n, which denotes 
the unit vector normal to dA pointing toward the region which we regard 
as exerting the force. This same notation also specifies the orientation of 
the element of area dA. 
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Upon applying to the tetrahedron of Figure 3 the condition that the 
total force acting upon a body in equilibrium must be zero, one finds that 

4 

the force density, Tn, may be expressed in terms of the force densities 
acting across the xy-, yz-, and sx-planes. If /, and n refer to the direc¬ 
tion cosines of n, one obtains for the components of Tn, 


.Yn=/.Y,+ ?wX, + «X, , 1 


2n = IZx “F fuZy “j“ ttZ. • J 





Fig. 2.—Illustration of traction 

Fig. 3.—Method of resolving traction across an arbitrarj- plane 

The components of force density AT,, . . . , are called “stress components." 
The number of independent stress components may be reduced to six 
upon applying to an elementar>' cube the condition that the total mo¬ 
ment acting thereon must be zero. One finds 

= , Xy =\\. (8) 

The independent stress components may therefore be taken as Xx, Yy, 

Zx, and Xy. The traction across an arbitrary plane may be e.xpressed 
in terms of these six stress components through equations ( 7 ) and (8). 

The numerical value of the stress components will depend upon the 
particular reference Cartesian axes. The relation between the stress com¬ 
ponents referred to different axes may be obtained by the use of equa¬ 
tions ( 7 ). We shall consider two sets of axes related by scheme ( 4 ). The 
quantity Xx> is obtained by letting n be directed along the x'-axis, multi¬ 
plying the first of equations ( 7 ) by h, the second by wi, the third by «i, 
and forming the sum. We obtain 

Xx- = l\Xx+ m\ Yy-\-n\Zx-y 2 F, + 2 n,hZx + llifn^Xy , 

with two analogous equations for Yy- and Z,^; and 

Xy' = lihXx+m\m2Yy-\-nxn‘iZt-\- »<2Wi) F, 

"F (m i/j “F H2/1) ”F Iztn^ Xy 

with two analogous equations for Y[' and Z^'. 


( 9 ) 
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C. RELATION BETWEEN STRESS AND STR.\IN 


We have seen that the deformation in any elementary region is com¬ 
pletely specified by the six strain components, exx, Cw, c,,, e.x, and 
We have further seen that the traction across an arbitrary plane passing 
through the elementary region is specified completely by the six stress 
components, Xx, Vy, Z*, and Xy. The generalized form of Hooke’s 

law states that the strain components are linear functions of the stress 
components and, conversely, that the stress components arc linear func¬ 
tions of the strain components. Thus, on the one hand, we have 


^xx = -^11 A'x + Ji2 l yH- SiiZ, + Fi + ■Tii^x + ' 

eyy= isi-^xH" -522 Ty “f* S 2 yZt -f ^24 F, -f ^ 26 ^,-f- S^Xy , 


^XV — •yBl-^’^x H“ -^62 I V + ^63^1 H" -^64 I x "I" ■^efr^x + ^66-^^ i ^ 


( 10 ) 


and, on the other hand. 


Xx = Cu Cxx*h CijCyy-f- CiiCit -|- Ch Cys + C 15 C^x + C16 ^xV » 1 
Fy = C21^XX+ C22^yyH" C23^xr+ ^ 24 + r2S^xx+ C%^iv t 


Xy — C61 ^xx "1“ ^63 H“ ^64 H“ ^65 “(“^66 • J 


( 11 ) 


The quantities in, , are called the “elastic constants,” the quanti¬ 

ties Cii, Ci2, . . . , the “elastic coefficients.” It may readily be verified that 
the matrix of the elastic constants, ||5||, and the matrix of the elastic 
coefficients, 
sense that 


, may be regarded as the reciprocal of each other in the 


c = I 1 


( 12 ) 


where j|l is the unit matrix; thus 


5v;:t 


and 


(13) 


iVi’l 


(14) 


The thirty-six elastic coefficients arc not all independent. In order to 
derive those relations that hold for all crystals, we observe that an in¬ 
crement of the strain-energy density, t/, is given by 

dU — X^de^x-h Fyrfcyy + Zx(icx, 4 - K.rfey,-j-Z*rfCxx-f XyrfCxv• ( 15 ) 
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Since the strain energ>* may be regarded as a single-valued function of the 
strains, it follows that 

dX, dVy 

d Cyy d €xz ' 

and hence, from equations (11), 



Similarly, one finds that, in general, 

C jk ~ Ckj • 


( 16 ) 


Upon forming the differential oi U — + V^uv + ■ • in place of 

Uy it follows that 

Sjk-Skj. ( 17 ) 

The matrices ||j|| and ||c|| are therefore symmetrical. These symmetry 
relations reduce the number of independent constants from thirty-six to 
twenty-one. 

The numerical value of the elastic constants or coefficients depends 
upon the particular Cartesian coordinates to which they refer. It is some¬ 
times desirable to express the elastic moduli with respect to one Cartesian 
system in terms of those referred to in another Cartesian system. This re¬ 
lation may readily be obtained by means of equations ( 5 ), (6), and ( 10 ), 
and of equations analogous to equations ( 9 ), which express the ...» 
stress components in terms of the X'x> stress components. These latter 
equations may be derived in the same manner as were equations ( 9 ). 
They are 

X. = llx',' + llYl’ + lX' + 2t2hr',' + 2hUZ',‘+2likX'.',etc. , 
and 

Xj, = /iWiXl'-|-/ 2 m 2 K/ 4 -/ 3 W 3 ^ 1 '+ K/ f 

+ (/3W1-I-/1W3) Zx' -j- (/1W2 + /2W1) Xv'• 


Several examples will be given for crystals of cubic symmetry, in which, 
as is shown later, the elastic constants are specified by the following 
matrix: 


Sn 

S 12 

S 12 

0 

0 

0 

Sx2 

Six 

S\2 

0 

0 

0 

S 12 

Si2 

Sll 

0 

0 

0 

0 

0 

0 

544 

0 

0 

0 

0 

0 

0 

■^44 

0 

0 

0 

0 

0 

0 

•S 44 
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a) As a first example, we shall obtain the constant s[i referred to the 
x'-, y'-, and s'-axes in terms of the constants referred to the x-, y-, and 
s-axes. According to equation (5), we multiply the first of equations (10) 
by l\, the second by m\, etc., and then form the sum, thereby obtaining e'^. 
We then transform the stresses in the right-hand side according to equa¬ 
tions (18). The coefficient of defined as s'u, is given by 




(19) 


b) As a second example, we shall find the elastic shear constant for the 
(110)[U0] shear of a cubic crystal. In this case our table of direction co¬ 
sines is as follows: 


X 


.f 


X 


V2 V2 ^ 


1 


0 


1 


V2 V2 


= 0 


0 1 


Upon applying this table to equation (6), we obtain 

— €-1 — Cyy . 


Now, when we express the strains in terms of the stresses by means of the 
|s|l matrix for a cubic crystal and then express the stresses referred to 
xyz-axes into stresses referred to x'y'z'-axes by means of equations (18), 
we obtain 


Cy'/ =* 2 (^ 11 - 512) Yi 


and hence 


^44 — 2 (ill — ^ 12 ) . 


( 20 ) 


c) As a third example, we shall find the shear constant corresponding 

to the (111)[211] shear of a cubic crystal. Our table of direction cosines 
now is: 



X 

y 

z 

x'[1121 

1 

VI 

1 

2 

V 6 

y'lllll 

1 

V3 

1 

V3 

1 

^/3 

1 

2'[1T0] 

1 

V2 

1 

“ V2 

0 
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Proceeding as before, we obtain 

Cz 1/ ! 2 ( “b 2 C,.) “b ( 2 €y^ I 

= i 4-2 5 ) A'i,- + -^= — s) (Z,-— x'^') , 

where 

We therefore obtain 

^,= 5l^. + 4(.„-5,J I. (21) 

A review of the methods of determining the elastic constants has 
been given by Hearmon.- 

D. EFFECT OF CRYSTAL SYMMETRY 

As previously mentioned, the values of the elastic coefficients vary 
with the orientation of the coordinate axes with respect to the crystal. If 
one of the coordinate axes coincides with an axis of symmetry of the crys¬ 
tal, the relative orientations of the coordinate axes and the crystal are un¬ 
changed by an appropriate rotation of the coordinate axes about the axis of 
symmetry. The elastic coefficients must therefore also remain unchanged. 
This invariance is possible only if they satisfy certain conditions given 
in the accompanying table for simple symmetry operations.^ 


RESTRICTIONS PLACED UPON ELASTIC COEFFICIENTS 

BY SYMMETRY RELATIONS 



Cvi 

^13 

0 

0 

Cie 

C|i 

C\2 


Cu 


0 

C\z 

Cri 

Cjs 

0 

0 

Ces 

Ci2 

C\i 

ClS 


C-ii 

0 



Cxt 

0 

0 


Cij 

C\2 

Cj3 

0 

0 

0 

0 

0 

0 

Cu 


0 

Cli 


0 

Cu 

0 

C^'a 

0 

0 

0 

Cis 

CiS 

0 



0 

0 

<•« 



c« 

Cu 

0 

0 

1 

: 1 

0 

0 

0 



C\\ ^C\t 






C2S 

Cu 

2 


s-axis has a twofold symmetry^ or reRec- r^axis has a threefold symmetry* 

tion symmetry exists upon xy-planc 


Ci\ 

Ci2 

^12 

0 

0 

1 

Ci€ 

1 

! ^11 

Cli 

Cu 

0 

0 

1 

0 

' C\2 

Cn 

ClJ 

0 

0 


Cu 

Cu 

C}2 

0 

0 



Cii 

c» 

0 

0 

0 


Cli 

Ci^ 

0 

0 

0 ' 

0 

0 

0 

C*4 

0 

0 

^ 0 

0 

0 

Cm 

0 

0 

0 

0 

0 

0 

c» 

0 

0 

0 

0 

0 

Cm 

0 ' 

Cie 


0 

0 

0 

Cm 

0 

0 

0 

0 

0 

Cii —C|J 

2 


^axis has a fourfold symmetry r<axis has a sixfold symmetry 

2. R. F. S. Hearmon, Rev. Afod. Phys., XVIII (1946), 409. 

3. Voigt, op.cit., pp. 581-89; A. E. H. Love, Mathematical Theory of Elasticity (Cambridge: 
Cambridge University Press, 1934), pp. 151-61; W. A. Wooster, Crystal Physics (Cambridge: 
Cambridge University Press, 1938), pp. 241-45. 
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The conditions satisfied by the elastic coefficients of a crystal having 
any type of symmetry may be obtained from this table. As an example, 
we shall consider a crystal with cubic symmetry. Here each of the three 
principal axes has a fourfold symmetry. Since the z-axis is no longer pre¬ 
ferred, we obtain the following matrix for the elastic coefficients: 


CRYSTAL OF CUBIC SYMMETRY 


1 

1 

! C\\ 

C\t 


0 

0 

I 

1 

0 


Cu 

Cii 

0 

0 

0 

\ C\t 

Cl2 

Cu 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

1 0 

0 

0 

0 

Cm 

0 

; 0 

1 

0 

0 

0 

0 

^44 


Because of the dominant role of the cubic system in metals, it is desirable 
to have the relations between the elastic constants and coefficients in this 
system. It may readily be seen that the equations for the c’s in terms of 
the 5 ’s are 

_ _ ^11 H~ ^12 _ 

(~ ^12) ( ^11 "h 2 512) 


C12 — 


“ 'S12 


(5n “ ^12) (^11 "h 2 5j2) ' 




_ 1 _ 

•^44 


From these equations the 5’s may be found in terms of the c’s 

C\\ + C12 


— 


(Cji “ C12) (cii + 2 C12) ' 



~ C12 _ 

^^12) (fit “h 2 C12) ’ 


1 

f 44 


These relations may be expressed more succinctly in terms of the two 
shear constants, 544, 2(511 — 512), and the compressibility, 3(5ii + 2512). 
These relations are 



1 

C11+ 2Ci2 


( 22 ) 


■Su + 2 5 i 2 = 


II 


LOW-TEMPERATURE ELASTIC CONSTANTS 

IN CUBIC METALS 


A. REVIEW OF DATA 


THE elastic constants and coefficients of those cubic metals which have 
been studied are given in Table 1 . A comparison of the values from dif¬ 
ferent sources has been given by Hearmon.^ 

The coefficient C44 may be interpreted directly in terms of structure. It 
is a measure of the resistance to deformation with respect to a shearing 
stress applied across the ( 100 ) plane in the [ 010 ] direction. In a simple 
cubic structure of hard balls this resistance would be zero. On the other 
hand, the coefficients Cn and cn have no such direct physical interpreta¬ 
tion. It is therefore desirable to find linear combinations of these coeffi¬ 
cients which have simple interpretations. Such linear combinations are 
(c\\ -|- 2ci 2)/3 and (cu — C\ 2 )/ 2 . The first of these combinations may 
readily be seen to be the bulk modulus and is therefore a measure of the 
resistance to deformation with respect to a hydrostatic pressure. The sec¬ 
ond of these coefficients may readily be seen, from equations (20) and (22), 
to be the resistance to deformation by a shear stress applied across the 

( 110 ) plane in the [llO] direction. This resistance would be zero for a 
body-centered cubic lattice of hard balls. In a discussion of the physical 
interpretation of the elastic constants, attention should therefore be fo¬ 
cused upon the following three coefficients: 



Cud- 2^12 

3 


C = C44 



These coefficients are given in Table 3 . 

The relative magnitudes of Ky C, and C are frequently of interest and 
may be expressed in terms of dimensionless ratios, which are of interest 
per se. In an elastically isotropic body, C and C' must be identical. We 


shall therefore define 





as the anisotropy factor. This ratio is likewise given in Table 3 . 

When a tensile stress, 5 , is applied along the x-axis, the resolved shear 
stress across the {100} planes in the corresponding <010> directions is 

I. R. F. S. Hearmon, Rev. Mod. Phys.y Vol. XVIII (1946), Table \'. 

16 
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zero. The resulting deformation is therefore specified completely by the two 
constants, K and C'. In order to see clearly the role that these two con¬ 
stants play, we consider the stress system as composed of two parts, 
namely, 


X.= Y. = Z, = J 



TABLE 1 


Elastic Data for Cubic Metals (Room Temperature) 


CSVSTAL 

TS'PE 

Metal 

Constants in 10”“ CuVD\-ne 

COEFPS 

. IN 10‘* Dyne/Cu* 

Rep. 

1 *" 


Su 

1 

^11 

Cli 

Cu 

b.c.c . 

Fe(a) 

0.757 

- 0.282 

0.862 

2.37 

1.41 

1.16 

1 

b.c.c . 

Na(2!0°K.) 

53.5 

-23.2 

20.37 

0.0555 

0.0425 

wSSm 

3 

b.c.c . 

K 

82.3 

-37.0 

38.0 

0.0459 

0.0372 

mm 

3 

her 

W 

0 257 

- 0 073 


5 01 

1 98 

1 51 

1 

b.c.c . 

3-brass 

3.88 

- 1.52 

V V V 

0.578 

0.52 

0.275 

1.73 

X 

4 

f.C.C . 1 

A1 

1.59 ; 

- 0.58 

3.52 

1.08 

0.622 

0.284 

1 

f r r 

Au 

2 33 

- 1 07 

2 38 

1 86 

1 57 

0 420 

1 

f.C.C . 

Ag 

2-32 

- 0.993 

2-29 

1.20 

0.897 

0.436 

1 

f.C.C . 

Cu 

149 

- 0.625 

1.33 

1.70 

1 23 

0.753 

1 

f.C.C . 

Pb 

9.30 

- 4.26 

6.94 

0.483 

0.409 

1 0.144 

2 

f.C.C . 

a-brass 

1.94 

- 0.835 

1.39 

1.47 

Ml 

' 0.72 

1 

f.C.C. 

CujAu 

1.34 

- 0.565 

1.51 

2.25 

1.73 

0.663 

5 


C(diamond) 

0.14 

- 0.043 

0.23 

9.2 

3.9 

4.3 

6 


NOTES FOR TABLE 1 

1. E. Schmid and W. Boas, Kristallplastizitit (Berlin: Springer, 1935), pp. 21, 200. 

2. E. Goens and J. Weerts, Phys. Zeitschr., XXXVII (1936), 321. 

3. S. L. Quimby and S. Siegel, Phys. Rev., LIV (1938), 293. 

4. W. A. Good, Phys. Rev., LX (1941), 605. 

5. S. Siegel, Phys. Rev., LVII (1940), 537. 

6. S. Bhagavantam and J. Bhimasenachar, Nature, CLIV (1944), 546. 


The first stress system results in the dilation specified by 


_ _ _ 5 

■ 


The second stress system results in the dilationless deformation, 



The first stress system gives rise to a transverse expansion, the second to 
a transverse contraction. The ratio of the net transverse contraction to 
the longitudinal extension is therefore a measure of the relative value of C 
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and of K. This ratio is called “Poisson’s ratio” and is denoted by ctioo. From 
the above sets of deformations, it is readily seen that 

3K-2C' 

6K-\-2C'' 

Values of this ratio are given in Table 3. Poisson’s ratio is positive or nega¬ 
tive according to whether K is greater or less than (2/ 3)C'. The former is 
the case in all metals so far examined. A mineral has been studied in which 
the latter is the case, pyrite.- The elastic constants and the K, C, and C' of 
pyrite are given in Table 2. 

When a tensile stress, S, is applied along the [111] axis, the resolved 
shear stress across the 1101) planes along the corresponding < 101 > direc¬ 
tions is zero. The resulting deformation must therefore be determined 


TABLE 2 

Example of Negative Poisson’s Ratio 



1 

Constants XIO*-* 

Corrrs. X 10 

-15 

Mineral 

FeSs. 

Ul 

0.289 

Sit 

0.044 

1 

su 

0.948 

K 

0.89 

C 

1.05 

c 

2.04 


* E. Schmid and W. Boas, KrisicUphjIiziW (Berlin: Springer, 19J5), p. 267. 


solely by the two coefficients K and C. From equation (19) the tensile 
strain parallel to the [111] axis is 



Since the total dilation is (l/3Ar)5, the strain transverse to the [111] axis 
is given by 


From these two equations one finds, as the ratio of the transverse contrac¬ 
tion to the tensile extension, 

SK-2C ,,,, 

The Poisson ratio, cm, may therefore be regarded as a measure of the 
relative values of K and C. This ratio is given in Table 3. According to 

2. The Poisson ratio of sodium chlorate was reported by \''oigt as being negative but has 
recently been found to be positive (see W. P. Mason, Phys. Rev., LXX (194<^, 529; and S- 
Bhagavantam and D. Suryanarayan, Phys. Rev., LXXI [1947], 553). 
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whether K is greater or less than (2/3)C, Poisson’s ratio, <riu, is positive 
or negative. From Table 3 it may be seen that /3-brass is the only crystal 
as yet examined for which this ratio is negative, i.e., for which a specimen 
subjected to a tensile stress along a < 111 > axis extends in the transverse 
direction. 


B. RELATION OF ELASTIC CONSTANTS TO INTER.-\TOMIC 

FORCES AND TO STRUCTURE 

In his classic analysis of the elasticity of cry'stals,^ Cauchy deduced cer¬ 
tain relations between the coefficients of elasticity. In this proof Cauchy 
assumed that all forces were central, i.e., acted along lines joining the cen- 

T.ABLE 3 


El.astic Behavior of Cubic Metals 
(Derived from Table 1) 


Meial 

Kx io-‘5 

cxio->* 

C'Xl0-»* 

A 


^lll 

Fe(a). 

1.73 

1.16 

0.480 

2.4 

0.37 

-f-0.24 

Na. 

0.0468 

0.0491 

0.0065 

7.5 

.43 

+ .11 

K. 

0.040 

0.0263 

0.0044 

6.3 

.45 

-f .23 

W 

3.00 

1.51 

1.51 

1.00 

.28 

-1- .29 

^'brass . 

0.36 

1.73 

0.093 

18.7 

.39 

- .39 

A1. 

0.77 

0.284 

0.230 

1.23 

.36 

+ .34 

Au. 

1.67 

0.420 

0.147 

3.9 

.46 

: + .39 

Ag. 

1.00 

0.436 

0.151 

2.9 

.46 

; + .23 

Cu. 

1.39 

0.753 

0.237 

3.3 

.42 

+ .27 

Pb. 

0.43 

0.144 

0.037 

3.9 

.46 

i -f .35 

a-brass. 

1.23 

0.72 

0.180 

4.0 

.43 

+ .26 

CuiAu. 

1.90 

0.663 

0.180 1 

3.7 

-34 

-j- .45 

C(diamond).. 

5.7 

4 3 

2.6 ! 

1 

1.6 

0.30 

-t-0.20 


ters of the atoms. Cauchy’s proof was not entirely satisfactory, and many 
alternative derivations have been presented, of which the best known is 
due to Voigt.The Cauchy relations, first explicitly stated by Voigt, are 

~ <^23 > ^55 = ^31 > ^66 ~ <^12 » 

Cm = C\i , Cc 4 = C25 j Cis = Cm • 

In the case of a cr>'stal with cubic symmetry, these relations reduce to 

^ 44 = ^^ 12 * 

Voigt’s derivation has recently been shown by Epstein^ to be not generally 
valid. From Epstein’s analysis one may deduce,® however, as had already 

3. A.L.Cauchy, “Delapressionde tension dansun systJme de points mat^riels " Exercicts 

de mathtmatique (1828). ’ 

4. W. Voigt, Lehrbuch der Kristallphysik (Berlin; Teubner, 1910), p. 608. 

5. P. S. Epstein, Phys. Rev., LXX (1946), 915. 

6. C. Zener, Phys. Rev., LXXI (1947), 323. 
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been shown by Love/ that the Cauchy relations are valid, provided that 
one additional assumption is made regarding the structure and provided 
that the ciy^stal is initially under no stress. This additional assumption is 
that every atom is at a center of symmetry. Only when each atom is at 
a center of symmetry is the displacement of each atom completely speci¬ 
fied by the macroscopic strain. This symmetry condition is satisfied by 
many ionic and metal crystals, namely, by simple cubic, body-centered 
cubic (b.c.c.), and face-centered cubic (f.c.c.) structures. 

In ionic lattices the dominant forces arise from the electrostatic inter¬ 
action of the ions. It is therefore to be expected that in those ionic crystals 
of appropriate symmetry the Cauchy relations will be approximately satis¬ 
fied. From Table 4 this is seen to be, in fact, the case. 

TABLE 4 

Test of Cauchy Relations for Ionic Lattices* 

(280* C.) 


1 

Elastic CoErnciENTS in Units 


or 10“ Dyxes/Cm> 


<11 

<44 

NaCl. 

0.130 

0.128 

KCl. 

0.060 

0.063 


• M. A. Duraad, Phys. L (1936), 449. 


Two dominant factors in metallic binding® are the distortion of the 
wave functions of the outermost electrons, giving rise to attraction, and 
the kinetic energy of these electrons, which increases with decreasing 
volume, giving rise to repulsion. Neither of these two types of energy may 
be regarded as resulting from forces acting along lines joining atoms, i.e., 
from central forces. One therefore anticipates that in metals the Cauchy 
relations will not be even approximately valid. This is indeed the case, as 
may be seen by reference to Table 1 (p. 17). 

An understanding of the characteristically metallic contributions to the 
elastic coefficients may be best gained through a study of sodium, the 
metal whose properties conform most closely to those predicted from 
quantum-mechanical computations based upon simple assumptions. If 
we neglect the overlapping of the inner cores and the van der Waals in- 

7. A. E. H. Love, Mathematical Theory of Elasticity (Cambridge: Cambridge University 
Press, 1934), Appen. B. 

8. For discussions of metallic binding see F. Seitz, Modern Theory of Solids (New York: 
McGraw-Hill Book Co., 1940), chap, x; and N. Mott and H. Jones, Properties of Metals and 
Alloys (London: Oxford University Press, 1936), chap. tv. 
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teraction between the inner cores, an approximation which introduces 
only a slight error in the case of the alkali metals, we obtain for the energy 
per atom® 

E = 4>{v)-X. (24) 

In this equation 0 is a function only of the atomic volume, v, and is inde¬ 
pendent of shear distortions, while the function X is primarily a function 
of shear distortions, being identically zero in the absence of such distor¬ 
tion. The binding energy and the bulk modulus, K, are thus completely 
determined by the function <^>(i)). This function may be found only through 
a detailed computation of the variation of the wave function with the lat¬ 
tice-spacing. Such a computation would have to be repeated for each metal 
studied. On the other hand, the shear coefficients are completely deter¬ 
mined by X alone. This function arises from the quadrupole-quadrupole 

TABLE 5 

CONTRIBUTIO.N OF X TO SHEAR COEFFICIENTS OF 

Monovalent Metals 
(I n Units of 10‘* Ergs/Cm’) 



f.c.c. 

b.c.c. 

C«4*. 

(cu—C i2)/2. 

40.5o-« 

4.52a-* 

16.8fl-« 

2.26fl-« 


* a: Lattice constant in units of cm. 


electrostatic interaction associated with the distortions of the individual 
atoms. In all monovalent metals of the same lattice type, X is the same 
function of shear strains and of lattice parameters. It thus need be com¬ 
puted only once for all monovalent metals of the same lattice type. Fuchs 
has performed these computations and has obtained the results given in 
Table 5. 

A comparison of the computed and experimental elastic coefficients for 
sodium are given in Table 6. The close agreement of theoretical and experi¬ 
mental values of the elastic coefficients is evidence of the essential correct¬ 
ness of the assumptions upon which the computations were based. 

In metals other than the alkalis the noncoulombic interaction of the 
ions gives an important contribution to the elastic coefficients. This non¬ 
coulombic interaction is of the Cauchy type, acting along a line joining 

the centers of the two interacting ions. The energy per atom due to this 
interaction is 

E='^m(r), 

K. Fuchs, Proc. Roy. Soc., London, CLI (1935), 585; CLIII (1936), 622; CLVII (1936), 
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where ir(r) is the interaction energ>' of a single pair of ions at a distance r 
apart and the summation is over all the neighbors of a particular reference 
atom. The contribution to the elastic coefficients, K, (cu — ci2)/2, and 
C 44 , will then be Nd-E;'ds-, where N is the number of atoms per unit volume 
and s is defined by 

5r3= r^s . 


6 {x, y, 



.r + y 
2 



t 


5 (x, y, z) = {y s, 0, 0) , 


respectively. It may readily be seen that the contribution to the two shear 
coefficients is 



TABLE 6 


Elastic Coefficients of Sodium 
(In Units of 10‘* Dynes/Cm*) 



Theoretical* 

Experiment alt 

(80* K.) 

K . 

0.062 

0.051 

Cai . 

.058 

.0.S9 

f f ^ 

0.0070 

0.0072 

^Id// . . 


• K. Fuchs, Proe. R«y- Soc., London, CLI (1935), 585; CLIII (1936), 
622: CLVII (1936), 444. 

t S. L. Quimby and S. Siegel. Phyt. Rn.. LIV (1938), 293. 


The noncoulombic ionic interaction is essentially short ranged. Only a 
slight error is therefore introduced by regarding the summation in equa¬ 
tion (25) as extending over only the nearest neighbors. We then see that, if 
the crystal type and shear are such that {dr/ds) is exactly zero for all 
closest neighbors, equation (25) becomes 



(26) 


Such is the case, for example, for a (100)[010] shear in a simple cubic 
structure, and for a (110) [llO] shearin ab.c.c. lattice. In these cases AC or 
AC' have the same sign as W\ being negative when W corresponds to a 
repulsive force. It may be noted that it is just these cases for which a 
packing of hard balls would offer no shear resistance. 

Several interesting conclusions may be deduced from equation (26). If 
the ions consist of closed shells, the interaction will certainly be repulsive 
and hence the contribution to AC or to AC' will be negative. If this nega- 
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live contribution is larger than the positive contribution due to X, the lat¬ 
tice will be mechanically unstable. One thereby deduces that metals such 
as copper would be mechanically unstable as b.c.c. structures. As will be 
pointed out in chapter iv, such instability would cause the lattice auto¬ 
matically to shear into a f.c.c. lattice. Such considerations lead one to 
exj)ect to find b.c.c. structures only when the ions contain incomplete in¬ 
ner shells. Aside from the alkalis, where the noncoulombic interaction is 
negligible, b.c.c. structures are found in V, Cr, I'e, Cb, Mo, Ta, W, Eu, 
which satisfy this requirement, and in Ba and Tl, which form exceptions. 
In alloys further exceptions are found in the /3-brass type b.c.c. phases. 
In these phases, however, the concentration range of the disordered b.c.c. 
phase rapidly narrows as the temperature is lowered, the phase disappear¬ 
ing if the critical temperature for ordering is not reached in time. It thus 
appears that these phases are almost, but not quite, mechanically un¬ 
stable. This conclusion is supported by the observation that in /3-brass, 
the only alloy of this type whose elastic coefficients have been measured, 
the (cu “ coefficients is extremely small, being only one-eighteenth 
that of Cu- The positive temperature coefficient of (cn — Ci>)/ 2 in /3-brass 
has likewise been interpreted*” in terms of equation (26). As the lattice 
expands with a rise in temperature, the negative part of this elastic coeffi¬ 
cient due to \V' decreases much faster than does the positive part, thereby 
resulting in a net increase. An interpretation is given in chapter iv of 
f.c.c. b.c.c. transitions based upon the low value of the (cn — ci 2)/2 
shear coefficient of b.c.c. structures. 

10. C. Zener, Phys. Rev., LXXI {1947), 846. 



Ill 

TEMPERATURE DEPENDENCE OF ELASTIC CONSTANTS 


THE elastic moduli of crystalline materials normally decrease with in¬ 
creasing temperature, although a few exceptions are known. At least part 
of this normal decrease is frequently attributable to some sort of relaxa¬ 
tion. An example is given in Figure 4 (after Ke).* In this case the relaxa- 



O 100 200 300 400 300 

TeHPenATuRe (*C) 

Fic. 4.—Illustration of effect of relaxation upon elastic modulus. (After K6) 


tion of shear stress by viscous slip across grain boundaries decreases 
the observed modulus by 33 per cent. The effect of relaxations upon 
the elastic moduli is discussed in detail in Part Two of this monograph. 
We are at present concerned with the effect of temperature in the ab¬ 
sence of any relaxation. In such cases the normal variation of the elastic 
coefficients with temperature is linear over a wide temperature range. 
This linear variation holds from well below the Debye characteristic 
temperatures to near the melting temperatures. Examples are illustrated 
in Figures 5 (after Quimby and Siegel)^ and 6 (after Rose,^ Durand,^ 

1 . T. S. K6, Phys. Rev., LXXI (1947), 533. 

2. S. L. Quimby and S. Siegel, Phys. Rev., LIV (1938), 293. 

3. F. Rose, Phys. Rev., XLDC (1936), 50. 

4. M. A. Durand, Phys. Rev., L (1936), 449. 
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TtMPEItftTURC ) 

Fig. 6.—Temperature variation of elastic coefficients of NaCl. (After Rose, Durand Hunter 
and Siegel.) ’ * 
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Hunter and Siegel’). The normal linear decrease of the elastic coefficients 
with temperature indicates that the elastic coefficients decrease in essen¬ 
tially a linear manner with thermal energy’. 

A discussion of the temperature dependence of the elastic coefficients 
would not be complete without an analysis of the difference between the 
coefficients measured under isothermal and adiabatic conditions. Of the 
three coefficients of a cubic crystal— 02 , and -the last is the same 
under both conditions, i.e., the application of a ( 100 )( 010 ] shear does not 
lead to a change in temperature. As we have found in previous discussions, 
here also a simplification is introduced by an analysis of (cn + 2 ci 2 )/ ^, 
(cn — Cv>), 2, and C 44 rather than of Cn, Cn, and C 44 - Since (cn — Ci 2)/2 as 
well as C 44 is identical when measured isothermally or adiabatically, as 
may readily be deduced from symmetry conditions, a discussion of the 
difference between the isothermal and the adiabatic coefficients reduces to 
an analysis of only the bulk modulus, K. 

The physical principles involved in the temperature variation of the 
elastic coefficients may most readily be understood from an analysis of a 
single particle in a potential valley'. We therefore start our analysis with 
such a simple sy’stem. 

SINCLE-P.VRTICLE MODEL 

The first example we shall study is that of a single particle in a sym¬ 
metrical anharmonic potential, as illustrated in Figure 7. Zero motion will 
correspond to 0° K., vibrational energy' corresponding to thermal energy'. 
At 0° K. the force constant, C, is given by 

F"( 0 ). 

Let us now consider that the particle has thermal energy, i.e., vibrational 
energy, which, when no external force is applied, has the value Eq. Let us 
now consider that an external force,/, is applied, the time during which 
this force builds up to its final value being very long compared with the 
period of vibration of the particle. As is well known, the vibrational energy 
no longer remains constant under these conditions but follows the 

equation 

£ (x) (.r) , 

where x is the mean position of the particle and oj is the angular frequency. 
This equation follows directly upon substitution of the approximate 

solution, 

X it) = A (/) cos /‘w (/) 


5. L. Hunter and S. Siegel, Phys. Rev., LXI (1942), 84. 
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[•'ic. 7.—Illustration of the influence of thermal energy upon adiabatic elastic coeflicients 
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into the differential equation for X(0> neglecting The total energy, E, 
may therefore be expressed as a function of x as 

£(.r) (2 7) 

0)0 

where Eo and coq refer to the vibrational energy and angular frequency 
when no force is present, respectively. Since w varies as the above 

equation becomes 

E{x) =V{x) . (28) 

A typical example of the influence of vibrational energy upon E{x) is 
given in Figure 7. In this example we have taken 

V (.r) = Fo^l — cos-^^^ (29) 

and 

£o= Fo. 

In this example the thermal energy reduces the adiabatic elastic constant 
£'^(0) by 50 per cent; reduces the range of stability, i.e., the range in 
which E"(x) is positive, by 60 per cent; and reduces the maximum elastic 
energy by a factor of nearly O.L.In the more general case, in which F(^) is 
given by equation (29) and £o has an arbitrary value, one finds for the 
adiabatic coefficient 

c= K"( 0 ) (30) 

The symmetrical potential of Figure 7 might be regarded as referring 

to the potential associated with a (100)[010] or a (110)[ri0] shear, i.e., 
with C 44 or (cii — Ci 2 )/ 2 . On the other hand, a potential associated with 
dilation, i.e., with the bulk modulus A*, must be unsymmetrical. If the 
origin of coordinates is chosen to be at the minimum of the potential, 
then F'(0) will be zero, but the higher odd derivatives, such as V"'{0), 
will, in general, not vanish. The adiabatic coefficient for an unsymmetrical 
potential will be given by the second derivative of E(x) in equation (28), 

X having that value acquired from thermal expansion. It may be obtained 
either by forming the average. 
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or by minimizing the free energ)% 

F= r(.v) 


kTlnkT 

hv{x) • 


(31) 


Upon neglecting terms containing kT to a power higher than the first, wc 
obtain from these two equations 


and 


1 K'" 

2 



(3 2) 



kT 


(33) 


respectively, the derivatives being taken at x = 0. These two expressions 
are, of course, equivalent. 

The adiabatic elastic coefficient, Cs, is obtained by taking the second 
derivatives of £(x) in equation (27). We obtain 

1 / 1 / " 

kTV- kT. (34) 


On the other hand, the isothermal elastic coefficient, Cj, is obtained by 
taking the second derivative of the free energy, as given by equation (31). 


We find 




B. REAL CRYSTALS 

It is customary to assume that the normal modes of vibration of a crystal 
interfere with one another only through thermal expansion. To this ap¬ 
proximation one may write, for temperatures above the characteristic 
temperatures, 

£= (3 5) 

i 



V-kT^ln 

i 


hvj ’ 


(3 6) 


In these equations the summation extends over all normal modes of vibra¬ 
tion, and the potential energy, T, and the frequencies are functions of the 
macroscopic distortions but are independent of the w/s and of T except 
implicitly through the dilation 6. Although the potential T is a sym¬ 
metrical function of the shear strains (100)[010] and (110)[110], nonethe¬ 
less a rise in temperature will influence the elastic coefficients and 
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(^■ii ~ C 12 ) 2 both through the influence of thermal energy per se and 
through thermal expansion. Thus, if s denotes the (100)[010] shear strain, 
the effect of thermal expansion will come about becau.se of the nonvanish¬ 
ing of 

The thermal dilation. $, is found by taking that value of d which mini¬ 
mizes the free energy. F. Upon letting p denote the geometrical mean of all 
the frequencies, we obtain from equation 36 


^ SNkT dlnv 

— ir-^- 


(37) 


Here N is the number of atoms per unit volume, and K is the bulk modu¬ 
lus. The derivative — dhip dS has been given the symbol 7 by Gruneissen 

dlitp 


Since v is defined as the geometrical mean of all the frequencies, 7 depends 
upon the volume variation of all the elastic coefficients, e.g., upon A', 
C 4 A, and (cu — Ci 2 )/ 2 . Now in b.c.c. lattices the volume variation of 
(cii ~ c\2)/2 is very small. In such lattices (cn — C\2)/2 consists of two 
terms." One term, arising from the electrostatic interaction of con¬ 
duction electrons and ions, varies as 1 , a*, where a is the lattice constant. 
The contribution of the second term, arising from the exchange interac¬ 
tion of adjacent ions, increases with increasing volume. Thus a frequency 
Pi which involves solely the elastic coefficient (cn — Ci 2)/2 in a body-cen¬ 
tered metal will vaiy' as 

dlnpi ^ , 

de 

As in the case of /3-brass, such frequencies may actually increase with in¬ 
creasing volume.® One therefore anticipates that 7 will be smaller for 
b.c.c. metals than for f.c.c. metals. Druyvesteyn'’ has already pointed out 
that such is the case, with the single exception of Ca. A comparison of the 
two types of cubic metals is given in Table 7. 

If we now let s refer to either of the shear distortions—(100) [010] or 
(li 0 )[Tl 0 ]—'or to a dilation, we obtain the corresponding adiabatic modu¬ 
lus by twice differentiating E in equation (35) with respect to 5, obtaining 

C,= V"+3NK-'y^^_kT + ’'-^kT. (38) 

6. E. GrUneissen, Attn. d. Phys., XXXIX (1912), 257. 

7. N. Mott and H. Jones, Properti^ of Metals and Alloys (London: Oxford University 
Press, 1936), p. 148. 

8. C. Zener, Phys. Rev.. LXXI (1947), 846. 

9. M. J. Druyvesteyn, Philips Res. Repts., I (1946), 77. 
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Here the first term refers to the elastic modulus in the absence of any 
vibrational energy. The second term is the change in the modulus intro¬ 
duced by thermal expansion. Finally, the third term is the direct effect of 
thermal energy, discussed at length above. 

As already mentioned, the second and third terms in equation (.^8) are 
due, directly or indirectly, to thermal vibrational energy’. On the other 
hand, special effects may cause the first term likewise to vary with tem- 


TABLE 7 


CO.MPARISON OF GrONEISSEN’S CON’STAN'T FOR 
b.c.c. AND FOR f.c.c. Metals 
(.After Druyvcstcyn*) 


f.c.c. Mctni 

1 

1 

y 

b.c.c. Mclal 

y 

Au 

3.03 1 

, Ta . 

' 1.75 

Pb . 

2.73 1 

1 Mo . 

1.66 

Ir. 

2,54 1 

: \v 

1 1.62 

Pt. i 

2.54 

Fe 

i 1.60 


2.40 

Nb. 

1.59 

Pd 

2.23 i 

K. 

1.34 

Rh.' 

2,20 

V 

1.30 

A1 . 

2 17 ' 

Na. . . 

1.25 

Cu. 

1.96 

Li. . 

1.17 

Ni . 

1.88 

Cr 

0 94 

Ca . 

1.14 ' 

' Ba 

0 86 


• M. J, Druyvesleyn, Philip^ Res. Repts., I (1946), 77. 


perature. Thus in /3-brass the disordering effect of temperature is associ¬ 
ated with a marked decrease in all the elastic coefficients.*” 

For a detailed account of the theory of the temperature dependence of 
the elastic constants and of the possible relation of this variation to the 
phenomenon of melting, the reader is referred to the work of Herzfeld and 
Goeppert-Mayer,” and Born and his collaborators.*^ 

10. W. A. Good, Phys. Rev., LX (1941), 605. 

11. K. Herzfeld and M. Goeppert-Mayer, Phys. Rev., XL\’I (1934), 995. 

12. M. Born etal.,Jour. Chem. Phys., YU (1939), 59\,Proc. Cambridge Phil. Soc..XXX\'l 
(1940), 160,173,454,466; XXXVII (1941),34.177;XXXVHI (1942), 61. 67, and 82; XXXIX 
(1943), 101, 104, 113; XL (1944), 151; R. Furth, Proc. Roy. Soc., London, CLXXXIII 
(1944), 87. 






IV 

MICROELASTICITY 


THE term “elasticity” as commonly used relates to the elastic response of 
a macrospecimen to applied stresses. In macrospecimens the strain never 
exceeds a value appreciably higher than 1 per cent before either plastic 
deformation or fracture ensues. The subject of elasticity is therefore con¬ 
fined to small strains; and, because of this restriction, Hooke’s law is gen¬ 
erally applicable, namely, stress and strain are proportional. On the other 
hand, it is commonly assumed that crystallographically perfect elemen¬ 
tary regions, of, say, a hundred unit cells to an edge, will not suffer plastic 
deformation or fracture until the stress exceeds the theoretical limit. In 
fact, plastic deformation will occur simply by the elementary region’s 
passing from one configuration of minimum energy to another by a homo¬ 
geneous strain. The term “microelasticity” will refer to the relation be¬ 
tween energy and strain of such an elementary region. 

The microelastic strain energy is a periodic function of the shear strain. 
As an example we shall consider the case of simple shear across the (111) 
plane of a face-centered cubic lattice. If the shear is in the [211] direction 
and is of sufficient magnitude, a twinned structure will result, as illustrated 
in Figure 8. In this figure the atoms in one (111) plane are denoted by 
crosses. The atoms in the second plane are denoted by circles. During the 
twinning deformation these atoms move from the positions denoted by 
open circles to the positions denoted by filled circles. Similarly, the atoms 
in the third plane move with respect to the atoms in the second plane, as 
these move with respect to those in the first plane, et cetera. During this 
twinning the elementary region passes from one minimum configuration 
to another. The essential periodicity of the strain energy with respect to 
such deformation is illustrated in Figure 9. In this fig^ure the contour lines 
are given for the energy of a microelement at 0° K. The coordinates may 
be regarded as the coordinates of a particular atom in the second plane, 
or perhaps simply as the two components of shear strain across the (HI) 
plane. 

In Figure 9 the .4-type positions correspond to a f.c.c. lattice and are 
the positions of minimum energy. Twinning corresponds to passing from 
one A position to an adjacent A position. Continued twinning, alternately 

in the [211] and [121] directions, will lead to a continuing deformation m 
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the [no] direction, as illustrated in Figure 10. Thus a continued plastic 
deformation may be represented in Figure 9 by a zigzagging along the 
valleys between the maxima of type C. Along this path subsidiary' minima 
may be present in positions of type B. These subsidiary minima are sepa¬ 
rated from the primary minima by divides in positions of type D. From 
the above description of plastic deformation it is evident that the resist¬ 
ance of an unhardened metal to plastic deformation is determined by the 
variation of the energy in going from one A position over a divide D, 
through a subsidiary minimum B, and over a second divide into an adja¬ 
cent A position. 



Fig. 8.—Atomic motions during twinning in a f.c.c. lattice. Plane of atoms is (111) direc 
tion of motion is [2111. 


We shall not further discuss the interesting subject of the relation of 
microelasticity to plastic deformation but shall restrict ourselves to a dis¬ 
cussion of the variation of the microelastic energy with shear strain. 

A microelastic energy plot may be constructed as in Figure 9 according 
to two schemes. According to one scheme, one takes as independent vari¬ 
ables all the six strain components. A diagram such as Figure 9 then corre¬ 
sponds to all the strain components held constant except the two shear 
strains across the (111) plane. According to the second scheme, the only 
strain components taken as independent variables are the two shear 
strains across the (111) plane, while all the stress components are held 
constant except the two corresponding to these two shear strain com¬ 
ponents. Thus, in passing from configuration A to 5, the distance be- 
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tween the (111) planes will remain invariant according to the first scheme 
but will change so as to maintain a minimum energy according to the 
second scheme. It is believed that this second scheme of representation is 
more useful, and it will therefore be adopted in the following discussion. 





+ + + 

D'°] , 

(III) plane 

Fig. 10.—Illustration of how an unlimited amount of shear strain may occur in a f.c.c. 
structure through a succession of twin deformations. 

Suppose that a shear twin proceeds only halfway, so that the atoms in 
the second plane of Figure 8 lie midway between the open and the filled 
circles. The resulting configuration is very close to that of a b.c.c. lattice, 
with its (110) plane and [iTOJ direction coinciding with the (111) plane 
and [211] direction of the original f.c.c. lattice, as may be seen by reference 
to Figure 11. If the b.c.c. lattice is mechanically stable, the lattice will au' 
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tomatically acquire the b.cx. structure, so that configuration B in Fig¬ 
ure 9 corresponds to a b.c.c. lattice. The free energy as a function of the 
twinning shear coordinate may therefore be represented as in Figure 12. 



Fig. U.—Illustration of readjustment which causes a f.c.c. lattice to pass through a h.c.c 
structure during a twin shear strain. 



Sheor stroin 


(III) [211] sheor sir gin of f. c. c. lor ties 

(110) [lio] shear stroin of b.c.c. lattice 
Fic. 12.—Illustration of microelastic strain cnerg>' 

The inner contour lines around configurations of type A in Figure 9 are 
circles, corresponding to the fact that the elastic constant associated with 
a shear across a (111) plane is independent of the direction of this shear. 
Its value is given by 

(3 9) 


^(in) — i ^44 + I • 2 (511 — ^ 12 ) . 
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On the other hand, the inner contours around configurations of type B are 
ellipses, corresponding to the fact that the shear elastic constant across a 
(110) plane is dependent upon the direction of the shear. Thus 


while 


‘(iio)(iioJ ^ (•^n •^ 12 ^ » 

(40) 

■^{110)1001) “ ^44 * 

(41) 


In metals the first constant is always considerably larger than the 
second, the ratio of the two being as high as 18 in /3-brass. The inner 



Fig. 13.— F—s curve for copper, illustrating instability of b.c.c, structure 

contours are therefore elongated along the axis joining the adjacent 
.4-type positions. 

From equations (39) and (40) an estimate of an energy curve such as 
that in Figure 9 may be obtained from a knowledge of the elastic constants 
of the f.c.c. and of the b.c.c. structures. An example will be given for cop¬ 
per. Here we know the elastic shear coefficients of the f.c.c. structure, 
namely, C 44 and (cn — c-i^/2 are 0.75 X 10^^ and 0.24 X 10^^ dynes/cm^, 
respectively. From the computations of Fuchs we may estimate the elastic 
coefficients for the b.c.c. structure. We obtain for {cn — CvD/l the value 
— 0.02 X 10*^ dynes/cm^. An energy curve satisfying the appropriate 
curvatures is given as Figure 13. The height of the potential barrier is 
less than RT^j where is the melting temperature. This relatively small 
height of the potential barrier separating two twin configurations suggests 
that the structure of the crystalline and of the liquid state may be essen- 
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tially the same, save that the lattice contains many local fluctuations of 
structure corresponding to twins. 

A study of the microelastic free energy plot of Figure 12 yields interest¬ 
ing conclusions regarding the possibility that a rise in temperature will 
induce a f.c.c. ^ b.c.c. phase transformation in those metals in which the 
b.c.c. structure is almost, but not quite, mechanically unstable. In such 
metals the free energy, when plotted as in Figure 12, will have a shallow 
minimum in the b.c.c. position at 0° K., a minimum which lies, however, 
above that at the f.c.c. position. As the temperature is now raised, a 
shallow minimum will be associated with an unusually large amplitude 
of vibration of the (110) [IlO] shear strain coordinate and hence with an 
unusually large entropy. An unusually large entropy implies, however, a 
rapid decrease of free energy with rise of temperature. One therefore 
anticipates that, as the temperature rises, the b.c.c. minimum will be 
lowered more rapidly than will the f.c.c. minimum and may eventually 
become even lower than the f.c.c. minimum. Such considerations have 
led to a search for a transition of the b.c.c. metals to a f.c.c. structure at 
low temperatures. In this search Barrett' has discovered a f.c.c. modifica¬ 
tion of lithium at liquid-air temperatures. 

Iron suffers a double transition—b.c.c.f.c.c. at 910° C. and back 
again at 1400° C. The upper transformation may be interpreted in terms 
of the foregoing discussion, namely, in terms of the low value of the 
(cii — Ci2)/2 coefficient likely to be associated with a b.c.c. structure. The 
lower transformation finds a ready interpretation in the fact that iron in 
the b.c.c. structure, but not in the f.c.c. structure, becomes spontaneously 
magnetized and that such magnetization is associated at low temperatures 
with a lowering of the free energy. According to this viewpoint, one would 
expect the lower transformation to occur at the Curie temperature. Ac¬ 
tually, this transformation does not occur until about 150° C. above the 
Curie temperature. It is possible that this delay of 150° C. in the trans¬ 
formation arises from the lowering in the free energy associated with the 
local correlation between neighboring spins. 

1. C. S. Barrett, Phys. Rev., LXXII (1947), 245. 
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V 

FORMAL THEORY OF ANELASTICITY 


A. GENERALIZATIONS OF ELASTIC EQUATIONS 

EARLY attempts were made, originally by O. IMeyer^ and W. Voigt,^ to 
generalize the equations of the classical elasticity theory' so as to include 
anelastic phenomena. The generalizations consisted essentially in regard¬ 
ing the stress components as a linear function of both strain and strain 
rates. Thus, for the simple case of tension, the tensile stress, tr, and the 
tensile strain, e, were assumed to be related by the equation 

ff = ae + jSe . (42) 


Solids which supposedly obey equation (42) are known as “Voigt 
solids.” As illustrated in Figure 16, such solids manifest an elastic after¬ 
effect. It is evident that such solids will also manifest internal friction. 

When a stress is suddenly applied to a Voigt solid, there is no instan¬ 
taneous strain; but, as may be seen from equation (42), the strain gradu¬ 
ally approaches an asymptotic value. Conversely, when the stress is sud¬ 
denly removed, the solid suffers no instantaneous recovery, but the strain 
gradually disappears. This behavior is illustrated in Figure 16. While cork 
behaves in just such a manner,* all metals manifest an instantaneous 
strain upon sudden application of a stress and an instantaneous recovery 
upon removal of the stress. 

When a Voigt solid is subject to oscillating stresses, the rate of dissipa¬ 
tion of energy is proportional to the square of the frequency of oscillation, 
as may be seen from equation (42). Lord Kelvin^ has observed that the 
rate of dissipation of energy increases less rapidly than the square of the 
frequency. It is therefore evident that the simple generalization of the 
elastic equations represented by equation (42) does not correspond to 
reality, 

A clue as to a more successful generalization of the elastic equations 
may be obtained by representing a solid by an appropriate mechanical 


1. “ZurTheorie der innere Reibung,*Vo«r. rein. u. angru). LXXVIII (18741,130. 

2. “^er innere Reibung fester Korper, insbesondere der Metalle,” d. PAy^.,'xLVII 

0 / 1 « 

3. E. Schmid, Gesiindheits-Ing., XLVI (1923), 69. 

1890) ^27(Cambridge: Cambridge University Press, 


41 



42 


ELASTICI rV AND AXELASTICITY OF METALS 


model, as was first done by Poynting and Thomson,^ and then deriving 
the equation for this model. A mechanical model which represents the 
Voigt solid is shown in Figure 14. This model consists of a spring and 
dashpot in parallel, the latter having the property that its rate of displace¬ 
ment is proportional to the force acting upon it. A model which does not 
have the defects of the Voigt model may be obtained simply by placing a 
spring in series with the dashpot, as illustrated in Figure 15. Such a sys¬ 
tem will have an instantaneous displacement when a force is suddenly 
applied. The magnitude of this instantaneous displacement is determined 



Fig. 16.—Mechanical bchav- Fic. 17.—Mechanical behav¬ 
ior of a Voigt solid. ior of a standard linear solid. 


solely by the spring constants. As the force is maintained, the force across 
the dashpot is gradually relaxed by deformation therein, resulting in a 
gradual increase in the observed over-all deformation. Conversely, when 
the force is suddenly removed, the springs will suddenly release some of 
their stored energy, resulting in a partial instantaneous recovery. The 
complete release of all the energy in the springs must await the gradual 
relaxation of force across the dashpot. Such a system therefore manifests 
the general features of elastic after effects in real solids, illustrated in Fig¬ 
ure 17. Further, the rate of energy dissipation by such a model during 
forced oscillation does not increase rapidly as the frequency increases. In 
fact, it may be seen that for large frequencies, when the displacement 
contributed by the dashpot is small compared with that contributed by 

5. J. H. Poynting and J. J. Thomson, Properties of Matter (London: C. Griffin & Co., 1902). 
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its companion spring, the rate of energ>' dissipation is independent of 
frequency. The equation corresponding to the mechanical model of Fig¬ 
ure 15 may be seen to have the form 

Ciff a-iff = bit • (43) 

Since this is the most general linear homogeneous equation in stress, 
strain, and their first time derivatives, a solid which obeys equation (43) 
will be called a “standard linear solid." 

While an elastic body subjected only to tensile stresses has only one 
constant, namely, the tensile elastic modulus, the corresponding standard 
linear solid has three essentially independent constants. Thus in the nota¬ 
tion of equation (43) the three independent constants may be taken as 
0*2 flj, bi ai, b- 2 'ai. In order to facilitate the understanding of the behavior 
of a standard linear solid, three new independent constants—r,, r,, Mn - 
will be introduced as follows: 

ff + T,(7 =Mn (« + Trf). (44) 

Here r, is the time of relaxation of stress under conditions of constant 
strain. Thus suppose that both £ and e are zero. Then equation (44) re¬ 
duces to 

a + Tiff — 0 , 

which has the solution 

0 - (/) = 0 - ( 0 ) . ( 45 ) 

On the other hand, suppose a strain, cq, is suddenly applied at / = 0. The 
stress then relaxes with the relaxation time r, to its equilibrium value 
Mnio. Thus the solution to equation (44) corresponding to this condi¬ 
tion is 

O'(/) -Mr€o-\- (ffo — Unto) , ( 46 ) 

where o-q is the initial value of the stress. Since the final value of the ratio 
of stress to strain —the value after all relaxation has occurred—is Mr, this 
quantity is known as the “relaxed elastic modulus." The constant is 
the time of relaxation of strain under conditions of constant stress. Thus, 
if both <7 and <r are zero, t obeys the equation 

€ + = 0 . ( 47 ) 

On the other hand, suppose a stress (Tq is suddenly applied at / = 0. The 

strain then relaxes with the time of relaxation to its equilibrium value, 

Mn^ao- Thus the solution to equation (44) corresponding to this condi¬ 
tion is 

c(/) =M«Vo-{- (€o-M«Vo) 


( 48 ) 
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From the three independent constants, t^, and Mr, certain new con¬ 
stants may be derived which have considerable physical significance. 
Suppose that in a very short time 5/, the stress receives a finite increment 
A(r. If we now integrate both sides of equation (44) with respect to time 
over the time interval 5/, the first term in each side approaches zero as the 
time interval 5/ is made smaller and smaller. We are left with the follow¬ 
ing relation between the increments of stress, Ao-, and of strain, Ac: 

r«A<7 —Mrt«^€ . 


The ratio A(r/Ae will be called the “unrelaxed elastic modulus” and will 
be denoted by Mu, since it gives the relation between changes in a- and 
in € which occur so rapidly that no relaxation has time to take effect. Thus 

Ao- =Aft/Ae , (49) 


with Mu given by the equation 



If 

T. 


(50) 


The deviation of the ratio Mr!M u from unity may be used as a measure 
of the relative change in stress or in strain which may occur through re¬ 
laxation. Examples of relaxation are given in Figures 18 and 19 for the 
cases in which Mr/Mu is 0.1 and 0.9, respectively. The first case is again 
illustrated in Figure 20, with time given on a logarithmic scale. In this 
figure it is particularly evident that the relaxation time of strain at con¬ 
stant stress is larger than is the relaxation time of stress at constant strain. 

Solids are frequently investigated by dynamical methods. It is therefore 
of interest to find the relationship between stress and strain in the stand¬ 
ard linear solid when these quantities are periodic. Toward this end we 
substitute in equation (44) the solutions 

< t { 1 ) = o-oC^', e(/) =eoe^S 


and endeavor to find the relation between <jq and co. We obtain 

(1 -j-fwr*) O’© —Mr (1 +ia>r<f) eo 

cq = SD?eo i 


(51) 


or 


where the complex modulus 9)2 is given by 



1 

1 d-tcoTf 



(52) 


The angle 6 by which strain lags behind stress is of particular interest, 
since its tangent is frequently used as a measure of internal friction. Upon 
observing that 

Imaginary part of 2)2 
° Real part of 2)2 ’ 



f 
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we obtain 


tan 8 = 


u> (r<r — T,) 

1 + w- (r<rr«) 


(53) 


In order to represent the variation of internal friction with angular fre¬ 
quency in as vivid a manner as possible, we shall introduce the geometric 
mean of the two times of relaxation, and r^, 

T = (54) 



Fic. 20.—Relaxation of stress and strain for special case of Ma }fv — 0.1. Ixigarithmic 
time scale. 


and the geometric mean of the two moduli, Mu and Mr. 

M = (MuMr) (55) 


Making use of equation (50), we then transform equation (53) to 


tan 



cur 

1+ (wr)*- 


(56) 


The first factor in equation (56) is essentially the relative difference in 
the relaxed and the unrelaxed elastic moduli. The second factor gives the 
frequency variation of tan 5. This factor has a maximum when the prod¬ 
uct, WT, is unity. When this factor is plotted against the logarithm of this 
product, it is seen to be a symmetrical function of wt, with the general 
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characteristics of an error curve. The maximum vaiue of this second fac¬ 
tor is , and hence 


(tan 5) 



It is also of interest to have some real measure of the ratio of stress to 
strain, a ratio which will be given the symbol The definition of is 
not unique. Thus we could define it as the magnitude of iOJ, as the real 
part of 9[)L or as the reciprocal of the real part of air'. The latter definition 
will be adopted simply because of the symmetrical expression to which it 
leads. According to this definition, is the ratio of the stress to that 
part of the strain which is in phase with the stress. From equation (52) 
we find that this definition leads to 



1 

a . 

1 


Upon using equation (50b we find may be written as 



Mu — M a 
1+ (C0f)2- 


In the limits of low and high frec|uency reduces to the following 
limits: 


(Mr , a>r« 1 , 

Mui — ( 

\Mu ) <jff^ 1 . 

The precise manner in which varies between these two limits is shown 
in Figure 21. The variation of is seen to be a maximum when the prod¬ 
uct WT is unity, just where tan 5 is also a maximum. It may likewise 
be seen from this figure that has effectively reached its two extreme 
limits at frequencies where tan 6 is still appreciable. This effective narrow¬ 
ness of the transition range of as compared to the range in which tan 6 
is appreciable, is a manifestation of the difference between the asymptotic 
expansions of these two quantities, as illustrated in the following equa¬ 
tions: 


M u — Mw — 


tan 5 = — 


My-MR ' 
(ur) 2 

M y—M R 




^ a>T» 1 , 


Mu~MR — {My — Mr) • (tor) ^ 

, M u’-M R , 

tan 6 =--• (tor) 


tOT« 1 
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B. Boltzmann’s superposition principle 

While the standard linear solid discussed above has certain general 
features in common with actual solids, it does not reproduce precisely the 
behavior of any real metal. In an attempt to formulate a theory which 
corresponds more closely to reality, two paths may be followed. On the 
one hand, one may attempt to describe the mechanical behavior by 
means of a differential equation, similar to equation (47) for the standard 



Fig. 21.—Frequency dependence of internal friction and elastic modulus 


linear solid but containing higher derivatives of stress and of strain. This 
procedure has recently been discussed by Alfrey.® On the other hand, one 
may abandon the attempt to derive a differential equation relating stress 
and strain and develop all the conclusions derivable from the assumption 
that the true relation between stress and strain is linear. This latter meth¬ 
od was introduced by Boltzmann^ and will be discussed in detail in this 
section. 

If the fundamental relations between stress and strain are linear in 


6. T. Alfrey, “Non-homogeneous Stresses in Visco-elastic Media,” Quart. Jour. Appl- 
Math., II (19^), 113; and “Methods of Representing the Properties of Visco-elastic Mate¬ 
rials ” im., Ill (1945-46), 143. 

7. L. Boltzmann, “Zur Theorie der elastische Nachwirkung,-” Ann. d. Phys., VTI (1876), 
624; also SUztingsb. Akad. IFw. Wien. Ber., LXX (1874), 275. 



FORMAL THEORY OF .ANELASTICITY 


49 


stress and strain and in their higher time derivatives, then the solutions 
satisfy the principle of superposition. Thus, let A(0 be the deformation 
produced by the force Fi{i) acting alone and be the deformation 
produced by the force ^ 2(0 acting alone, then this principle states that 
Di{t) + DiQ) is the total deformation produced by the forces Fi(f) and 
7 ^ 2(0 acting together. E.vpressed in more formal language, the superposi¬ 
tion principle states: 

If Di (/) is produced by Fi (/) 

and Di (/) is produced by F 2 (/), 

then Di (/) -|- Do (/) is produced by Fi (/) + F2 (/) . 

Conversely, the deformation may be regarded as the independent variable, 
the applied force as the dependent variable. For example, the deformation 
may be varied according to a predetermined schedule, and the force meas¬ 
ured which is necessary to produce this deformation. In this case the su¬ 
perposition principle states: 

If Fi (/) is required for Di (t) 

and Fi (/) is required for (t ), 

then Fi (/) -{- Fi (/) is required for Di (/) + Do (/). 

The superposition principle may be formulated in a still more concise 
manner. Toward this end we define the quantity 5(/), which will be called 
the “creep function,” as the deformation resulting from the sudden appli¬ 
cation at ( = 0 of a constant force of magnitude unity. Now, if during the 
time interval from t to t dt the force changes from F to F + Fdt, we 
may consider that during this time interval a constant force has been ap¬ 
plied of magnitude Fdt. According to this vie^vpoint, the deformation at 
any instant / is the result of a continuous series of constant forces previous¬ 
ly applied. Thus, according to this viewpoint, the superposition principle 
states: 

D(t) = J\(t-l')FU')dt'. (57) 

Conversely, as before, we may regard the deformation as a specified func¬ 
tion of time. Toward this end we define the quantity J{t), which will be 
called the “stress function,” as the force which must be applied in order 
that the deformation may suddenly change at / = 0 from zero to unity 
and remain unity thereafter. By applying the same reasoning as that used 
for equation (57), we derive 

FU) - f f{t~nD{ndF. 


(58) 
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The anelastic behavior of a solid may be regarded as completely spe¬ 
cified either by b{t) or by/(/). In order to obtain a relation between these 
two quantities we shall apply equation (57) to the case in which 

ro , / < 0 , 

DU) =< 

ll , />0 , 

and hence where also 

fO , /< 1 , 

FU) = 

’/(/) , ^> 0 . 


After / = 0. the left-hand side of equation (57) is equal to unity, while 
the right-hand side may most conveniently be separated into two terms. 
The first term is the integral over the time interval in the immediate 
vicinity of ^ = 0 during which D{t) is changing from 0 to 1. The second 
term is the integral after D{t) has attained the value unity. We thus find 

1 = 5(/) /(O) -h f‘ dU-n findi'. (59) 

This equation relates the two quantities 5(/) and /{I). The precise evalua¬ 
tion of one function in terms of the other would be ver>' difficult, requiring 
iterated numerical integrations. Certain useful relations may, neverthe¬ 
less, be derived without such numerical integrations. The two factors of 
the integrand in equation (59) are plotted in Figure 22 as a function of t' 
for typical examples. We see that 8(t — t') differs appreciably from 5(/) 
only when / is verj^ small. It will therefore be advantageous to write 

5 (/-/') as 5(/) - (5(/) - 5(/-C) . 

The first term then contributes the major part of the integral, and, further¬ 
more, its integration may be carried out directly, leading to the following 
relation between the creep and stress functions: 

5 (,) /•(/) = IH- /'*{ 5{/) - 5 U-n \ findi '. (60) 

•'o 

Now the first factor in the integrand is always positive, and the second 
factor is always negative, so that the creep and stress functions satisfy 
the following inequality: 

S(i) fit) . (61) 

It may be seen from equation (60) that the equality sign is valid at verj' 
short times. If 8{t) approaches a limiting value, the equality sign is also 
valid at very long times, the inequality being appropriate only for inter¬ 
mediate times. 

.As an e.xample of the above discussion, we shall evaluate 5{t),/it), and 
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the product for the standard linear solid. From equations (49) and 

(50) we obtain 

a (l) =m 7 - -Ma^) (62) 

/ (/) =MH-\- {Mv~M h) c“' ^ (63) 

These functions and their product are plotted in Figure 23 for the special 
case of Mv, Mr — 10. The product h{t\f{t) for this particular case reaches 
a minimum value of 0.65. 




Kici. 22.—Interpretation of integrand of equation (19) 


Anelastic experiments are frequently made in which a periodic displace¬ 
ment or force is applied. We shall now examine how the results of such e.\- 
I)eriments arc related to the creep and relaxation experiments. 

In order to analyze the force required for a periodic deformation, we 
shall find it convenient to re-write equation (58) in another form. We first 
integrate by parts, obtaining 

FU) =/(0) Z)(f) + /*' /(/-/')/>(/') (64) 

and then substitute t for / — obtaining 

Fit) =/(0) />(/) + rf{r) D{t-r)dT. 

Jo 


(65) 
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From this equation we find that the force required for the periodic defor¬ 
mation, 

D {t) Z?o sin a>t , 


is given by 


F (/) = Do I / (0) d-y' / {t) cos (at dt\^ sin wt 

— Do I y f it) sin (at dt\^ cos cat. (66) 



.1 I 10 

Time (lo Units of %) 

Fig. 23.—Example of relation between creep and stress relaxation 

From equation (66) we see that a periodic displacement requires a peri¬ 
odic force. Conversely, a periodic force gives rise to a periodic displace¬ 
ment. The corresponding stress-strain diagram is therefore an ellipse, and 
the relation between force and deformation may be represented by a com¬ 
plex modulus 9D?. The real part of this complex modulus will be called the 
“dynamic modulus” and will be denoted by Af(a>). 

From equation (66) we may derive, in certain important cases, the rela¬ 
tion between the dynamic modulus and the stress function. The im¬ 
portant case for which this relation may be derived is that in which 
d Inf/d In t is small compared with unity and varies only slightly with 
In t. An example of a stress function satisfying this condition is given 
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in Figure 24. In this figure d In J/d In t has the constant value 0.20, 
and / decreases from 0.8 to 0.5 as t increases from 0.1 to 1,000. The 
salient features of such a plot of f{t) on a linear time scale are independent 
of the time units. Thus, irrespective of whether the time scale extends from 
0 to 10, from 0 to 100, or from 0 to 1,000 seconds, / drops 30 or 40 per 



Fig. 24.—Illustration of a typical stress relaxation plotted on a linear scale 


cent during the first tenth of the time interval and 8 per cent, at most, 
during the remaining nine-tenths of the interval. From equation (66) the 
dynamic modulus is given by 

jW (w) = / (0) + /" / {t) cos (at dl . 

•'0 

When integrated by parts, this equation becomes 

M i<a) = J' f ^ sin X dx , 
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The contributions to the integral from the ranges O-tt, 7r-27r, 27r-37r, . . . , 
are of alternate signs and of slowly decreasing magnitude. A very good 
approximation to the sum of this alternating series is given by the inte¬ 
gral extended over the first half of the first range. Thus 


MM 



sin .V dx. 


As may be seen from Figure 24, a good approximation to this integral is 
obtained by regarding/( a:' co) as a constant, with its value in the center of 
the range over which the integration extends. We are therefore led to the 
following relation between the dynamic modulus and the stress relaxation 
function: 

.UM=/Q, (67) 


where P is the period of oscillation, Iir w. 

From equation (66) we may also derive, for the above-discussed case, 
the relation between the stress function and the internal friction. As has 
previously been mentioned, the internal friction is the tangent of the 
angle 5 by which the deformation lags behind the applied force in periodic 
oscillation, and it is, when small, proportional to the relative energy dis¬ 
sipated per cycle. The tangent of the angle by which the deformation lags 
behind the force is the ratio of the coefficient of cos ui to the coefficient 
of sin (jjt in equation (66), and hence 



If we now write the integrand of the numerator as {df/d In t) (sin 
the first factor may be regarded as essentially constant. Integration over 
the remainder of the integrand gives the numerical factor ir/2. Now the 
major part of this integral, 87 per cent, comes from the integration over 
the range 0 < wf < The essentially constant factor d//dlnt can 
therefore reasonably be assigned its value in the middle of this range. 
We thus obtain the following relation between the internal friction and 
the stress function: 


tan 5 M = — 




(69) 


Equations (60) or (61), (67), and (69) constitute three independent re¬ 
lations between the following four functions: creep, stress, dynamic 
modulus, and internal friction. The equality sign in equation (61) and 
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equations (67) and (69) are good approximations for the cases typified by 
Figure 25. Three more approximate equations could therefore be derived 


from the above, one between d/(cj) and 8{l) by substitution of for /"in 
equation (67); one between tan 5 (w) and 5 by substituting 5“' for f in 


equation (68); and, finally. 

tan 6 {.(jo) = 


TT ({ InM tojl 

2 (i III u) 


(70) 


obtained by combining equations (67) and (69). A beautiful experimental 
demonstration of the validity of equations (67)-(70) has been presented 
by Ke*^ in his study upon the viscosity of grain boundaries. 



»00 100 10 I 

•— rMOUCKV (OTOUS/MC.) 

Flo. 25.—Internal friction measured by Benncwitz and Rdtger. The calculated values of 
the maxima are marked bv vertical lines. 


C. RELAXATION SPECTRUM 

A particularly simple example of anelasticity was discussed in Sec¬ 
tion A. In this example a certain fraction of the initial stress relaxes in an 
exponential manner when the strain is maintained constant. Conversely, 
when the stress is maintained constant, the strain approaches an asymp¬ 
totic value in an exponential manner. In chapter vii we shall see that many 
examples occur in nature in which the creep and stress relaxations obey 
these exponential laws. One is therefore tempted to regard all examples of 
anelastic phenomena as the superposition of elementary processes in which 
the stress relaxes exponentially. This point of view was adopted many 
8. T. S. Kd, Phys. Rev., LXXI (1947), 533. 
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) ears ago by J. J. Thomson^ and by E. Wiechert.*® In previous discussions 
of this point of view it has been assumed that the stress which relaxes is 
only a small percentage of the original stress. The recent observations of 
W. A. West’^ show that the major fraction of the stress may, in fact, relax. 
It is therefore necessaiy- to re-examine this viewpoint. 

In discussing various types of stress relaxation, it is desirable to have 


some dimensionless measure of the total stress relaxation. Expressed in 
terms of the unrelaxed and relaxed moduli, one could use (Mu- Mu)/Mu 
or (Mu — Mu), Mr or (Mu — MR)/y/MuMu- In Section A of chapter 
vii we shall see that, in the case of simple exponential relaxations, the sec¬ 
ond of the above measures of relaxation is most appropriate, since it may 
readily be expressed in terms of the physical constants of the system. We 
shall therefore define the relaxation strength of a simple e.xponential re- 
la.'cation by the quantitv 


(71) 


In most experiments the pertinent modulus is either Young’s modulus £, 
the rigidity modulus G, or the bulk modulus K. The measured relaxation 
strength will, in general, depend upon which modulus governs the rela¬ 
tion between stress and strain. A subscript will therefore be given A to de¬ 
note the type of experiment to which the relaxation strength refers. Thus 



(72) 

(73) 

(74) 


Generally no relaxation occurs under a pure hydrostatic pressure; thus in 
these cases is exactly zero. When Ar is zero, a relation exists between 
Ae and Ag. From the relation*- between E, G, and K, namely, 


one derives 



9 KG 
3A" + G’ 



9. “On Residual Effects,” in his Applications of Dynamics to Physics and Chemistry (Lon¬ 
don; Macmillan & Co., 1888), chap. viii. 

10. “Cber elastische Nachwirkung” (inaugural dissertation [Konigsbergi. Pri., 1889]); also 
“Gesetze der elastische Nachwirkung fur konstant Temperature,”/I Phys.,h (1893), 335. 

11. “Elastic After-effects in Iron Wires from 20® to 550° C.,” Trans. CLXVTI 

(1946), 192. 

12. R. V. Southwell, Theory of Elasticity (Oxford: Clarendon Press, 1936), p. 307. 
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Upon using the following formula relating Gt, and Eu in terms of Pois¬ 
son’s ratio, (t: 

Ev = 2{\ + <T)Gvr 

we obtain 



3 

2 (1 + cr) 



(75) 


For the usual range of a, Ac is about 15 per cent larger than A^. 

The relaxation strength of a single crystal will depend upon its crystal¬ 
lographic orientation, as well as upon the type of deformation to which it 
is subjected. If we assume that a cubic crystal undergoes no relaxation 
under a pure hydrostatic pressure, it will have two relaxation strengths, 
corresponding to the relaxation of the two shear coefficients, Cu and 
{cn — Cn)/ 2 . These two relaxation strengths will be denoted by 6 and 
6', respectively. In order to compute Ag and Ac in terms of these two 
relaxation strengths, we need the following formulae for the tensile and 
torsion moduli:*^ 

“ ^11 “ I 2 (ill — 512) “ 544] ^ , { 75 a) 

G-* = 544 + 2 [2 (j„- Sn) - Sii] <t> » ( 75 ^>) 

where 

j i 2 { 2 2 i 2 2 

<p = + yzyz + 73 ^ 1 » 

72, 73 being the direction cosines of the specimen’s axis with respect 
to the principal axes of the crystal. We shall also need the following formu¬ 
lae relating the elastic coefficients (c’s) to the elastic constants ( 5 ’s); 

1 


^44 = 

C 44 ’ 

(75c) 

5u “ 512 = 

1 

{75d) 

Cll — Cit' 

■511 + 2 5i2 = 

1 

(7Se) 

C 1 I+ 2Cj2 ’ 


equations which are identical to equation (22) of Part I. From equations 
(75a)-(75c) we find directly that, for cubic crystals, 



3 


(C3'-c'a) 0 


<^44 (Cll+ C12) 
2 (cn + 2 C 12 ) 


iC~C')<f>' 


(75/) 


Art — 


C'5 + 2 (Ca'-C'5) <f> 
C'+2(C-CO0 


(75p) 


1935^ p ^2?' (Berlin: Springer, 
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where C and C denote the two shear coefficients, namely, 

C-= C44 - 

Cii C12 


In deriving equations (75/) and ilSg), we have made use of the equation 


6 = 


(Cii) a — (^44) tf (S44) R {■y^ 4 ) V 


(<■ 44 ) R 


{ ^ 44 ) u 


and of a similar equation for 6'. 

If the total relative stress relaxation of a specimen is much less than 
unity, the stress relaxations due to different sources may be regarded as 
independent of one another. A rela.xation spectrum may then be defined 
unambiguously as the plot of the distribution function, A(r), against In r, 
where A(t) - d In t is the contribution to the total relaxation strength of 
those simple exponential relaxations whose times of relaxation lie within 
the range d In r at r. The distribution function, A(r). satisfies, by defini¬ 
tion. the equation 

5 = /" A (r) d III r . (76) 

•/ — 00 

The relative stress which yet remains to be relaxed by the simple exponen¬ 
tial relaxations in the range d In r at r is 

A (t) e~‘^'d In t . 

The distribution function. A(t). may therefore, at least in principle, be 
determined from stress-relaxation experiments through the integral 
equation 

/« 7-. (77) 

%/ _m 


a {co f 


Again, the contribution to internal friction arising from the simple ex¬ 
ponential relaxations with relaxation times in the range d In r at r is, 
from equation (53), 

UlT 


A(r) 


1 -{-cu^r® 


d In 


T . 


The distribution function may therefore, in principle, be determined from 
internal-friction measurements through the equation 


tan 5(a)) = J" A(r)-j 


(jiT 


+ 0)^2 


d In T . 


(78) 


As has been mentioned above, the distribution function, A(r), may be 
computed in principle either from a knowledge of the stress-relaxation 
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function, by means of equation (77) or from a knowledge o! 

the internal friction, tan 6, as a function of frequency by means of equa¬ 
tion (78). In each case the determination of A(r) involves the inversion of 
an integral. This inversion is practicable only when the stress relaxation 
and the internal friction are known analytic functions of time or fre- 
quenc)', respectively. In such cases standard inversion methods may be 
adopted, as has been recently discussed by Gross.'^ In such inversions 
small errors in the stress relaxation or internal friction functions would be 
reflected as relatively large errors in the distribution function. A(t). 

Because of this extreme sensitivity of A(t) to small errors in observa¬ 
tion, it has not seemed to the author profitable to express experimental 
results in terms of this distribution function. Rather, it has appeared 
wiser to express experimental results directly as tan 5. A plot of tan 5 
against the logarithm of time may then be regarded as the relaxation 
spectrum, as is illustrated by the Frontispiece. 

14. B, Gross, "On Creep and Relaxation.” Jour. Appl. Phys., X\ I1! (1947). 212. 
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A. RELATION BETWEEN MEASURES OF INTERNAL FRICTION 

IN THE determination of the relaxation spectrum, namely, in the deter¬ 
mination of the internal friction, tan 6, as a function of frequency, the 
angle of phase lag S could be measured directly. This method must be 
used if the internal friction is comparable to, or greater than, unity. In the 
past only experiments have been reported in which the internal friction 
is small compared to unity. In such cases a variety of indirect methods 
may be employed to determine the internal friction. 

One indirect method of measuring the internal friction is to observe the 
amplitude of vibration under conditions of forced oscillation, the fre¬ 
quency of the impressed force being slowly varied while its amplitude re¬ 
mains constant. The amplitude of vibration will be a maximum when the 
impressed frequency is equal to a critical resonance frequency, p, of the 
specimen and decreases essentially to zero as the impressed frequency 
either increases or decreases away from this critical frequency. Let Av be 
the change in impressed frequency necessary to change the amplitude 
from half-maximum on one side of the maximum to half-maximum on the 
other side. It will then be shown below that, provided that tan 6 is a small- 
order quantity, at the resonance frequency 

tan 6 = . (79) 

V3v 


A second indirect method of measuring internal friction involves the 
measurement of the energy, AE, dissipated per cycle of vibration. When 
tan 5 is a small-order quantity, one may define the energ>' of vibration, E, 
as the strain energy where the strain is at a maximum. When the strain 
is not at a maximum, part of the energy E is in the form of kinetic 
energy. The ratio AE/E is commonly called the “specific damping capac¬ 
ity.” It is shown below that, in the case of small tan 5, 


1 AE 


( 80 ) 


A third indirect method of measuring internal friction is to observe the 
gradual decay in amplitude of vibration during free oscillation. The loga¬ 
rithmic decrement (log. dec.) is then defined as the natural logarithm of 
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the ratio of successive 
small tan 5, 


amplitudes. It is shown below that, 


tan 8 = 


_ log. dec. 

TT 


in the case of 
(81) 


It is evident that the above relation must become invalid for large values 
of the internal friction, for then the internal damping of the specimen itself 
will overdamp the vibration, resulting in aperiodic motion. The author* 
has shown that critical damping occurs when 5 is 53°, provided that the 
specimen has only a single pair of relaxation times, r, and t^. In practice 
it is often convenient to obser\'e the time, t„, required for the amplitude 
of vibration to decay to one-«th of its original value. If v is the frequency 
of vibration, it is then evident from equation (81) that 


tan 8 =-. 

TTUTn 


(82) 


We shall now derive the theoretical relations given in equations (79)- 
(81). In order to derive the first of these equations, we must write the 
equation of motion of the specimen as a whole. This is particularly simple 
when an auxiliary inertia member is used of such magnitude that the 
inertia of the specimen itself may be neglected. Let I be the inertia of the 
auxiliary member, D the deformation of the specimen, / the “force” with 
which the auxiliary member acts upon the specimen, and F the external 
“force” which acts upon the auxiliary member. Thus in the case of tor¬ 
sion, I is the moment of inertia of the auxiliary member, D is the angle of 
twist, and / and F are torques. The equation of motion of the auxiliary 
member is then 



(83) 


In the case of periodic motion we may consider D,f, and F to vary with 
time as 

= Z7oe'"', etc. ( 84 ) 

Deformation and force may then be related by 


F = mD, (85) 

where the complex modulus iOi may be written as 

902 =^(1+itan 5). 

Upon substituting equations (84)-(86) in equation (83), we obtain 





_ M _ 


C. Zener, Jour. Appl. Phys., XVIII (1947), 1022. 


(86) 
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The amplitude of the deformation is a maximum at the resonance angular 
frequency, 



This amplitude is half its maximum value when 

1 — ^ = + Vs tan 6. 

The difference of the two values of oj which satisfy this equation. Aco, is 
seen to be given by 

Au) 

— = vs tan d . 
wn 


in the case in which tan d is small. This equation is identical with equa¬ 
tion (79). In cases in which an au.xiliary inertia member is not employed, 
the proof of equation (79) follows in a similar manner, except that a dif¬ 
ferent interpretation- must be assigned to the quantities D,f, and I. In 
particular, D is the coefficient of a normalized characteristic displacement 
vector in the expansion of the displacement vector in terms of these char¬ 
acteristic vectors. Analogous interpretations are assigned to / and F, 
while I becomes the density of the material. 

We shall now derive equation (80). The rate at which energ>' is dissi¬ 
pated per unit volume is <ri. The energ>' dissipated per cycle throughout 
the specimen is therefore 



Ujion setting 

<7 = <70 cos ul , 

€ = Co COS (w/ — 5) . 

we obtain 


A£ — ir sin 6 fioffud v . 


(87) 


On the other hand, when tan 6 is small, we may take the energ>’ of vibra¬ 
tion per unit volume as Woe©- The energy of vibration of the whole speci¬ 
men is therefore 

E = 5 / (7o€o d V . ( 88 ) 


Upon comparing equations (87) and ( 88 ), we obtain 


E 


= 27r sin 


b 




which is identical with equation (80) in the case in which tan b is small. 


2. C. Zener, Phys. Rev., LII (1937), 230. 
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Equation (81) may be derived most readily directly from equation (80). 
In free oscillation the logarithmic decrement for the energ>' of vibration 
is simply AE/E. Since both AE and E are proportional to the square of 
the amplitude of vibration, provided that this amplitude is small, the 
logarithmic decrement for amplitude (log. dec.) is just one-half that for 
energy; whence 


1 AF 

log. dec. = - ~ 


B. .ME.\SUREMENT OF INTERNAL FRICTION 

The earliest measurements^ of internal friction were made by the 
method of free decay. The specimen formed the suspension element of a 
ballistic galvanometer, and the logarithmic decrement was obtained by 
observing the rate of decay of the amplitude of vibration. This early meth¬ 
od has recently been revised and used with great success by Kc.^ When 
the specimen is in the form of a wire, the period of oscillation may readily 
be made of the order of a second, and so the obsei^^ations may be made 
visually. When specimens of the general size of a standard tensile si>eci- 
men are used, the periods of vibration are correspondingly shorter, so that 
photographic or other nonvisual means of recording must be adopted. 

Such equipment was first used by FdppP and has recently been improved 
by Norton.® 

The free oscillation method may also be employed without the use of an 
auxiliary inertia member. In this case the specimen is excited in one of its 
normal modes of vibration, and the decay of amplitude is then observed. 
In such cases the frequency is considerably higher than when an auxiliary 
inertia member is employed, and hence direct visual observation may no 
longer be employed. For merely relative values of the internal friction, 
the time of decay may be determined remarkably accurately by the ear,' 
provided, of course, that the frequency is in the audible range. This audio 
method has been used with considerable success by Waller.' Electrical 
methods are commonly employed accurately to observe the decay of 
oscillation. Thus by measuring the time t„ required for the amplitude to 
be reduced to one-«th of its original value, the internal friction may be 
computed from equation (82). On the other hand, electrical methods have 

3. W. Weber, Poggendorfs Ann., XXXIV (1837), 247. 

(1948)^281*^*’ XIX 

5. O.Foppl, Ber.d.H'erkstofansschussesd. Ver. denlsc/,. PJsmhultenUnlf.Ko. 36 1923 

o. J. Norton, Rev. Set. Instruments, X (1939), 77. 

7. Mary D. Waller, Proc. Phys. Soc., L (1938), 144. 
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been described® in which one measures directly the number of oscillations 
required for the amplitude to be reduced by a predetermined factor. 

Rather closely related to the method of free oscillation is the method 
in which a high-frequency pulse is introduced at one end of a bar, and the 
amplitude of the pulse is observed as it suffers successive reflections at the 
ends. This method has been employed by Huntington.^ On the other 
hand, this method suffers from the danger that the loss of intensity of 
the wave packet may be due largely to scattering rather than to true ab¬ 
sorption, as, for e.xample, scattering by the random orientation of crystal¬ 
lites observed by Mason. 

When resonance frequencies are employed, the time of decay in free 
oscillation is frequently too short to measure conveniently. In such cases 
it is customary to measure the width at half-maximum of a resonance 
curve. The internal friction is then computed by equation (79). The reso¬ 
nance and decay methods may be regarded as complementary to each 
other. The former method is more applicable to those cases in which the 
width at half-maximum is comparatively large and becomes difficult as 
the internal friction is reduced so low that the half-width is too small. On 
the other hand, the decay method is more applicable when the time of 
half-decay is comparatively long. It becomes difficult when the internal 
friction increases to a value such that the time of half-decay is too short 
to measure accurately. 

The free decay and resonance methods of measuring internal friction 
have been especially popular, since they require only relative, rather than 
absolute, measurements of amplitude of vibration. Various electrical and 
magnetic methods have been employed in these measurements. Wegel and 
Walther" introduced the electromagnetic method of induced eddy cur¬ 
rents. In this method the driving end of the specimen is subjected to both 
a static and an osciUating magnetic field. The osciUating field induces eddy 
currents in the specimen, which, in turn, are pulled and pushed by the 
stationary, inhomogeneous, magnetic field, thereby exerting an oscillating 
force upon the specimen. At the detecting end the specimen is subjected 
likewise to a stationary magnetic field. Vibrations therefore induce eddy 
currents at this end. The field of these eddy currents is then 
detected by coils surrounding this end. This electromagnetic drive 
and detection through eddy currents has been used by the author in 

8. F. Forster and H. Breitfeld, Zeitschr.f. Mftallk., XXX (1938), 343. 

9. H. B. Huntington, Pkys. Rev., LXXII (1947), 321. 

10. W. P. Mason, Jour. A.S.A., XIX (1947), 464. 

11. R. Wegel and H. Walther, Physics, VI (1935), 141. 
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a number of investigations.'- A stronger coupling of the specimen 
with the driving and detecting system is obtained if one employs pole 
pieces glued onto the specimen rather than eddy currents. Such pole pieces 
have been employed by the author'^ in cases in which the internal friction 
was unusually high. .\n ingenious method has been adopted by Forster" 
for exciting and detecting oscillations. There, as in most methods, the speci¬ 
men is slung horizontally by two wires. In most methods the two ends of 
each wire are attached to rigid supports. In Forster’s method one end of 
each wire is attached to the movable membrane of an electromagnetic 
drive or pickup system. The vibrations are thereby excited and detected 
through the two supports. Still another method of exciting and detecting 
vibrations has been adopted by Quimby and his school.'^ In this method 

the specimen is cemented onto a quartz crystal, which is then excited pie- 
zoelectrically. 

The measured value of the internal friction is always raised above its 
true value by the dissipation of energy to the surroundings. One source of 
energy dissipation is acoustical radiation. Such dissipation can always be 
reduced to negligible values by placing the specimen in an evacuated 
chamber. The author has found such a precaution necessary when the 
internal friction is as low as 10"^ 

As mentioned above, in most electrical methods of excitation and de¬ 
tection the specimen is slung horizontally on two wire supports. If the 
specimen suffered no motion at the supports, the dissipation of energy 
through the supports could be reduced essentially to zero. However, owing 
to the Poisson contraction, the surface of the specimen undergoes some 
motion even at the nodes during either transverse or longitudinal vibra¬ 
tion. When proper precautions are taken, such dissipation of energy is 
however, very small, as is evidenced by the fact that measured values of 
tan 5 in transverse vibrations have been as low'* as 2 X 10“*. In order to 
avoid all loss through the supports, Frommer and Murray'^ have de¬ 
veloped special supports for use in torsional vibration. 


12. H. Randall, F. C. Rose, and C. 2^ner. Pkvs Rev I VT noWi n w j 

R. H. Randall. Trans. A.I.M.E., CXXXVII (194^, 41. ^ 

13. C. ^ner, D. Van Winkle, and H. Nielsen, Trans. CXLVII (1942) 98- and 

C. Zener, Trans. A.I.M.E., CLII (1943), 122. ^8, and 

14. F. Forster, Zeilschr.f. Metallk., XXIX (1937), 109. 

15. S.L. Quimby, /?«»., XXV (1925), 558, and XXXIX M9391 70 i. • 

Rrv XLIV (1933), 116; S. Siegel and S. L. Quimby, 

Cooke, Phys. Ra,.. L (1936), 1158; T. A. Read, Phys’. ' 

16. C. Zener, H. Clarke, and C. S. Smith, Trans. A.I.M.E., CXLVII (1942), 90. 

17. L. Frommer and A. Murray, Jour. Inst. Metals, LXX (1944), U. 
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In the dynamical method the internal friction of a single specimen may 
be measured over a range of frequencies by working at the fundamental 
frequency and at its first few overtones. The relation between these fre- 
fiuencies is trivial in the case of longitudinal and torsional vibration. The 
relation is given for the case of transverse vibrations in Table 8. 

TABLE 8 

Freque.vcies in Tr.\ns- 

VERSE ViBR.ATlONS* 

« ‘■it ‘■o 

0 1 

I 2.76 

2. 5 41 

8.94 

4. 13.37 

* A, Kalahn, Hanlbu'ft tier Phy'tik, 
iOl-2. 


In this dynamical method of free oscillation it is usually desirable to 
minimize support losses by suspending the specimen at its nodes of vibra¬ 
tion. The formula for the position of these nodes is trivial in the case of 
longitudinal and torsional vibration. The position of the nodes, x, in the 
case of transverse vibration is given by Table 9 in terms of the length of 
specimen, /. 

TABLE 9 

Position of Nodes in Tr.ansverse Vibration (j; /)• 


n 

1 

1 0 

1 1 

1 

' 2 

4 

1 

i i 

4 


0.224 

0 132 1 

0.094 

i 0.074 

i 0.060 


0.776 

.5 ! 

356 

1 .277 

.227 

X l 

1 . . . ) 

0.868 

.644 

. 5 

409 

1 \ 

\ 

0 906 

1 

' 753 

.591 




1 

0 927 

.774 

1 

! 

• 

• 


• 

0.940 


* A. Kaliihn. llaitdbucfi der Physik. VIII. 201-2 


C. VARLATION OF FREQUENCY 

In determining the relaxation spectrum, one must measure the internal 
friction over a frequency range sufficiently broad to include at least one 
relaxation band. As may be seen from Figure 21, the frequency range 
must cover at least two cycles of 10 in order that the internal friction in a 
band may be reduced on each side to 20 per cent of its maximum value. 
Such large frequency ranges have been employed.** Measurements over 
a large frequency range entail considerable experimental difficulty, and 

18. K. Bennewitz and H. Rotger, Pkys. Zeilschr., XXXVII (1936), 578; and Zeitschr.f. 
tech. Phys., XIX (1938). 521. 
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it is usually desirable to circumvent this difficulty through use of some 
indirect method. 

As discussed in chapter v, in certain circumstances the internal friction 
may be obtained as a function of frequency by a single experiment upon 
stress relaxation at constant strain. The relation between tan 6 and stress- 
relaxation measurements is given by equation (69). This method has been 
used with success by Ke'® and is applicable whenever the rela.vation band 
is much broader than for the case of a single relaxation time. 

The stress-relaxation method has the inherent defect that determina¬ 
tion of tan 6 entails differentiation of an experimental curve, namely, the 
stress-relaxation curve. Differentiation of an e.xperimentally determined 
curve necessarily introduces loss of accuracy as well as loss of detail. It 
is therefore usually desirable to employ an alternative indirect method. It 
is sometimes possible to predict, from theoretical considerations, that the 
internal friction will be a function of frequency v and of a time of relaxa¬ 
tion r in the form 

tan 5 = / (vt) , 

where / is some unknown function. This function may be determined as 
well through a variation in r as through a variation in v. Thus a plot of 
tan 6 versus log r is identical, aside from a horizontal shift of the abscissa, 
to a plot of tan 6 versus log y. 

The above indirect method has been used'^** to study the relaxation spec¬ 
trum associated with the stress-induced intercrystalline thermal currents 
in polycrystalline specimens. As discussed in Section A-4 of chapter vii, 
the time of relaxation r may in this case be set equal to d', D, where d is 
the mean grain diameter and D is the thermal diffusion coefficient. An 
effective frequency range covering five cycles of 10 was obtained by vary¬ 
ing d, i.e., by examining specimens of different mean grain diameter. 

This indirect method has also been applied’’ to several cases in which 
the time of relaxation r is suspected to vary with temperature according 
to an Arrhenius equation, namely, as 


One thereby obtains a plot of tan 6 versus in v by making measurements 
at one frequency over a range of temperatures and by plotting the results 
versus {U/RT). The heat of activation H is obtained by making observa- 


19. T. S. Ki. Phys. Rev., LXXl (1947), 533. 

20. R. H. Randall, F. C. Rose, and C. 2^ner, P/iys. Rev., LVI (1939), 343 . 

Rhysica, IX (1942), 862; C. Zener, Trans. \f F pi rr rioiii 

.Rct^LXXI (1947), 34;L. Dijlcstra, 11^19471 3S7 T S Kfe' 

Phys. Rev., LXXI (1947), 533, and LXXII (1947), 41. U> 4 /).s.->/, I. b. K 6 . 
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tions at two frequencies. Thus, suppose tan 5 versus l/T curves are ob¬ 
tained at the two frequencies n and V 2 . Then, if tan 6 is indeed a function 
of the product vt where r obeys an Arrhenius equation, these two curves 
may be superimposed by a horizontal shift. If A{X/T) is the magnitude of 
this shift, then the heat of activation H is determined from the relation 

Inih/fi) = (H/R) -Ad/r). 

As an example, an effective frequency range of three cycles of 10 will be 
traversed in raising the temperature from 0° to 60° C. when the heat of 
activation is 20,000 cal/gm/mole. 




VII 

PHYSICAL INTERPRETATION OF ANELASTICIT^' 


A. HOMOGENEOUS RELAX.A.TION 

IN NO real metal is the strain in the preplastic range a function of stress 
alone. It changes also with temperature, with composition, and in some 
metals with magnetic field or perhaps degree of order. In nonconducting 
solids the strain may also be changed by an electric field. The dependence 
of strain upon thermodynamically specifiable variables other than stress 
is a common source of anelasticity. 

As an introduction to a general discussion, the well-known example of 
thermoelastic coupling will be reviewed in detail. The results of this re¬ 
view will be applicable to thermodynamical variables other than tem¬ 
perature. 

For the present we shall consider that the strain depends only upon the 
two variables, stress and temperature. The strain, referred to zero stress 
and to a standard temperature as the reference configuration, will then 
be given by 

e XAT . ( 89 ) 

In this equation the relaxed modulus, Mr, is also known as the “isother¬ 
mal modulus,” the coefficient X is the linear thermal expansion coefficient, 
and AT is the deviation of the temperature from the standard tempera¬ 
ture. Obviously, a relation cannot be obtained between stress and strain 
unless the variation in temperature is specified by an independent equa¬ 
tion. Now two phenomena induce temperature changes. The first of these 
is diffusion, that is, the equalization, or relaxation, of temperature fluctua¬ 
tions. In many important cases, discussed in detail below, the change in 
temperature due to diffusion may be expressed to a very good approxima¬ 
tion by the following simple relaxation equation: 

= T 'AT. ( 90 ) 

The quantity r will be known as the “relaxation time.” The precise value 
of the relaxation time will depend upon the restrictions placed upon stress 
or strain during relaxation. If the stress is rriaintained constant and the 
strain relaxes with temperature, we shall have one relaxation time. On the 
other hand, if the strain is maintained constant and the stress relaxes with 

69 
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temperature, we shall have another relaxation time. These two values will 
be denoted by and t,, respectively. The second phenomenon that af¬ 
fects temperature is a change in strain which arises from sources other 
than the temperature change itself. Thus, if a rise in temperature induces 
an increase in length, conversely an adiabatic increase in length must be 
associated with a decrease in temperature. This second change in tempera¬ 
ture may therefore be written as 


where 



(91) 


Upon combining equations (90) and (91), we obtain 



(92) 


where, in accordance with our prior definition of the two types of relaxa¬ 
tion time, T, denotes the relaxation time for constant strain. 

An equation governing the relation between stress and strain is now de¬ 
rived by eliminating AT between equations (89) and (92). One thus ob¬ 
tains 

M ft e -\-MuTt e = <T -j- , (93) 

where the quantity A/tr is defined by 

Ma= (i+>^y)Mft (94) 

and, from the form of equation (93), denotes the unrelaxed modulus. One 
may therefore interpret r, as the relaxation time for stress relaxation, as 
well as for temperature relaxation, at constant strain. It therefore seems 
that r, may be the relaxation time for strain relaxation, as well as for tem¬ 
perature relaxation at constant stress. If this is true, then equation (93) 
may be written as 

jV/r ( e r^e) = 0-+ , (95) 

and the relation between the two times of relaxation is 

T<r _ Mu 

The identity of the r, defined in this manner with the original definition of 
relaxation time for temperature at constant stress may be verified by the 
elimination of e and e between equations (89) and (95), leading to 

^ AT = — Tw ^AT — yM u ■ 
at 


(96) 
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We conclude that, whenever the temperature fluctuation is governed 
by equation (92), the thermoelastic coupling changes the solid from a 
Hooke solid into a standard linear solid, whose anelastic properties were 
discussed in detail in chapter v, section A. Thus the dependence of internal 
friction and elastic modulus upon frequency is given by Figure 21, while 

the relaxation of stress at constant strain, or of strain at constant stress, 
is given by Figures 18-20. 

Considerable physical insight into the reason why the internal friction 
has a maximum in the vicinity of the angular frequency w 1, r and ap¬ 
proaches zero asymptotically on either side of this maximum may be 
gained from an analysis of the precise mechanism whereby energy is dis¬ 
sipated as heat. If a quantity of heat, 5Q, is added to a unit volume at 

temperature T, the entropy, 5, of the unit volume increases by the 
amount 



The total increase in entropy during a cycle of vibration is therefore 

where the summation extends over all increments during one cycle of vi¬ 
bration. Once a steady state is reached, the specimen suffers no net in¬ 
crease in entropy, the net flow of heat away from the specimen exactly 
balancing the dissipation of energy within the specimen. If is the mean 
temperature of the specimen, and 6T is its deviation from the mean tem¬ 
perature at any one instant, then the equilibrium condition, 

A5 = 0 , 

may be re-written, with the help of equation (97), as 


AQ= ST 

where the net outflow of heat per cycle is 


(98) 


AQ 2 50 . 


I'rom equation (98) we see that the energy dissipated per cycle contains 
a summation of products of the type 50 • 5r. In the very high frequency 
range, where the vibration is essentially adiabatic, the energy dissipated 
is very low because of the smallness of the first factor, 50. In the very 
low frequency range, where the vibration is essentially isothermal, the 
energy dissipated is very low because of the smallness of the second factor, 
ST. Only in the intermediate frequency range, where the vibration is 
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neither essentially adiabatic nor isothermal and where both 8Q and dT 
are sensibly different from zero, is the dissipation of energy appreciable. 

From equation (94) we obtain the following formula for the relaxation 
strength, as defined in equation (71): 

A = \y. (99) 

An alternative expression will be derived in terms of quantities which are 
more directly measurable than y. If we regard <r and T as the independent 
variables and e and S as the dependent variables, then 

de=M'^^da-\-\dT , (100) 

dS = \dtx-\-iidT, (101) 

where 

..(gV (102) 


The symmetry^ of the coefficients follows from the fact that the increment 
of the free energ>' density is a perfect differential. Thus from 


it follows that 


(f (£ -* <7 e — TS) = — ed <T — SdT , 



(103) 

(104) 


Upon setting dS equal to zero in equation (101), 


we obtain for M's 


and therefore 



(105) 


Still another expression for the relaxation strength may be derived by 
observing from equations (100) and (101) that 

\dSA M* 


Substitution of this relation in equation (105) leads to 


I. THERMODYNAMICAL VARIABLES 

As previously intimated, temperature is only one of several thermo¬ 
dynamical variables in which fluctuations are introduced by strains and 


Pm^SlCAL INTERPRETATION OF ANELASTICITY 73 

for which the relaxation of such fluctuations gives rise to anelasticity. In 
the present section the various types of thermodynamical variables are 
discussed. In the following sections the various mechanisms of relaxation 
are analyzed. 

As a guide in the formulation of a general theory' of relaxation, we shall 
apply thermodynamical considerations. If a system of unit volume is 
completely specified by stress and temperature, the increase in energy as¬ 
sociated with a change in state is given by 

dE== ade-\-TdS , (106) 

where S is the entropy of the system, i.e., the entropy of an amount of ma¬ 
terial which has unit volume at zero stress and at a standard temperature. 
Now T may be regarded as a thermodynamical potential, since a system 
is in equilibrium only when everywhere T has the same value. On the 
other hand, the entropy per unit volume may be regarded as a thermo¬ 
dynamical density, since the entropy of a system is equal to the entropy 
of its component parts. Now when more variables are necessary complete¬ 
ly to specify the system, it is possible to generalize equation (106) to the 
form 

dE = <Tde-\-'ZjAjdai , (107) 

where all the A’s may be regarded as thermodynamical potentials and the 

a*s as thermodynamical densities. Equation (89) may then be generalized 
to 


e — M R <T -f- 2yXjA A /. 


(108) 


In this equation each AA,- is the deviation of the corresponding potential 
from a standard value, e is the strain referred to the standard state in 
which o- = 0 and all the potentials have their standard values, and the co¬ 
efficient X/ is defined as 



(109) 


the partial derivative denoting that the stress and all the potentials other 
than Aj are to be maintained constant during the differentiation. 

In many cases thermodynamical potentials other than temperature 
obey a relaxation equation such as equation (92), i.e., 

dLAj -1 

-ry AAy-yye. (HO) 


The relaxation of such thermodynamical potentials contributes to the 
anelasticity of the specimen in precisely the same manner as does the re- 



74 


ELASTICITY AND ANELASTICITY OF MET.ALS 


laxation of temperature fluctuations. In particular, the relative difference 
between the associated relaxed and unrelaxed moduli is given by 

‘ =X>T>, (111) 

or, in analog)’ to equation (105), by 


where 



(lllfl) 

(lllA) 


In deriving this equation it has been implicitly assumed that the same 
thermodynamical potentials or densities, other than Aj or a>, are held 
constant in the measurement of as in AIa - 

In a specimen which contains a mixture of elements all in solid solution, 
strain induces changes in the partial thermodynamical potentials of each 
constituent. These partial thermodynamical potentials, sometimes called 
“chemical potentials,” will be referred to unit volume rather than to a 
gram-mole, as is more usual. They are defined as follows: Let G be the 
Gibbs free energy of an amount of material which contains unit volume 
under standard conditions, i.e., at zero stress and at standard composition. 
Thus 

G=E- ea-TS . (112) 


Then the partial thermodynamical potential, 6’/, will be defined as 



( 113 ) 


where «/ is the number of atoms of type j and where, in the differentia¬ 
tion, <r, T, and the number of all other types of constituent atoms are to 
be maintained constant. Upon solving equation (112) for E and differen¬ 
tiating, we find, therefore, 

dE = ade + TdS + 'LGjdm . (114) 


The G;’s and the deviation of the «>’s from the standard concentrations 
may therefore be regarded as thermodynamical potentials and densities, 
respectively. The Gj at each part of a specimen will be changed by 
strain. Since equilibrium is attained only when the Gy’s are constant 
throughout the specimen, a fluctuation of the Gy’s will gradually be equal¬ 
ized by diffusion. Thus relaxation of the fluctuations in the Gy’s will be 
accompanied by the segregation into the regions under extension of those 
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atoms which increase the lattice constant and, conversely, the segregation 
into those regions under compression of those atoms which decrease the 
lattice constant. 

In ferromagnetic materials the coupling between strain and external 
magnetic field is well known. In many cases of importance the external 
magnetic field is essentially linear, and the axis of the specimen is parallel 
to this field. If the transverse dimensions of the specimen arc very small 
compared with the wave length of the strain within the specimen, then 
the macroscopic magnetic field within the specimen may be regarded 
everywhere as parallel to the axis of the specimen. Application of the Max¬ 
well equation. 


cur! 11 = 0.4W+ C-* 

at 


(115) 


to this case shows that under equilibrium conditions, when i and dD dl 
are both zero, the macroscopic magnetic field strength. H, is constant 
across the specimen and is equal to the external magnetic field strength. 
The macroscopic magnetic field strength, H, may therefore, under the 
above conditions, be regarded as a thermodynamical potential. In order to 
find the associated thermodynamical density, we note that changes in the 
magnetic energy density are given by 

^dE) iiinKIivlic ~ "i IIdB , 

47r 


where B is the magnetic flux density. Thus, if H is regarded as a thermo¬ 
dynamical potential, B Air must be regarded as the associated density. 
The relaxation of the magnetic strength fluctuations induced by strains 
is one important cause of anelasticity. 

One cannot, strictly, speak of the “temperature” of a specimen unless 
all parts of the specimen are in thermal equilibrium, i.e., unless the tem¬ 
perature is uniform throughout. Likewise, one cannot even speak of the 
temperature of an element of volume unless all degrees of freedom within 
the element are in thermal equilibrium. In the previous discussion of 
temperature relaxation and of the anelasticity associated therewith, it was 
implicitly assumed that all degrees of freedom within each element of 
volume were in thermal equilibrium and hence that a temperature could 
be associated with each element. In many important cases such an assump¬ 
tion is not valid. Rather it is necessary to divide the normal coordinates 
into two groups, the normal coordinates within the same group being in 
thermal equilibrium. Anelasticity then arises through the relaxation of the 
temperature difference between the two groups of coordinates. One group 
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contains the coordinates associated with thermal vibration of the atoms. 
The other coordinates may determine the distribution of atoms on the 
lattice positions in case more than one type of atom is present. Such 
coordinates are known as “ordering” coordinates. In the case of molec¬ 
ular crystals the second group of coordinates may refer to the coordinates 
of the molecules with respect to their center of gravity. With each group 
of coordinates there is associated a corresponding entropy density. 

II. ORTHOGONAL THERMODYNAADC POTE>mALS 

In the above discussion of the anelasticity associated with relaxation, 
it was explicitly assumed that the thermodynamic potentials satished a 
simple relaxation equation of the type of equation (HO). Such an assump¬ 
tion is probably valid for the relaxation associated with the establishment 
of thermal equilibrium between the various coordinates of an element of 
volume or for the relaxation of a temperature fluctuation between a speci¬ 
men and the surrounding medium when the specimen may be regarded as 
at a uniform temperature throughout. In those cases, however, in which 
the rela.xation relates to fluctuations in thermodynamic potentials between 
different regions of the specimen, the simple relaxation equation cannot be 
applied without further modification. 

One possible modification is to replace the thermodynamic potentials, 
which are functions of position, with an infinite series of normalized orthog¬ 
onal functions, the latter satisfying the appropriate diffusion equation 
and boundary' conditions. This procedure, originally developed by the au¬ 
thor,* is described below. An alternative procedure has been described by 
Pasler." The coefficient of each of the above normal orthogonal functions 
is a function of time and satisfies a simple rela.xation equation of the 
type of equation (110). These coefficients will have the dimensions of 
thermodynamic potentials and will therefore be called “orthogonal ther¬ 
modynamic potentials.” Each type of relaxation, be it thermal, concen¬ 
tration, or magnetic, will then involve the simultaneous relaxation of all 
associated orthogonal thermodynamic potentials. When the total relative 
relaxation is small compared with unity, i.e., when 

A«1 , 

the internal friction may be expressed as the sum of the contributions from 
the individual potentials. As an introduction to the general analysis of 
this method, an example of thermal relaxation will be presented. 

1. C. Zener, “General Theorj- of Thermoelastic Internal Friction,” Pkys. Rev., LIII (1938), 

90. 

2. M. Pasler, “Zur Theorie der thermische Dampfiing in fester Korpcr,” Zeilschr.f. Phys., 
CXXII (1944), 357. 
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For this example we shall consider the internal friction arising from re¬ 
laxation of temperature across a specimen vibrating transversely. The 
two equations which must be solved simultaneously are 

c = + XAT (116) 

and 

DV'-AT-ye. (117) 


The second equation replaces the simple relaxation relation of equation 
(92). In order to derive equations in which the spatial coordinates are 
eliminated, we expand e, a, and AT in a complete set of normal orthogonal 
functions. Thus, 

C =2jey(/) Uj , (118) 

=2ja>(/) Uj, (119) 

A7' = 2yr>(0 Uj. (120) 

These functions will be chosen as the characteristic solutions of the differ¬ 
ential equation 

( 121 ) 

and of an appropriate boundary condition. Such solutions exist only when 
T has discrete values, known as the “characteristic relaxation times.” 

In the particular case under consideration, strain may be regarded, to a 
close approximation, as a function only of the transverse dimension of the 
specimen. Equation (121) may therefore be replaced by 

(■^+’-‘)f/ = 0. (122) 


the coordinate x being taken as this transverse dimension. The boundary 
condition will be taken to correspond to no heat flow across the boundaries. 
If the transverse width of the specimen is a, the boundary condition will 
then be 


i£. 

dx 


0atX = ±|. 


The characteristic normalized solutions and relaxation times for this prob¬ 
lem are 

C/t = 'v/?sin I (2k + \)^\ 

a { a ) 

Tk fl* ’ 


and 
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respectively, where k is any one of the sequence of numbers 0 , 1 , 2 , 

3 f « • • a 

If we now multiply equations (116) and (117) by U^, integrate over the 
volume of the specimen, and make use of equation ( 121 ), we obtain the 
following two equations: 

Ck = M R ctc + Wk (123) 

and 

Tk = - rVlk-yCk . (124) 

Each orthogonal thermodynamic potential, Tk therefore satisfies a simple 
relaxation equation. 

Now the rate of dissipation of energ>' by the specimen is given by 

E = ^aedv , (125) 

the integral extending throughout the volume of the specimen.. Upon mak¬ 
ing use of equations (118) and (119), we arrive at 

E = Z,a,es. (126) 

If we now suppose that c, a, and T vary periodically with time and if wc 
take the time average, we obtain, uj)on using equation (123), 





(127) 


Taking a» as the 
tion (124) 


angular frequency of vibration, wc derive from equa- 




(128) 


In the case in which Mr/M u is only slightly less than unity, the internal 
friction may be written as in equation (80), where A£ is the energy loss per 
cycle and E is the average energy of vibration. Now 

SE = — E. (129) 

O’ 


and, when Mr/M u is only slightly less than unity. 


E = MnSe^dv=MRZie]. 

Upon combining equations (127)—(130), we obtain 




' i l-fa^r? 
27 ? 

f 


(130) 


tan 5 = X 7 


(131) 
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If we now make use of equation (94) and define the quantity 


e- 

t- 


(132) 


we transform equation (131) into 


tan 6 = AsS,7.- 


a>r 


1+0)272 • 


(133) 


By the manner in which they are defined, the quantities /; automatically 
satisfy the following condition: 

2/>=l. (134) 

The evaluation of internal friction due to relaxation may thus be re¬ 
duced to the problem of finding the characteristic relaxation times, 
and the coefficients, />. The characteristic relaxation times have already 
been found for the special case under consideration. The coefficient/> is 
given by 




/.= 


/ 2 \l /2 

( - ) / .V sin (2 7 + 1) (ir.v/a) dx 

\a/ J-a/i 


(135) 


/ X'dx 
= 967r-M2i+l)-'. 

The first few values of the coefficients are given by 

/o = 0.986, /, = 0.012, /2 = 0.0016 . 

In this particular case, only a very slight error would be made by setting 
/o equal to unity and all the following/’s equal to zero. The thermal relaxa¬ 
tion across a specimen vibrating transversely may therefore be regarded 
as essentially a simple relaxation with 

Relaxation lime = a^lir'^D\ a = Width. 


III. GENERAL THEORY OF RELAXATION BY DIFFUSION 

In the preceding section the concept of orthogonal thermodynamic po¬ 
tentials was introduced, and it was shown how the internal friction arising 
from diffusion could be computed by use of such potentials. In the present 
section the results of such computations in particular cases are sum¬ 
marized. 


It has been shown that each orthogonal thermodynamic potential re¬ 
laxes in a simple manner and therefore contributes an internal friction 
which varies with frequency in the manner shown in Figure 21. The con¬ 
tribution from all the orthogonal potentials may therefore be expressed as 
in equation (133), namely, 


tan 6 = A/j 



OinTn 


(136) 
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Each coefficient/„ may be regarded as the relative strength of the asso¬ 
ciated orthogonal potential. The sum of all the relative strengths is unity, 
as indicated by equation (134), The relative strengths are computed as 
follows. Let eo (x, y, s) be the strain when the amplitude of the vibration 

is at a maximum. Further, let U„ (x, y, z) be the wth normalized orthog¬ 
onal potential. Then 


_ (fl/neodv)^ 
feldv 


(137) 


The orthogonal potentials are the characteristic solutions of the equation 

(VV^+T-^) (138) 


and of the associated boundary conditions. The characteristic relaxation 
times, r„, are those values of t for which such solutions exist. 

In one important class of diffusion problems all terms in the summation 
in equation (136) are negligible compared with the first. Such was the case 
in the example cited in the previous section. In this class of problems the 
diffusion is essentially a volume diffusion, there being no diffusion across 
surfaces at which the properties of the medium suffer a discontinuous 
change. 

In a second important class of problems, successive terms in the sum¬ 
mation in equation (136) converge very slowly, particularly at high fre¬ 
quencies. In such problems the slow convergence is due to the presence of 
surfaces across which diffusion takes place and across which the physical 
properties of the medium change discontinuously. In the absence of all 
diffusion the thermodynamic potentials would suffer a discontinuous 
change across these surfaces during vibration. Therefore, at increasingly 
high frequencies, where diffusion per half-cycle becomes less and less, the 
gradient of the potentials across these surfaces becomes steeper and steep¬ 
er. At increasingly high frequencies the diffusion therefore becomes lim¬ 
ited essentially to the immediate vicinity of these surfaces and is normal 
thereto. Two alternative methods may be used for computing the internal 
friction in those cases of high frequency when the convergence of the sum¬ 
mation in equation (136) is poor. One may convert the summation into 
an integral, thus 

tan d = A^//„ , . dn . (139) 

On the other hand, one may take cognizance of the fact that under these 
conditions of high frequency the diffusion is essentially confined to the 
discontinuous surfaces and must therefore be proportional to the area of 
such surfaces per unit volume. Thus the only computation required is the 
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evaluation of the dissipation of energy per unit area of such surfaces. This 
computation is given in the appendix to this section. The result is 

tan 3 = (140) 

The precise value of the numerical coefficient /S depends upon the type of 
system. If the only surface is the boundary of the specimen, then 

^ _ ( Cq) surface average 

(ej) volume average' (141) 

If, on the other hand, internal boundaries are present across which prop¬ 
erties change discontinuously, such as grain boundaries, the evaluation 
of ^ requires special consideration. A particular case is examined below. 

a) DIFFUSION BETU'EEN SPECIMEN AND SURROUNDING MEDIUM 

i. Longitudinal vibration of specimen of circular cross-section} —The rel¬ 
ative strength,/„, is given by 



where q„ is the ;ith root of 

J<i{q) - 0 . 

The characteristic relaxation time, r„, is given by 

_ 

where a is the radius of the specimen. From the formula^ 

32’ 

one obtains the asymptotic relation, 

D 

tan 5 = — 

while from equations (140) and (141) one obtains the asymptotic relation, 

(143) 

These asymptotic relations were originally derived by Kersten."^ 
ii. Longitudinal vibration of reeds. —The relative strength is given by 

/ - -- ^ 

{2n-\-\)V 

1010 ^ “General Theor>' of Macroscopic Eddy Currents,” Phys. Rev., LIII (1938), 

4. G. N. Watson, Theory of Bessel Functions (Cambridge, 1922), p. 502. 

5. M. Kersten, Zeitschr.f. tech. Phys., XV (1934), 463. 
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The characteristic rela.\ation time. r„. is given by 

4a2 


( 2 « + 1 ) 2 ^ 2 /)’ 

where a is the transverse width of the reed. 

From the formula*^ 

one obtains the asymptotic relation, 

,an 5=A..3^. a,«^, 

while from equations (140) and (141) one obtains 

tan 

\la-w/ or 


(144) 


(145) 


iii. Transverse vibration of specvnen viith circular cross-section- 
relative strength./„, is given by 

f 

where is the «th root of 

JAq) = 0 . 

The characteristic rcla.xation time, r„, is given by 


-The 


a- 


Tn - 


From the formula 






1 

192 ’ 


one obtains the asymptotic relation, 


tan 5 = A 


a-o> 


B 


1\D' 


while from equations (140) and (141) one finds 


tan h 






D 


a- 


(146) 


(147) 


iv. Transverse vibration of reeds .—The relative strength,/„, is given by 


6. E. P. Adams, Sniitksonuin Xfathematical Formulae (Washington. D.C.: Smithsonian 
Institution, 1922), p. 140. 

7. Watson, op. cit., p. 502. 
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while the characteristic relaxation time, t„, is given by 

a~ 

where 2a is the transverse width. 

From the formula’^ 



one obtains the asymptotic relation, 

a-u D 

tan 5 = w«—. (US) 

i ^ u a- 


From equations (140) and (1411 one finds the asymptotic relation 

/ D D 

tan 6 = Aa-* 3 (—-) , w»—. ( 149 ) 

\Za~o}/ U‘ 


b) M XCROSCOPIC DIKFl'SION WITHIN SHhXIilES’ 


When diffusion is confined to the interior of the specimen and when 
the diffusion is essentially macroscopic in nature, the following summaiy- 
shows that no appreciable error is made by retaining only the first term 
in the infinite summation in equation (136). The thermodynamic poten¬ 
tial giving rise to diffusion may thus be regarded as obeying a simple re¬ 
laxation equation. As we have previously seen, the analysis for the cou¬ 
pling between strain and stress, on the one hand, and the thermodynami¬ 
cal potential, on the other hand, may in this case be carried out exactly, 
irrespective of the magnitude of Mv/Mr. According to equation (56), the 
internal friction may therefore be written as 


where 


tan 


M v — M H iOT 


T = (r,T,) . 


On the other hand, the elastic aftereffects may be obtained precisely 
through the solution of equation (93) or equation (95). Therefore the re¬ 
laxation of stress following sudden application of a constant strain is given 
by equation (46), the relaxation of strain following the sudden application 
of a constant stress is given by equation (48). 

i. Longitudinal vibralion.—fo — 1, precisely. 

_ D 

To ^ 2 , 

where C is the velocity of propagation. 

8. .\dams. ol>. cif. 
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ii. Transverse vibration of reeds.^ —The relative strengths are given by 

/„= 967r-‘‘(2« + l)-S 

the first few values of which are 0.986, 0.012, 0.0016. The longest char¬ 
acteristic relaxation time is given by 

To = ^a^Tr~^D~^ ; 

successive values are given by 

7 - - 

” (2tt + l)2- 

iii. Transverse vibration of specimens with circular cross-section}^ —The 
relative strengths are given by 

where is the wth root of 

J\{q) =0. 

The first few values are 0.988, 0.010, and 0.0015. The characteristic re¬ 
laxation times are given by 

_ 0 ^ 

the first few values of which are 1.85 aW, 0.221 aW, and 0.117 a}D. 


C) INTERGR/VIN DIFFUSION 

In polycrystalline specimens of noncubic metals the elastic anisotropy 
and at least partial random orientation result in a fluctuation of stress and 
strain from grain to grain. Application of a stress will therefore result in 
a fluctuation of the thermodynamic potentials from grain to grain, a fluc¬ 
tuation which is relaxed by intergrain diffusion. Such intergrain diffusion 
will give rise to characteristic anelastic effects. 

No exact computation can be made for the difference between the re¬ 
laxed and the unrelaxed moduli associated with intergrain diffusion. At 
best, an estimate may be obtained. Suppose, for this purpose, that the 
stress in each grain is uniform and is equal to the stress in every other 
grain and is a uniaxial tensile stress. The rise in thermodynamic potential 
due to the adiabatic application of this stress is, in every grain, equal to 



9. C. Zener, “Theory of Internal Friction in Reeds,” Pkys. Rev., LII (1937), 230. 

10. C. Zener, “General Theory of Thennoelastic Internal Friction,” Phys. Rev., LIII 
(1938), 90. 
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From the equations 

d € — -f- Xrf.-l , 

da = \d<j nd A , 


which are analogous to equations (100) and (101), we obtain 



and therefore 





The coefficient X will be dependent upon crystallographic orientation in 
noncubic metals. Therefore w'ill vary from grain to grain. The devia¬ 
tion of from the mean is, according to the above assumption of uni¬ 
form stress, given by 

= “■ (X — X) . 

We now compute the strain in each grain which will result from the equal¬ 
ization of potential fluctuations. The strain in each grain will be given by 

Ae= -(|j) = X(X-X)m-V . 


The over-all strain resulting from the equalization of potential differences 
between adjacent grains will be the average of this strain, namely, 


Ae= (X2-X^)m-V . 

The change in the modulus due to local rela.\ation betw'een grains may 
therefore be written as 


where 


Mr' -Mv ^ 



2 


X2-X 



(150) 

(151) 


The first factor in equation (150) is the change in AT”* which would result 
from complete relaxation of all changes in the thermodynamic potential. 
The second factor, R, is that fraction of this change which takes place by 
intergrain diffusion, that is, by relaxation of the fluctuations in the po¬ 
tential which occur from grain to grain. 

In metals of trigonal and hexagonal symmetry X varies with orienta¬ 
tion as 


X = Xii cos® 0 + Xi sin® 6 . 


(152) 
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Substitution into equation (151) leads to 


^ —^l) ■ (5 Xii + 4Xn Xx + 8X|) ’ . 


A similar analysis for polycrystalline specimens of cubic metals gives a 
zero fluctuation of thermodynamic potentials between adjacent grains, 
since in cubic metals the coefficient X is independent of crystallographic 
orientation. A more detailed analysis, in which we no longer assume a uni¬ 
form stress, shows that even in cubic metals a sudden application of stress 
will result in fluctuations of thermodynamic potentials between adjacent 
grains, the relaxation of which will give rise to anelastic effects. Thus the 
constraints imposed by neighboring grains prevent uniaxiality of stress 
within the grains. As an example, consider a grain boundary normal to the 
tensile a.xis. The Poisson contraction normal to this axis will be different 
for the two grains on either side of the boundary. At the boundary, how¬ 
ever, they must contract equally. This boundary constraint is accom¬ 
panied by transverse stresses at and in the vicinity of the boundaries. 

I he inequality of these transverse stresses results in an inequality in the 
thermodynamic potentials and hence in diffusion accompanied by anelas¬ 
tic effects. Any estimate of the magnitude of the difference between the 
relaxed and the unrelaxed moduli will lead to an expression of the type in 
equation (150), where now the factor R is the relative mean square devia¬ 
tion of some function of the elastic constants. For ease of computation 
we shall take the following: 


R^ 

Upon using the relation” 


{(£-1 _(£-!) 




av 


(153) 



where yi, y^, and 73 refer to the three direction cosines of the axis of the 
specimen with respect to the crystallographic axes, one obtains 


where the quantity 

-(jii— iis)— J44 /1CC1 

q = 2 - 

is a measure of the elastic anisotropy of the grains. The quantities g and R 
for a number of cubic metals are tabulated in Table 10. 

Irrespective of whether the polycrystalline specimen is of a noncubic 

11 . £. Schmidt and W. Boas, Kristdlplastizitat (Berlin: Springer, 1935), p. 23. 
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or of a cubic metal, at sufliciently high frequencies dilTusion will be con¬ 
fined essentially to the immediate vicinity of the interfaces between 
grains, across which the elastic constants change discontinuously. In this 
case an approximation to the internal friction arising from this diffusion 
may be obtained by use of equation (140). The grain-boundary area per 
unit volume is inversely proportional to the mean linear dimensions of the 
grains. We thus obtain 

(D 


1/2 


tan 


G.S. 


(156) 


where G.S. denotes any reasonable measure of grain size, e.g., the grain 
size determined by conventional methods. 


TABLE 10 

InTERGR.MN REL.WATION CONST.AN'TS 

FOR Cubic Met.als 


MetttI 1 

1 

1 

1 

g 

R 

1 

^ brass . 

1 2.0 

0.18 

Na . 

! 2.5 

.10 

K.. ... 1 

2 44 

. 15 

Pb. 1 

2.2 

10 

a-brass (r0-.10) | 

2.1 

i 091 

Cu 

1.95 

071 

.\u 

1.8.? 

.067 

Ag 

1.80 

.00,? 

Fe.. 

1.00 

040 

A1. . 

0.51 

.0024 

W. 

0.00 1 

1 

1 • 

1 0.0000 

1 

1 


APPENDIX TO SECTION A, SUBSECTION III 


In this appendix the computation is given for the dissipation of energy attending 
the diffusion across a boundary, the diffusion distance being so small that the diffusion 
currents may be regarded as everywhere normal to the boundary. 

The medium will be regarded as bounded on one side by the plane = 0 and ex¬ 
tending indefinitely along the positive axis. The energy, AE, dissipated per cycle per 
unit surface is then given by 

AE = — I ae dx . (i) 

The stress and strain are related by the equation 


de=M-^d<T + \dA , 


where the potential A is given by the diffusion equatioi 



dA 

dt 



dM 

dA-2 



and by the boundary condition 
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Corresponding to the condition of periodic vibration, we shall set 


e = eo sin w/ . 

Then, by combining equations (i) and (ii), we obtain 

2 


(v) 


A£= - f .1 e d. 

w J<s 


(vi) 


In order to evaluate this integral, we must now solve equation (iii). Toward this 
end we introduce the auxiliary function. V{x, /), which satisfies the equation 




dt 


a.v2 


and the boundarv conditions 


r(o,/) = 0, 

V (.V, 0) = 1. .V > 0 . 
The solution of equation (iii) may then be written as 


.1 = U{x,t-n\-ye{n\dt' 

_ CO 


Upon transforming 


T = l-r , 


this equation becomes 


.1 


= —yf i’ {x. t) e {I — t) dT . 
•'0 


(vii) 


When wc now substitute equation (v) in equation (vii), multiply by e, and take the 
time average, we obtain 


Ae= - 


yoj-e' '** 


• ^00 

- / 6’ (.V, T) cos u)T dr . 

*'0 


Upon substituting this equation into equation (vi), we obtain 


A£= ~ irywell^IJ J T'’(.v, t)cos<j)T dr dx 


(viii) 


It is now necessary to substitute in this equation an explicit expression for V{x, r). 
This function is given in all standard treatises on heat conduction.** It is 


t/2 y Dt 


e~^d0. 


(ix) 


We first integrate with respect to t. After integrating by parts, we obtain 

f U {x, t) cos u}T dr = f sin ut dr . 

-'o 2 V tD -'0 

12. W. E. Byerly, An Elementary Treatise on Fourier Series (New York: Ginn & Co., 1893), 
Ex. 5, p. 83. 
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We next integrate with respect to .v, obtaining 


Jo Jo J J T-^^-s\n i>}T d 


© 1/2 


- 3/2 . 


and therefore 


//)\l/2 

lE=TT\y(^j-J Mh^'o. 


If we now denote by fi the total surface area and by V the total volume, the inter 
nal friction arising from this diffusion is 

tan (iA/fffS-K) 


TT 


or 




(x) 


IV. RELAXATION BY THERMAL DIFFUSION 

Thermal conduction provides the best-known mechanism of relaxation 
in metals. A rise in temperature at constant pressure is always accompa¬ 
nied by an increase in volume. Conversely, the adiabatic application of a 
tensile stress gives rise to a reduction in temperature and therefore tends 
to cause heat to flow into the specimen from its surroundings. As the tem¬ 
perature change is gradually relaxed, the specimen suffers a further slight 
increase in length. 

The thermal relaxation strength (Es-Er)/Er is given, according to 
equation (105), by 



where a is the linear thermal e.vpansion coefficient and is the specific 
heat per unit volume. The numerical value of Aj? for the common metals is 
given in Table 11. The time of relaxation for the establishment of tempera¬ 
ture equilibrium is given by 



(157a) 


where d is comparable to the distance that heat must flow for the estab¬ 
lishment of temperature equilibrium and D is the thermal diffusion co¬ 
efficient given by 
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where a is the thermal conductivity, C,., as above, is the specific heat per 

unit volume. The room-temperature values of D for the common metals 
are listed in Table 11. 

1 he importance of thermal conduction in damping sound waves was 
recognized by Kirchhoff’^ as early as 1868. In all types of longitudinal vi¬ 
brations, including sound, the diffusion distance is comparable with the 
wave length, X. Upon utilizing the relation between X. the frequency v, 
and the velocity of sound C 

C = v\ . 

T.VBLIC tl 

Thermal Rklax.-\tio\ Coxstaxt.s 


Metal 

D 

( S <|. Cm / Scc > 

\ ^ 

1 

Cd. 

0.46 

0.010 

Zn. 

0.41 

.0088 

Mg 

0.60 

.0050 

Al. 

0.88 

.0046 

Be . 

0 .. S 3 

.0046 

Sn. 

0.40 

.0040 

Brass ( 70 - 30 ). . . 

0.38 

.0036 

Ag. 

1 74 

.0034 

Cu. 

1.2 

.0030 

.\i. . 

0 l.> 

()029 

Vb 

0.24 

.0025 

I-e. . 

0.20 1 

.0024 

IM 

0.26 

.0020 

Sb. 

0.12 

.0018 


1.1 

.0017 

PI. . . 

0 26 

.0015 

Bi , . 

0.065 

.0014 

W 

0.61 

.00078 

Rh 

0.29 

.00069 

Ta . 

0.22 1 

1 

0.00030 


we find from equation {\57a) that the ma.ximum damping, which occurs 
when the period of vibration is comparable with the time of relaxation, 
appears for the wave lengths 


^ C' 


A precise computation gives 


X = 27r 


D 


for the wave length at maximum damping. Upon taking for air the value 
0.18 cmVsec for D and 33,000 cm/sec for C, we obtain for the optimum 
wave length 3.5 X 10“® cm. A similar computation for the common metals 


13. G. Kirchhof!, “Ober den Einfluss der Warmleitung in einen Case auf die Schallbe- 
wegung,” Poggetuiorf's Ann., CXXXIV (1868), 177. 
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leads to even smaller values for the optimum wave length. Artificially pro¬ 
duced longitudinal waves have, therefore, veiy little damping arising from 
the flow of heat along the temperature gradients produced by the waves. 
The waves are essentially adiabatic. 

Many cases arise, however, in which temperature fluctuations are jires- 
ent during vibration across distances which are ver>^ short compared with 
the wave length of a longitudinal wave of the same frequency. The first 
example studied arose in transverse vibrations. Here temperature relaxa¬ 
tion occurs through the flow of heat across the specimen. The appropriate 
formulae for this example are given on page 84. Here the relaxation 
has essentially a single time of rela.\-ation, so a plot of internal friction 
versus log frequency should lead to a symmetrical bell-shaped curve, as 
in Figure 21. The first recorded observations on the internal friction dur¬ 
ing transverse vibration, by Bennewitz and Rdtger,'' are reproduced in 
Figure 25. The vertical lines in this figure correspond to the positions of 
the maxima computed according to the theoretical formula, 


(diameter)^ ’ 


The agreement of the computed with the observed maxima is surprisingly 
good, particularly in view of the fact that Bennewitz and Rotger were un¬ 
aware of this formula. After learning of the theoretical interpretation’*’ of 
their original experiments, they performed very careful experiments’’’ upon 
one metal, German silver, in order to check the theory in all details, includ¬ 
ing the height of the maximum and the shape of the internal friction 
versus log frequency curve. Their results arc presented in Figure 26. Near¬ 
ly perfect agreement between theory and experiment was obtained. 

The success of the theory of damping of transverse vibrations by trans¬ 
verse thermal currents led to the development of the theoiy’' of damping 
by intercrystalline thermal currents and then to its experimental verifica¬ 
tion.The theory has been reviewed on pages 84-87. The transition from 
the case of isothermal to the case of adiabatic vibration, with respect to 


14. K. Bennewitz and H. Rotger, “fiber die innere Reibung fester Kbrper; .^bsorptions- 
frequenzen von Metallen in akustischen Gebiet,” Fhys. Zeitschr., XXXVTI (1936), .S78 

LlS aW7): Demonstration of Thermoelastic 

«-8-'''n.ingungon,-- 

(1938) Theory of Thermoelastic Internal Friction,” Phys. Rev., LIU 

18. R. H. Randall, F. C. Rose, and C. Zener. “Intercrvstallinc Thermal Currents as a 
Source of Internal Friction,” Phys. Rev., L\'I (1939). 343. ' currents as a 
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adjacent grains, was achieved by varying the grain size as well as the fre¬ 
quency of vibration. Thus the degree of adiabaticity can depend only upon 
the following quantities; frequency (p), grain size (G.S.), and thermal dif¬ 
fusion coefficient (D); and the only dimensionless combination of these 
variables is »'(G.S.)^ D. As shown in Figure 27, when the observations at 



V (c.ftec.) 


Fig. 26. —Experimental check of theorj’ of damping of transverse vibrations by thermal 
currents. Example of German silver, after Bennewitz and Rotger. Curve is given by theory 
with no arbitrar)- parameters. 


all frequencies and for all grain sizes are plotted against this parameter, 
the observations all lie upon a common curve. As anticipated, the internal 
friction is seen to have a maximum when this parameter is of the order of 
magnitude of unity. 

The maximum of the internal friction due to interciystalline thermal 
currents will be greater, the greater the elastic anisotropy of the individual 
crystallites. The effect of this elastic anisotropy enters through the factor 
R in Table 10. As a check upon the theoretical dependence of the magni¬ 
tude of the internal friction upon the anisotropy of the crystallites, a com¬ 
parison has been made of polycrystalline brass and aluminum. The in- 


93 


PHYSICAL INTERPRETATION OF ANELASTICITY 

tcrnal friction of the latter was found'® to be of an order of magnitude less 
than that of the former, in agreement with Table 10. 

Finally, experiments have been performed-" to check the theoretical pre¬ 
dictions concerning the precise manner in which the internal friction aris¬ 
ing from intercr>’stalline thermal currents varies with frequency of vibra¬ 
tion and with grain size. It was demonstrated in section A-III of this chap- 





Fio. 27.—Internal friction of polycrj-stalline brass, (.\fter Randall, Rose, and Zener) 


ter (pp. 80-81) that in the nearly adiabatic high-frequency case the ther¬ 
mal currents are confined to the immediate vicinity of the grain bound¬ 
aries and that in this case 


Internal friction'’^. 

The observed linearity of internal friction with in the nearly adiabatic 
case of large grain size and high frequency is demonstrated in Figure 28. 
On the other hand, in the opposite extreme case of nearly isothermal vi¬ 
bration, with small grain size and low frequency, the total amount of heat 


(1940) ^ Friction of Aluminum,” Phys. Rev., LVIII 

.1.Size," Trans. 
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which flows between adjacent grains per half-cycle is nearly independent 
of frequency. The thermal current is therefore proportional to the fre- 
(ILiency, and the rate of degradation of energy' by thermal flow is propor¬ 
tional to the square of the frequency. The mechanical energy dissipated 
per cycle and hence the internal friction are therefore in this extreme case 
proportional to the first power of the frequency. The experimental demon¬ 
stration of this deduction is given as Figure 29. 


Harmonic 



Fig. 28. —Example of (frequency) ^ dependence of internal friction in c.xtrcmc adiabatic 
case. (.After Zener and Randall.) 


We shall now return to the subject of the damping of longitudinal 
waves by the flow of heat from regions of compression to regions of dila¬ 
tion. It was pointed out that the optimum wave length, Xo, for maximum 
damping was extremely short and that therefore such damping is of no im¬ 
portance in artificially produced waves. The possibility arises, however, 
that this damping may be of importance to the Debye heat waves. This 
possibility is examined below. 

Shortly after introducing his now well-known theory' of specific heat. 
Debye^* pointed out that the thermal waves in nonconductors must be 

21. M. Planck, P. Debye, el al., Vortriige iiber die kinetiseke Theorie der Malerie und der 
Eleclrizitdt (Leipzig: Teubner, 1914). 
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coupled with each other in order that the thermal conductivity may be 
finite. If there were no such coupling, the thermal conductivity would be 
infinite. The coupling resides, of course, in the variation of frequency with 
dilation. No unique method exists, however, for using this coupling to 
find out how the waves interact. Debye introduced the viewjwint that the 
thermal waves give rise to fluctuations in density and that it is these dens¬ 
ity fluctuations which scatter the heat waves in a manner analogous to 

Harmonic 



Frequency, in kilocycles per second 

• 29.—Example of linear dependence of interna! friction upon frequency in extreme 

isothermal case. (After Zener and Randall.) ^ exiremc 

that by which fluctuations in the density of the air scatter light waves. 
R. Peierls22 has pointed out that this viewpoint leads to an intrinsic dififi- 
culty, namely, the scattering cross-section of the density fluctuations de¬ 
creases with the frequency of the waves as and therefore the mean free 
path would be proportional to u-*. Since the frequency density of the 
waves varies as only the integrand for the mean free path will involve 
i'--. The Debye vie^vpoint would therefore lead to an infinite thermal con¬ 
ductivity arising from the low-frequency thermal waves, 

22 . “ZurkinetischeTheoriederWarmleitunginKrista!le,”^««.d. Phys., Ill ( 1929 ), 1055 . 
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The above-mentioned difficulty in the Debye viewpoint disappears as 
soon as we recognize that the low-frequency, long-wave-length thermal 
waves may be damped by the heat currents which they generate. For fre¬ 
quencies much less than the optimum frequency, vq, the internal friction 
due to the thermoelastic coupling varies as v. The mean free path of these 
waves will therefore vary as The integrand for the average mean free 
path in the long-wave-length region, therefore, does not diverge. 


V. RELAXATION BY ATOillC DIFFUSION' 


In an unstressed solid solution, equilibrium is attained only when the 
solute atoms are randomly distributed. When stress fluctuations are pres¬ 
ent within a system, equilibrium conditions no longer correspond to a uni¬ 
form distribution of solute atoms. Thus, if an increase in concentration 
is attended by a lattice e.xpansion, then under equilibrium conditions a 
dilation will be attended by an increase in concentration. Fluctuations 
in dilation will therefore tend to induce fluctuations in concentration. 
In other words, fluctuations in dilation will induce fluctuations in the 
partial thermodynamic potential, Gk, of the solute k, fluctuations which 
are rela.xed only by atomic diffusion. That such diffusion will lead to an- 
elastic effects was first recognized by Gorsky,-^ who formulated the com¬ 
plete theory for the effects of this diffusion. No experimental study has as 
yet been made on anelastic effects resulting from atomic diffusion of either 
a macroscopic or a microscopic type. 

A characteristic feature of relaxation by atomic diffusion is the relative¬ 
ly long time of rela.xation and the rapid dependence of this time of relaxa¬ 
tion upon temperature. The relaxation time, r, will be given appro.ximate- 


ly by 



(159) 


where d is the mean distance over which diffusion must take place for rc- 
la.xation and D is the atomic diffusion coefficient. This diffusion coefficient 
may be written as ^ 

(160) 

where a is the distance that an atom moves in an elementary diffusion 
act and tq is the average time interval between elementary diffusion acts. 
Combining the above two equations, we obtain 

r^(^yr„. (161) 


23. W. S. Gorskv. “On the Transitions in the CuAu Alloy. III.” Phys. Zeitsebr. Sow., 
(1936), 77. 
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The rapid temperature variation of the relaxation time, r, comes about 
through the rapid temperature variation of tq. According to all theories of 
diffusion, the mean time interval, to, depends upon temperature essen¬ 
tially as exp (H. RT), where the quantity H is known as the “heat of ac¬ 
tivation” and R is the gas constant. According to the well-known Lang- 
muir-Dushman-^ theory, the coefficient is given by approximately {hN/II) 
where h is Planck’s constant and N is Avogadro’s number. Thus 

hN 

(162) 


The known values of the heats of activation of elements in metals are 
given in Table 12. As an example of the exceedingly long times required 
for an elementary diffusion act at room temperature, the following exam¬ 
ples will be quoted for the room temperature values of tq as derived from 
equation (102): 1.6, 2 X 0.5 X 10'^ seconds for // = 20.000, 30,000, 
40,000 cal/molc, respectively. 

The relative change in the elastic modulus associated with relaxation 
of all changes in the partial thermodynamic, Gt, is, according to equation 

(111a). 



Here wt is' the number of solute atoms of type k per unit volume. When the 
solute atoms form an interstitial solution, the partial thermodynamic po¬ 
tential will be defined in the standard manner, namely, as 



(164) 


where G is the free energy per unit volume. A different definition will be 
adopted in the case of solutes which form substitutional lattices, since in 
such cases diffusion of the solute atoms necessitates a counterdiffusion of 
the solvent atoms. Here the partial derivative in equation (164) implies 
the change in G induced by the addition of one solute atom with the simul¬ 
taneous subtraction of a solvent atom. After using the transformation, 

d e _ Be duk 
^k~'^k'dG,' 

24. S. Dushman and I. Langmuir. Phys. Rev.. XX (1922), 113. 



TABLE 12 


Heats of Activation of Atomic 
Diffusion Coefficients* 


Solvent 

Solute 

n (Cal/Molc) 


1 

fAu 

29,800 


1 

Cd 

22,350 

Ag. 

i 

■ Cu 

24,800 


1 

Sb 

21,500 


1 

Sn 

V 

21,400 

.\ 1 . , 

/Cu 

34.900 


iMg 

38,500 



:au 

51,000 

Au 

1 Cu 

27.400 


Pd 

37.400 


:ai 

40.200 


1 Be 

27.900 

Cu 

'Cu 

57.200 


'Si 

39.950 


(Zn 

42,000 

Fe (a). . 

' (Ct 

18,000 


m 

20,000 


Fe§ 

77,200 

Fe(7). , .. . 

/cii 

32,000 


IN 

34,000 


Fe§ 

48,000 


C 

25,000 


0 

29.000 


f. 

Au 

13,000 

Pb. 

Bi 

18.600 


: 

Hg 

19,000 

1 

Pb 

27,900 

< 

a 

5n 

24,000 


fn 

i 

VIo 

80,.500 

\v. . 

4 

4 

rh 

120,000 


1 

\ 

u 

100,000 


• Unless otherwise stated, data are from R. M. Barrer, Diffu¬ 
sion in and through Solids (New York: Macmillan Co., 19-11). 

t J. L. Snoek, Physica. VIII (1941). 71!. 

JT. S. Ke. h/flals Tech. (1943). (In Press.) 

S C. Birchcnall and R. Mehl, J. App. Phys., XIX (1948), 217. 

II C. Wells and R. Mehl. Trans AJ.if£.. CXL (1940), 294. 

ITT. S. Kfi. Phys. Rep. (1948). (In Press.) 
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and equation (164), we may re-write equation (163) as 



(1 65) 


In the case of a dilute solution or, more precisely, when the total free 
energ>^ minus the free energ>- of mixing is a linear function of the number 
of the various types of atoms, the denominator of equation (165) may lx* 
further reduced. In this case 

= kT } In iih - hi (.Y - m) 1 , (1 66) 

where .V is the total number of atoms per unit volume in the case of a sub¬ 
stitutional solute, the total number of interstitial positions in the case of 
interstitial solutions. Upon defining the atomic concentration, A', as 


we obtain 


and therefore 



dn- 


kT 

X 


A'-'(l - A')-' 


A'(l - A) EafdeV 
kT ■ X \dxl- 


(167) 


(168) 


The second factor is of the order of magnitude of the strain energy intro¬ 
duced by a single solute atom. 

For purposes of computation it is convenient to replace in erpiation 
(168) by ipiYo/lV), where p is the density, is Avogadro’s number, and 
IV is the atomic weight. Then, upon replacing A^^ by R, the gas constant, 
wc obtain 



= A(1 - A) 


IF£t/ 

pRT 



(169) 


We shall now estimate the maximum value of that might be expected. 
The product A(l —A) has the maximum value of -J. For the third factor 
we shall take copper as a typical solvent and shall consider 700° C. as a 
typical temperature, obtaining 210 for this factor. In estimating the maxi¬ 
mum value of the last factor, we combine the approximate relation 
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where u-. and ax are the atomic radii of the solute and solvent atoms, re 
spectively, with Hume-Rothery’s rule, 


and obtain 0.02 as the maximum value of the last factor. Combining the 
above approximations, we have 

-ii- < 1 . (170) 

The low atomic diffusion rates render it unlikely that examples will ever 
be observed in which the modulus is appreciably relaxed through macro¬ 
scopic atomic dilTusion. The relaxation occasioned by diffusion between 

adjacent grains is smaller than that for macroscopic diffusion, by the ratio 
R given in Table 10. 

In many cases the assumption of a dilute solution cannot be used. In 
such cases the relaxation strength is more conveniently written as 



Two extreme cases discussed below deserve special study. In one the cur¬ 
vature d'G dX‘ is much smaller than that given by the assumption of a 
dilute solution; in the other it is much larger. 

Certain pairs of metals are completely soluble in each other at elevated 
temperatures, while a limited region of immiscibility exists below a certain 
critical temperature. An example is given in Figure 30. Now the existence 
of a single phase throughout the entire concentration range implies that 
the curvature of G(X} is everywhere positive—in other words, that a given 
line can be tangent to G(X} at only one point. Such is the case in the lower 
curve of Figure 31. On the other hand, the existence of a two-phase region 
implies that a line can be drawn which is tangent to G(X) at two points, 
as illustrated by the upper curve of Figure 31. It is readily seen that im¬ 
miscibility implies that, over a limited concentration range, d'G dX- is 
negative. At the critical temperature it is therefore necessary that d'G dX‘ 
be zero at a certain concentration. For a critical concentration the relaxa¬ 
tion strength associated with atomic diffusion therefore approaches in¬ 
finity as the temperature is lowered to the critical temperature. Mayer-' 
has, in fact, argued that d-G/dX' will be zero over a finite range of con¬ 
centration and of temperature. 

25 W. Hume-Rothen-, G. W. Mabbot.and K. M. Chaitnel-Evans, Phil. Trans. Roy. Soc., .1, 
CCXXXIII (1934). 44.' 

26. W. Hume-Rothery, The Structure of Metals and Alloys (London, 1936), p. 59. 

27. J. Mayer and M. Alayer, Slaiislical Mechanics (New York: John Wiley & Sons, Inc., 
1940), pp. 312-13. 
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In the second extreme case the curvature d'G/dX'^ is much larger than 
would be given by the assumption of a dilute solution. This curvature is 
always anomalously large in those phases which have a narrow-solution 
range. As may be seen from Figure 32, a narrow-solution range at all tem- 



FlC. 30.—Example of a critical temperature for inimlscibility 



Fic. 31.—N'arialion of free energy near critical tem[>crature 

peratures necessarily implies a large curvature. Relaxation effects asso¬ 
ciated with atomic diffusion are therefore anomalously small in such 
phases. 

VI. RELAXATION' BY MAGNETIC DIFFUSION 

It has been recognized for over a century that a coupling exists between 
the magnetic and the mechanical properties of metals. Joule-® reported in 

"On the Effects of Magnetism upon the Dimensions of Iron and Steel 
Bars,” Phil. Mag., Vol. LXXVI (1847). 
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1847 that the magnetization of iron and steel resulted in changes in dimen¬ 
sions; and in the same year Matteuci-® reported that changes in dimensions 
resulted in a change in magnetization. Through this coupling all types of 
magnetic relaxation give rise to anelasticity. 

The type of magnetic relaxation most amenable to theoretical treat¬ 
ment is that which results from the macroscopic electrical eddy currents, 
which tend to shield the interior of a ferromagnetic material from changes 
in magnetic flux. Thus, suppose a tensile stress is suddenly applied to a 
very long and thin magnetized specimen, thereby suddenly changing the 
intensity of magnetization. This change in magnetization induces eddy 



Fig. 32.—Illustration of how a large curvature a*6*/a.V*is related to a narrow range of 
solubility. 


currents at the surface of the specimen of such a magnitude that the mag¬ 
netic flux within the specimen is unchanged. In other words, the magnetic- 
field intensity induced by the eddy currents is of just such a magnitude 
that the change in magnetic flux density, 

is exactly zero. The eddy currents gradually diffuse inward from the sur¬ 
face, thereby allowing magnetic flux to diffuse inward, at the same time 
reducing the magnetic-field strength to its original value. In such a case 
the magnetic-field strength, H, may be regarded as a thermodynamic po¬ 
tential, the magnetic flux density as a thermodynamic density in the sense 
of Section I of this chapter. 

The justification of the above viewpoint may be obtained directly from 


29. Mattcuci, Compi. reiid. Acad. Set., 1847. 



103 


PH\'SIC\L INTERPRETATION OF ANELASTICITY 

Maxwell’s equations. Neglecting the polarization current density com¬ 
pared with the true current density witliin the metal, one obtains 

dB 108 

^ = OTwi(172) 

where K is the electrical conductivit.v in ohm-' cm.-' When the stress, and 
therefore B and H, suffer small periodic fluctuations, we may set 

dB^= (C|irfo-„+Ci;(/o-j„+ . . .) , 

with analogous equations for dB^ and dB,. If the stress components are 
uniform over the specimen, as in tension, or vary, at most, in a linear man¬ 
ner, as in the torsion of a cylindrical rod or in bending, the substitution of 

these equations in equation (172) and the use of the condition that the 
divergence of B is zero lead to 

(173) 

where the magnetic diffusion coefficient. D, is given by 

^ = 0 74 ) 

According to this equation, instantaneous changes in stress leave the mag¬ 
netic flux density unchanged, the latter changing only by the diffusion of 

flux from the outside, the diffusion coefficient being given by equation 
(174). 

The change in energy per unit volume associated with a change in 
magnetic flux density \sH • (dB Ttt). Therefore, in the caseof a long, thin 
specimen, when H and B are parallel to the axis of the specimen, if H is re¬ 
garded as a thermodynamic potential, B/^tt must be regarded as the cor¬ 
responding density. Hence the equations relating strain, e, and magnetic 
flux density, B, to tensile stress, a, and magnetic-field strength, H, are 

de=E-^da-\-\dn , 
dB — 4T\d<T -p fidH . 

The strength of the relaxation due to magnetic flux diffusion, (Eb-Eu)/ 
En, is computed from these equations as 

As = . (175) 

The magnitude of the relaxation strength, A^, is a function of the state 
of magnetization. It is zero at zero magnetization, for at zero magnetiza¬ 
tion a stress can induce no change in magnetization. The relaxation 
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strength is also zero at saturation magnetization; for. at saturation, stress 
can likewise produce no change in magnetization. The maximum value 
of Aj.; will be in the vicinity of remanence magnetization. Becker and 
Ddring-'^® have shown that, at remanence. 

(176) 

(^i 


where X, is the tensile strain which accompanies a change in magnetiza¬ 
tion from zero to saturation, and (r, is a measure of the residual micro¬ 
stresses, which may have arisen either through prior plastic deformation 
or through the presence of impurity atoms. In pure well-annealed ferro¬ 
magnetic materials the only residual stresses arc those arising from mag¬ 
netostriction. In this case 


and therefore 


ffi ^ \Ei’ ; 


(177) 


No observations have been made on the magnitude of associated 
with the macroscopic diffusion of magnetic flux. Such observations would 
involve a comparison of the elastic moduli measured under conditions of 
essentially constant B, such as at high frequencies or upon thick speci¬ 
mens, with the values measured under conditions of essentiallv constant 
//, such as at low frequency or upon thin specimens. 

The effective time of relaxation, t. for magnetic diffusion is given ap¬ 
proximately by the equation 



( 178 ) 


where d is a. transverse dimension of the specimen, and the magnetic dif¬ 
fusion coefficient, D, is given by equation (174). In this equation, // refers 
to the reversible permeability, which ranges from 20 for the common 
steels to 14,000 for carefully prepared high-purity iron. Over this range of 
reversible permeability, D changes from 40 to 0.05 cm- sec. In the high- 
frequency dynamical methods commonly used to determine the elastic 
moduli of steels, the measurements are therefore made essentially under 
conditions of constant B. The internal friction arising from the macro¬ 
scopic diffusion of B is at a maximum when the period of vibration is com¬ 
parable to T, as given by equation (178). The precise manner in which the 
internal friction varies with frequency was first computed by Kersten^' 

30. R. Becker and \V. Doring, Ferromagnethmus (Berlin. 1939; New York: Edwards Bros., 
1943), p. 378. 

31. M. Kersten, “Zur Deutung der mechanischen Dampfung ferromagnetischen Werke- 
stoffe bei Magnetisierung.” Tech. Phys. XV (1934'l, 463. 



105 


PH\'SIC-VL IXTERPRETATIOX OF ANELASTICITV 


for the case of longitudinal vibration in specimens of circular cross-section 
and by the author^- for other examples. The formulae have been given in 
Section Ilia (pp. 81-83). 

Honda^^ discovered what he called, and what is now commonly known 
as, the A£-etTect. The term refers to the difTercnce in the elastic 

modulus of a specimen in the demagnetized state and at magnetic satura¬ 
tion. In the demagnetized state the elementary magnetic domains are free 
to change their direction of magnetization under the action of an applied 
stress and to grow at the expense of neighboring domains having a differ¬ 
ent orientation. At magnetic saturation such changes are no longer pos¬ 
sible. A complete theory of the A£-effect was developed by Akulov and 
Kondorsky.^^ Their formula is 

■^-.3 ^ ^ 

^ "S Xo^ioo . 


Here xu is the magnetic susceptibility at zero magnetization, Xioo is the 
longitudinal magnetostriction for a single crystal saturated in the (100| 
direction, and J, is the saturation magnetization. In their derivation of the 
above equation, Akulov and Kondorsky assumed constant stress through¬ 
out the specimen. Cooke and Brown^^ have subjected this formula to a 
careful experimental test. They find that their observed value of lE/E 
lies between that given by the above formula, in which constant stress was 
assumed, and a formula obtained under the assumption of constant strain 
throughout the specimen. review of the many e.xperiments ujjon the 
A£-effect has been given by Becker and Doring.^® They have shown that 
in soft magnetic specimens, where the only residual stresses are due to 
magnetostriction, A£, £ is of the order of magnitude of 0.3. Values as high 
as 0.17 have been observed in nickel.*' 

As previously mentioned, relaxation by macroscopic diffusion of mag¬ 
netic flux does not occur either in the demagnetized slate or at magnetic 
saturation. Hence such diffusion does not contribute to the A£-effect. If 
however, the elastic modulus is measured as a function of magnetization, 
from zero magnetization to saturation, the macroscopic diffusion of mag- 

1010^ ^ Zener, “General Theory of Macroscopic Eddy Currents," Phys. Rev., LIII (1938), 

^ i Kusakatee, "Veranderung des Elastizitatskoeffizienten 

380 Substanzen mfolge von Magnclisierung," Phys. Zeitschr., Ill (1901-2), 


34. Akulov and Kondorsky. Zeitschr. f. Phys. LXXVHI (1932), 801; LXXX\' (1933), 661 

35. W. T. Cooke and \V. F. Brown, “Variation of the Internal Friction and Elastic Pnn 

slants with Magnetization in Iron,” Phys. Rev., L (1936), 1158, 1165. ^ ” 

36. Op. cit. 

xciv (1935 )“,’ Ph^-. 
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netic flux will contribute to changes in the modulus at some intermediate 
range of magnetization. 


VII. RELAXATION* OF ORDERED DISTRIBUTION'S 


nie cr\ stals, such as ionic crystals, the position of each type of 
atom or ion is uniquely determined. Thus in CsCI, the Cs ions are located 
at the lattice points of one simple cubic lattice, the Cl ions are at the lattice 
j)oints of a second simple cubic lattice. In other crystals,of which j3-brass of 
the composition CuZn is an example, each type of atom manifests a prefer¬ 
ence for certain lattice positions below a certain critical temperature. Thus 
below 470° C. the copper atoms in CuZn have a preference for the lattice 
positions of one simple cubic lattice, the zinc atoms have a preference for 
the lattice positions of the othersimple cubic lattice. Such solids are said to 
be “ordered” below their critical temperature. Further examples of ordered 
metals, as well as a review of their properties, have been given in several 
recent reviews.^^ Ordered structures have two properties which render 
them of special interest to the present discussion: (1) below the critical 
temperature the degree of order is incomplete, the equilibrium order be¬ 
coming more complete as the temperature is lowered; and (2) a change in 
lattice dimensions accompanies changes in degree of ordering. As was first 
pointed out by Gorskya change in lattice dimensions by an externally 
applied stress will therefore change the equilibrium order, and the time 
delay associated with the establishment of equilibrium will lead to an- 

4 

elastic effects. 

The quantitative analysis of the anelasticity associated with ordering 
may most conveniently be approached by considering the entropy of a 
unit volume to consist of two parts. The first part pertains to the uncer¬ 
tainty in the positions of the individual atoms with regard to their respec¬ 
tive equilibrium positions and will be referred to as the “lattice entropy,” 
Sl- The second part pertains to the uncertainty in the equilibrium posi¬ 
tions of the atoms and will be denoted by 5.u. This second part is common¬ 
ly referred to as the “entropy of mixing.”^® Corresponding to these two 
entropies, there are two temperatures, Tl and T,/, defined as follows: 



38. Borelius, Zeiischr.f. Klectrochemie, XL\' (1939). 16; and F. C. Xi-t and W. Shockley. 
“Order-Disorder Transformations in .Alloys,” Rev. Mod. Phys., X (1938). 1. 

39. W. S. Gorsky, “Die elastische Nachwirkung in geordneler CuAu Legierung,” Phys. 
Zeitschr. 5ow.. VIII (1935). 562. 

40. R. H. Fowler and E. Guggenheim. Statistical Thermodynamics (Cambridge; Cam¬ 
bridge University Press, 1939), p. 163. 
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where H is the heat content function per unit volume. The entropy of 
mixing, can increase or decrease only through a transfer of heat from 
the vibrational coordinates to the coordinates of mixing or vice versa. The 
anelasticity associated with ordering may therefore be regarded as arisin 
from the flow of heat between the vibrational coordinates and the mixin 
coordinates. In computing the rela.xation strength associated with order¬ 
ing, due regard must therefore be given to the conditions which sat¬ 
isfies. 


Consider, as a first example, that 7T remains essentially constant, either 
through maintenance of thermal equilibrium with the surroundings or 
through maintenance of thermal equilibrium between different parts of 
the specimen, as in transverse flexure. Then, from equation (105). 




(179) 

where 

fdS.,r\ 

varuX 

(180) 

and 



/ d€\ 



II 

< 

(181) 


An alternative expression is obtained from equation (105a), namely, 



(1K2) 


Those parts of the thermal expansion and of the heat content which are 
associated with disordering vary with temperature in approximately the 
same manner. An example is given in Figures 35 and 34 for the alloy 
CuaAu, the data for the specific heat being taken from Sykes and Jones" 
and those for the thermal expansion coefficient from Nix and MacNair.^- 
The parallel behavior of the anomaly in thermal expansion and in heat con¬ 
tent suggest that the derivative {de/dSAi)^ is essentially independent of 
temperature. We shall therefore assume 



where ^e is the strain associated with the order-disorder transformation 
and ASm is the corresponding change in entropy. Thus 




(183) 


^ ' “The Atomic Rearrangement Process in the CooDcr-Gold 
Alloy Cu,Au.” Proc. Roy. Soc., London, A, CLVII (1936); 213. ‘ ^ 

42. K C. Nix and D. MacNair, “A Dilatomctric Study of the Order-Disorder Transforms 
tion m CuAu Alloys,” Pkys. Rev., LX (1941), 320. i^Jsoraer iranstorma- 


tc tc 
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The entropy change. SSm, may be readily estimated from theoretical con¬ 
siderations. Thus let the alloy have the composition ArB,. If a unit volume 
contains atoms, the entropy of mixing in the disordered state, regard¬ 
ing all atoms of a like kind as indistinguishable, is klnlNl /{N—pN)\ 
qX ):!, where p and q are the ratios r/{r+s) and y/(r+5), respective¬ 
ly. In the completely ordered state the entropy of mixing is zero. One 
therefore obtains 


Thus 


A5.U = — I pin p-\- gin q]. 


l0.64.V^,{r, s) = (1, 1) ,j 
/0.56A'^^,(r, s) = (1,3).) 


(184) 

(18:>; 



Equations (183) and (185) will now be applied to the case of CuaAu. For 
this case Ac is 0.0019,'*^ and. from equation (185) ASyf is 0.56 Xk. We find 

Ab = 0.004, when hm = • 

Just below the critical temperature, is at least thirty times as large as 
Ml- Thus just below the critical temperature Ag rises to a value as large 
as 0.1. 

As previously mentioned, equation (183) is applicable to those cases in 
which the temperature of the lattice vibrations, Tg, remains essentially 
constant. These conditions are most readily obtained by flexure of thin 
specimens, where the fluctuations in Tg are rapidly equalized by trans¬ 
verse thermal conduction. It is, of course, necessaiy that the reIa.xation 
time for equalization of Tl across the specimen be short compared with 


43. Ibid. 
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the relaxation time for the equalization of and of 7'/.. Elastic after¬ 
effects under these conditions have been reported by Gorsky^' for CuAu. 
Gorsky did not find that this alloy behaved as if it had a single time f re¬ 
laxation. and he believed a second type of relaxation to be present. Since 
CuAu forms a tetragonal lattice in the ordered state, with a tetragonalitv 
ratio differing only slightly from unity, it is possible that this second type 
of relaxation is associated with the movement of twin interfaces. 

In many cases the ordered structure has a tetragonal or a rhombo- 
hedral symmetry-, although the disordered state has a cubic symmetry. 

The alloy CuAu is an example. In these cases Ac’ in equation (183) must 
be replaced by 


(Ac)= = 

' ' rtV 


1 

5 


Ac^ 


+§-i 



Fig. 34.—Thermal expansion cocfticicnl of CuaAu. (.\flcr Xi.x and Mac.Xair) 


where Acn and Ac^ are the strains introduced by ordering parallel and 
normal, respectively, to the principal axis. In the case of CuAu, these 
strains are 0.07 and 0.035, respectively.^^ The quantities (Ac)|v and there¬ 
fore Ajs are hence two orders of magnitude greater than in the case of the 
cubic-ordered alloy, CusAu. 

The preceding discussion has been limited to the case in which the tem¬ 
perature associated with the lattice vibrations, Tz., remains essentially 
constant. In this case A^ refers to the complete relaxation of all fluctua¬ 
tions in the temperature associated with ordering, 71,,. We shall now con¬ 
sider the case in which no heat flows in or out of an element of volume. 
The only exchange of heat will then be between the mixing coordinate and 
the lattice vibrations. In order to compute the associated relaxation 

44. Gorsky, Phys. Zeitsclir. Sow., VUI (1935), 562. 

45. Ibid. 
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strength, we start from the following equations, which relate the depend¬ 
ent variables e, 5^, and Sl with the independent variables tr, and Ti^\ 

de=E-^da-\r\^.dT^t + \LdT^ , | 

dSyf = \^ida -{- fi^/dTM , !> ( 186 ) 

dS L ~ 'Ktdff ^i,dT /, . j 


In writing these equations it has been assumed that r,, is of 

negligible magnitude. This assumption seems reasonable, since Tl can in¬ 
fluence Ssf only through the thermal e.xpansion associated with a change 
in T We are now interested in the relative difference in Young’s modulus 
in the two cases in which S^t and Sl are both held constant and in which 
the sum of 5,/ and Sl is held constant and Tm and Tl are kept equal to 
one another by a flow of heat between the mixing coordinates and the lat¬ 
tice vibrations. In other words, we shall now compute 








\de/s:^f+SL’ t^-Tl 


A straightfonvard computation based upon the set of equations (186) 
leads to 

l(^) -(#) (■ ('»') 


Tm 


In order to see the effect upon of the restraint of Tl, equation (187) 
should be compared with equation (182). When y.M = mlj the second 
factors differ by a factor of 2. At least in the case of CujAu, the third factor 
in equation (187) is two orders of magnitude less than the corresponding 
factor in equation (182). This result is a consequence of the close identity 
of the two terms in the third factor of equation (187). In order to interpret 
this result, we observe that 

\d5.wA. T, V Aw- Tl V 


( 188 ) 


•5.V ■ 5 L 


which follows directly from equation (186). The close identity of the two 
terms in the last factor of equation (187) therefore implies that when a ten¬ 
sile stress is applied to the metal under strictly adiabatic conditions with 
respect to all types of heat transfer, Tm and Tl are increased by essen¬ 
tially the same amount. 

Observations have been made by Koster^® upon the internal friction of 

46. W. Roster, Zeitsckr.f. Elekfrochemie, XLV (lOSQ). 31. 
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/3-brass under conditions where + 5^. was essentially constant. His ob¬ 
servations show a maximum near the critical temperature. One must sur¬ 
mise from these results that the two derivatives de, dS,^f and de/dS^ in 

equation (187) in no wise cancel each other in /3-brass as they do in 
CusAu. 

The only study which has been made of the time of relaxation for the 
establishment of equilibrium order is that of Gorsky*' upon CuAu. He 
finds this time of relaxation to obey the relation 

7 ^ cQ/fiT , 

where the heat of activation, Q, has the same value as for macroscopic dif¬ 
fusion, namely, 27,400 cal mole. At 275° C., which is about 100° C. be¬ 
low the critical temperature, this time of relaxation is 10 seconds. 


VIII. RELAXATION OF PREFERENTUL DISTRIBU¬ 
TIONS INDUCED BY STRESS 


a) GENERAL THEORY 

The atoms in a disordered solid solution of cubic structure are distrib¬ 
uted isotropically in the absence of an externally applied stress. Thus, in a 
stress-free solid solution the distribution of solute atoms is isotropic with 
respect to each individual solute atom, and the nearest solute atom has an 
equal probability of lying along any of the crystallographically allowed di¬ 
rections. In body-centered cubic (b.c.c.) structures the interstitial posi¬ 
tions, namely, the centers of the faces and edges, have tetragonal sym¬ 
metry, the tetragonal axis pointing along one of the three principal axes. 
In a stress-free b.c.c. lattice containing interstitial solute atoms, these 
atoms are distributed isotropically, i.e., the principal axis of the position 
occupied by a particular solute atom has an equal probability of being 
parallel to any one of the three principal axes. 

The equilibrium distribution of solute atoms is no longer necessarily iso¬ 
tropic in the presence of a stress. Thus in the case of a face-centered cubic 
(f.c.c.) lattice the axis joining nearest neighbors is along a < iI0> di¬ 
rection. If the solute atoms are larger than the solvent atoms, two neigh¬ 
boring solute atoms whose axis is along the [110] direction will produce 
an anisotropic distortion, in the sense that the lattice will be locally ex¬ 
tended along the [110] direction more than along the directions normal 
thereto. A tensile stress along the [110] direction will therefore induce a 
preferential distribution in which more pairs of solute atoms have their 


Cutb-^‘-I'kris.allcn von 
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axes lying along the [110] axis than along any of the other five < n0> 
axes. Similarly, a tensile stress along the [100] axis of a b.c.c. lattice with 
interstitial solute atoms will induce a preferential equilibrium distribution 
in which each solute atom has a greater probability of being in an intersti¬ 
tial position whose tetragonal axis is along the [100] axis than along 
either of the other two < 100 > axes. The continual striving of the metal 
to maintain an equilibrium distribution will give rise to anelasticity. In 
particular, the associated internal friction will be a maximum when the 
time of relaxation for the establishment of equilibrium is comparable 
with the period of vibration. 

W e shall now compute the relaxation strength associated with the es¬ 
tablishment of the equilibrium preferential distribution. Toward this end 
we consider that a tensile stress, u. is applied along one of the crystallo¬ 
graphic axes which may be preferred, e.g., along a < I10> axis in a f.c.c. 
lattice. We shall denote by n the number of elements per unit volume 
which may acquire a preferential orientation, e.g., number of pairs of 
solute atoms in the case of substitutional solid solutions, number of solute 
atoms in the case of interstitial solution in b.c.c. lattices. The quantity tip 
will denote the excess number of elements, per unit volume, in the pre¬ 
ferred orientation over the number in this orientation when the distribu¬ 
tion is random. 

The tensile strain, €. and the energy difference, U, of an element in a 
preferred and a nonpreferred orientation may then be considered as func¬ 
tions of <r and Up. Thus 



Under ecjuilibrium conditions. will be proportional to U ; thus 



(189) 


where is a numerical coefficient whose precise value will have to be eval¬ 
uated for each particular case from statistical considerations. The relaxa¬ 
tion strength may now be obtained through elimination of tip and U from 
the above three equations. We obtain 





(190) 
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An estimate may usually be made of the derivative (de dnX from ex¬ 
perimental data, while no such estimate is available for the derivatives 
(dU/ and {dU dnp\. It is therefore desirable to transform these last 
two derivatives. The first is transformed by forming the perfect differ¬ 
ential, 

(/(£ —ter) = —6(Ia-\-UdUp, 

where E is the energ>^ per unit volume. From this equation we obtain 



(191) 


The second unknown derivative. (dU will be estimated from the 

following considerations. It is anticipated that both <7 and influence C 
primarily through their effect upon the tensile strain, e. We therefore ex¬ 
pect the change in U induced by diip to be of the same order of magnitude 
as the change in (j induced by a stress increment Sa, which gives rise to 
the same strain increment, Se, as does 6«p. We therefore set 



(192) 


where uka. numerical coefficient of the order of magnitude of unity. This 
numerical coefficient corresponds to the Weiss factor in the theory of the 
paramagnetism of ferromagnetic materials. 

Upon substituting equations (191) and (192) in equation (190), we ob¬ 
tain 

T 

where 

This equation for Ae is, of course, valid only as long as T is above the criti¬ 
cal temperature, T^, for a self-induced preferential distribution. The rela¬ 
tion of r, to V and To must be obtained for each particular case through 
statistical mechanics. 


) 


(194) 


b) FXAMPLE OF SUBSTITUTIONAL SOLID SOLUTION 

A relaxation has been observed in a-brass by the author,^* for which he 
had no interpretation at the time that the measurements were made. The 
internal friction of an a-brass single crystal vibrating at a frequency of 
about 600 cycles/sec was found to rise steadily as the temperature was 

C1943).S‘22^'"‘'‘'' Crystal.- Trans. CLII 
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increased to 425° C., and then to decrease as the temperature was further 
increased. The observations are reproduced in Figure 35. Since zinc atoms 
are larger than copper atoms, it appears ver>" likely that this anomalous 
behavior of the a-brass cr>'stal was due to the anelasticity associated with 



Fig. 35.—Internal friction due to stress induced by preferential orientation of axes joining 
pairs of solute atoms. (After Zener.) 

the relaxation of the preferential distribution induced by the stresses of 
vibration. The maximum internal friction occurred at that temperature 
at which the time of relaxation of the preferential distribution was equal 
to the period of vibration. This viewpoint is further strengthened by the 
fact that the observations agree very well with a formula for a single time 
of relaxation with a heat of activation, namely, 

tan 5 = 


(195) 




where 
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P = -L 
To 


(196) 


and r is the period of vibration. The full curve in Figure 35 is the plot of 
equation (195). The estimated heat of activation //, 33.000 cal/mole, is 
consistent with the heat of activation for the dilTusion of zinc in copper ob¬ 
served by Rhines and Mehl.^« This heat of activation they found to de¬ 
crease from 40,000 to 31,000 cal mole as the zinc concentration increased 
from 10 to greater than 16 per cent. 

In a later paper'*^ the author has shown that the above interpretation is 
consistent with the observed rela.xation strength of the a-brass crystal, 
namely, 0.025 at 400° C. Since the full relaxation curve was observed at 
only one frequency, it is not possible from these e.xperiments to reach any 
conclusion regarding the value of the critical temperature, 

It is expected that the same type of anomalous anelastic behavior found 
for a-brass crystals will be found in all solid-solution crystals in which the 
solute and solvent atoms have different diameters, as is usually the case. 
Particularly striking crystallographic effects will be found in body-cen¬ 
tered crystals. Here the axes joining closest neighbors will be along the 
< 111 > directions. The angles which each of these directions makes with 
a < 100> axis are identical. Hence a tensile stre.ss api)Iied along a < 100> 
direction will not induce a preferential orientation of pairs of solute atoms 

and hence will not be associated with the above-described type of anelas- 
ticity. 


c) EXAMPLE OF INTERSTITIAL SOLUTION 


The first recorded anomaly in the internal friction of iron as a function 
of temperature was given by Cantone^* in 1895. During the ensuing half- 
century this anomaly was reported by many investigators, in apparent 
Ignorance of the previous work. In Figure 36 is reproduced the observa¬ 
tion of Woodruff^- (1903) upon the influence of temperature on the time 
of decay of the vibrations in steel tuning forks. Since the time of decay is 
inversely proportional to the internal friction, the time of decay has a 
minimum where the internal friction has a maximum. Woodruff’s ob¬ 
servations were the first to demonstrate that the temperature at which the 
internal friction is a maximum depends upon the frequency, being higher, 


51. M. Cantone, Nttovo cimento, I (1895), 165, 205; IV (1896), 270, 354. 

0903)1*321"'^ Temperature upon a Tuning Fork," Phys. 
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the higher the frequency. A very extensive study was made by Robin^^ in 
1910-11, typical results being reproduced in Figure 37. Robin was the 
first to demonstrate that the anomalous internal friction was greater in 
low- than in high-carbon steel and that the anomaly disappears upon 
plastic deformation, reappearing again after annealing. Several later au¬ 
thors have reported less detailed observations upon the anomaly.'"^ 

Mild steel also shows an anomalous magnetic aftereffect which cannot 
be understood from magnetic considerations alone, as was first demon- 



Temperoture ‘C 

Fig. 36. —Internal friction of steel tuning forks. (After Woodruff.) / « 128 c>’cles‘'sec 

strated by Ewing (1885).“ In careful observations upon both the magnetic 
and the mechanical aftereffect in carbonyl iron, Richter,^''’ demonstrated 
the very close relation between the two effects. An example of this relation 
is illustrated in Figure 38. Richter demonstrated that the magnetic and 
mechanical aftereffects have the same heat of activation, namely, 19,600 
cal/mole, and that the magnitude of both aftereffects has a temperature 
coefficient of —0.003 near room temperature, a temperature coefficient 


53. F. Robin, “Report on the Wear of Steels and on Their Resistance to Crushing," Iron 
and Steel Inst., Carnegie Schol. Mem., II (1910), 259; “Ph6nomene de [’extinction de son dans 
le fer,” Cmnpt. rend., Acad. Sci., CL (1910), 780; and “The \-ariations in the .\coust.cal Prop- 
erties of Steel with Changes in Temperature/^ Iron on4 Sl<€l Inst, CarnegU Sckcl. .l/w., Ill 
(1911), 125. 

54. K. Jokibe and S. Sakui, Phil. Mag., XLII (1921), 397; A. Esau and E. \o\zi.Zeitschr.f. 
tech. Phys., XI (1930), 78. 

55. J. Ewing, Phil. Trans. Roy. Soc., London, CLXX\T (1885), 554. 

56. G. Richter. “Cber die mechanische und magnetische Xachwirkung des Carbonylei- 
sens,” Ann. d. Phys., XXXII (1938), 683; and chapter on “MagneU^he und mechanische 
Nachwirkung,” in R. Becker (ed.), Probleme der technischen Magneiisierungskurve (Berlin. 
Springer, 1938; New York: Edwards Bros., 1944). 
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which is consistent with a variation as 1 T. Of great signihcance was Rich¬ 
ter’s observation that a strong magnetic field suppressed the magnetic 
aftereffects but had no influence upon the mechanical aftereffects. It was 
therefore evident that their common origin was of a mechanical rather 
than of a magnetic nature. This conclusion invalidated Snoek’s original 
hypothesis*' that the elastic aftereffect had its origin in magnetostriction 
and in the concentration of solute atoms along the boundaries of the mag¬ 
netic domains. The essential role that interstitial solute atoms play in the 
aftereffects was found by Snoek^'* when he removed all traces of carbon 
and of nitrogen, thereby removing both the magnetic and the elastic after¬ 



effects. When less than 0.01 per cent of either carbon or nitrogen was re- 
I>laced, the aftereffects returned. Two years later he^^ developed the the¬ 
ory of the preferential distribution of interstitial carbon and nitrogen 
atoms. 


In his theory Snoek assumes that the interstitial positions occupied by 
iron and nitrogen atoms are at the centers of the cube faces and at the 
centers of the cube edges, which are crystallographically equivalent. 'Phis 
assumption as to the location of the dissolved carbon and nitrogen atoms 
is consistent with X-ray evidence,*® but is inconsistent with the viewqjoint 
that interstitial solute atoms must be located at the largest interstitial 

57. J. Snoek, “Time Effects in Magnetisation,” Physica, V (1938i, 663. 

591^^' ■*' "^^lechanical After-effect and Chemical Constitution,” Physica, VI (1939;, 


Parser, H. Li,,son and Andre.v M. 











119 


PH\^S1CAL INTERPRETATION OF ANELASTICITV 

holes. Thus the interstitial position at 0) corresponds to a consider¬ 
ably larger hole than the position (i, 0, 0) assumed by Snoek. When the 
criterion of minimum strain energy is employed to decide the location of 
the interstitial solute atoms, it is found that the (A. 0, 0) type positions are 
indeed correct. Although the iron lattice is cubic, the local symmetry of 
these interstitial positions is tetragonal. Thus the interstitial position 
(4, 0, 0) is closer to the two neighboring carbon atoms along the [100] axis 
than to the four other neighboring carbon atoms. The (100| direction may 
therefore be regarded as the tetragonal axis of the (^, 0. 0) interstitial posi¬ 
tion. Similarly, every interstitial position has associated with it a tetrago¬ 
nal axis which is parallel to one of the three principal < 100 > axes. When 
an iron crystal is under no stress, all interstitial positions are equivalent, 
and the interstitial solute atoms are therefore randomly distributed. When, 
however, a tensile stress is applied along a principal axis, the equilibrium 
distribution will be one in which those interstitial positions are preferred 
whose tetragonal axes are parallel to the tensile axis, since an interstitial 
solute atom distorts the surrounding lattice in an asymmetrical manner, 
causing a greater dilation along its tetragonal axis than normal to this 
axis. Conversely, a tensile stress along the tetragonal axis will tend to in¬ 
duce a preferential distribution of interstitial solute atoms in those inter¬ 
stitial positions with tetragonal axes parallel to the tensile axis. 

The Snoek theory of the anomalous behavior of mild steel is thus a spe¬ 
cial case of the general theory developed in section a (pp. 111-13). It is 
therefore anticipated that the relaxation strength will be given by equa¬ 
tions (193) and (194). In applying these equations, we shall initially as¬ 
sume that the tensile stress is along the [100] axis. When we take cogni¬ 
zance of the fact that there are twice as many nonpreferred as preferred 
positions, we find the numerical coefficient defined in equation (189), to 
be (2/9). It remains to evaluate To. This quantity may be obtained from 
the data relating changes in the lattice parameters of martensite with car¬ 
bon composition.®^ From these data one deduces that 



where 5Cp is the change in concentration, referred to total number of iron 
atoms, of the carbon atoms in the preferred interstitial positions, — 5C 
the change in the nonpreferred positions. Using these data, one obtains 

To= 1190X, , 

HiU Book Co"^,l93^V2^2' ^ ' McGraw- 
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where A', is the concentration of carbon atoms, in weight per cent. For 
concentrations of carbon obtainable in solid solution in mild steel, is 
very small, less than 50° K. Equation (193) thus reduces to 

, IIQOX, 

=- - -. ( 198 ) 

Aa an example, for iron with 0.01 weight per cent of carbon, the above 
formula gives the value of at room temperature as 0.040. The experi¬ 
mental value of Afi- for these conditions has been found by Dijkstra and 
reported by Polder^^* as 0.043. The verification of the temperature depend¬ 
ence of At- as given by equation (198) is contained in the work of Richter.’’^ 
The most convincing verification of the Snoek theory has been given by 
Dijkstra through his observations upon the crystallographic dependence 
of the relaxation strength in single crystals. If a tensile stress is ap¬ 
plied along a <111> axis, no interstitial positions become preferred, 
since all the tetragonal axes are equally inclined toward the tensile axis. 
One should therefore anticipate no relaxation due to carbon or nitrogen 
atoms. This has been shown by Dijkstra to be indeed the case. 

It would be anticipated that the anelastic effects associated with carbon 
and nitrogen in iron would be associated with all interstitial solutions in 
b.c.c. metals, provided that the solute atoms are sufficiently large to pre¬ 
fer the (j, 0, 0) type positions over the (^, j. 0) type positions. Such effects 
have been described in detail by Ke®’*“ for the case of carbon, nitrogen, 
and oxygen dissolved in tantalum, and have likewise been observed in 
columbium'^^^ and in tungsten.®^*^ 

Recent observations by Dijkstra®^ indicate the potential usefulness of 
anelastic measurements in studying the metallurgy- of interstitial solutions 
in b.c.c. metals. Since the relaxation strength, or any quantity related 
thereto, such as the internal friction at the optimum temperature, is a 
measure of the concentration of solute atoms in solid solution, its measure¬ 
ment may be used to follow the rate at which the solute atoms come out 
of solution and, finally, even to measure the solution pressure in equilib- 

62. D. Polder, “Theory of the Elastic .After-effect and the Diffusion of Carbon in .Alpha- 
Iron,” Philips Research Repts., I (1945), 1. 

63. Loc. cil. 

64. L. J. Dijkstra. “Elastic Relaxation and Some Other Projjerties of the Solution of 
Carbon and Nitrogen in Iron,” Philips Research Repts., II (1947), 357. 

64a. T. S. K$, “Internal Friction in the Interstitial Solid Solutions of C and 0 in Tanta¬ 
lum,” Phys. Rev., 1948. 

64b. L. J. Dijkstra, personal communication. 

64c. J. Snoek, personal communication. 

65. L. J. Dijkstra, A.I.M.E. (in press) (1948). 
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rium with the precipitate. In such a study of nitrogen in iron Dijkstra has 
not only measured with a precision heretofore unattainable the solution 
pressure of iron nitride from the eutectoid temperature down to 200° C. 
but has also discovered a new phase and has measured its solution pres¬ 
sure, a phase which is metastable with respect to iron nitride. The exist¬ 
ence of such a metastable phase is consistent with contemporary dilato- 
metric measurements by Cohen and Averbach.®^* 

In our prior analysis it has been assumed that the applied stresses arc 
so low that the degree of preferential distribution is a linear function of 
stress. This restriction will now be removed. If « is the total number of 
interstitial atoms per unit volume and Hp is the number in the interstitial 
positions with tetragonal axes parallel to the local tensile axes, we may 
define the order parameter, P, by 


When defined in this manner, P varies from 0 to 1 as the distribution 
changes from complete randomness to preferred positions for all the inter¬ 
stitial atoms. If (J denotes the decrease in potential energy which an in¬ 
terstitial atom suffers in changing from a nonpreferred to a preferred in¬ 
terstitial position, then it may readily be deduced that 


where 


P =: — _ \ 



( 200 ) 

( 201 ) 


It now remains to deduce how the energy change, i/, depends on the stress. 
We assume that U is independent of temperature. Then, at the absolute 
zero temperature, 

d{U - ea) ^ - eda-UdUp , (202) 


where V is the energy per unit volume. Since the left-hand side is a per¬ 
fect differential, 



(203) 


If we now assume that 1/ is a linear function of stress, we obtain 

0/ = ^ <j . 

N dCp 

65a. M. Cohen and B. L. Averbach, A.I.M.E. (1948; (in press). 


( 204 ) 
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In the case of carbon, we obtain, in view of equation (197), 

VE 

RT 

where e is the elastic strain. At room temperature this equation becomes 

.V = 3 8 6 e . 

In Figure 39 the degree of order, P, is plotted as a function of x and of e. 
It is seen that a large reduction of the anelastic effects by residual stresses 


Etestic Stroin et R.T. 



Fig. 39.—Degree of order of carbon in iron induced by strain 

would require these stresses to give an elastic strain of about 0.01, that is, 
stresses of the magnitude of 300,000 psi. 

It has previously been mentioned that in mild steel the concentration of 
carbon in solid solution is always so small that Tq and hence also probably 
Tr are less than 50° K. In mild steel we may therefore neglect the interac¬ 
tion of the carbon atoms upon one another. The author®^ has previously 
suggested that it is just this interaction which gives rise to the tetragonal- 
ity of quenched steel, a body-centered structure in which all the carbon 
has been retained in solution, a structure known as “martensite.” It looks 
as though in martensite the critical temperature for self-induced prefer¬ 
ential distribution, Tc, is above room temperature. It is therefore of impor- 

66. C. Zener, “Kinetics of Decomposition of Austenite," Trans. AJ.AI.E., CLXYII (1946), 

sso. 
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tance to obtain the relation between and p\nd Toward this end we 
write the expression for the free energ>'. F, of a unit volume of unstressed 
iron. Upon letting M denote the total number of interstitial positions, n 
the number of carbon atoms, the number in the preferred type in posi¬ 
tions, and, finally, Xtip the excess of tip over n 3, we obtain 




dltr. 


(A«„) - kTlu j - , 

' n„\ (U — ;/„)! 




(it - lip) \{2M - It d-llr,) 


The equilibrium value of Slip is determined by the condition 



When, as is always the case in actual specimens, the total number of car¬ 
bon atoms is very small compared with the number of iron atoms, i.e.. 
H if/ « 1. this equilibrium condition reduces to 


where 




.r = kT III 




(205) 


(206) 


When T is sufficiently large, equation (205) has no solution other than 
X = 0. As T is decreased, this equation attains two additional roots. The 
critical temperature, T,, is determined by the condition that at the 
transition from the random distribution to the preferred distribution is 
attended by no change in free energy. This equilibrium condition is satis¬ 
fied when the integral of the left member of equation (205) with respect 
to X from X = 0 to the third root is equal to the corresponding integral of 
the right member. This root cannot be found graphically merely by plot¬ 
ting the right member of equation (205). As shown by Figure 40, such a 
plot does not have sufficient dispersion in the region of interest. Upon tak¬ 
ing cognizance of the fact that areas remain unaltered by simple shears, 
we see that the requisite dispersion can be obtained by a suitable simple 
shear parallel to the vertical axis. Such a sheared plot is shown in Figure 
41, from which plot a graphical solution for Xc gives 

Xf = 0.34 , 



— It 


(207) 
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Upon combining equations (192), (194). and (207) and upon using our 
previously found value of 5 for we obtain 

T, = l.()5rr„. (208) 

Upon returning to equation (193). we see that the maximum theoretical 
value of the relaxation strength, namely, the value at 1 \, is given by 

(A) 
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Shear Stress 

Fic. 42.—Illustration of movement of twin interface. Twin interface is a (Oil) plane of 
tetragonal lattice. * 

Below the critical temperature, relaxation will come primarily through 
the movement of interfaces between regions in which the preferred direc¬ 
tions for the tetragonal interstitial positions are along different principal 
axes. A possible example of such an interface is illustrated in Figure 42, 
where the interface is in a (Oil) plane. Residual stresses will impede the 
movement of these interfaces in a manner analogous to that in which they 
impede the movement of the boundaries of magnetic domains. 

The above interpretation of martensite as an ordered structure is based 
upon the assumption of an interaction between the solute carbon atoms of 
sufficient magnitude to induce an ordered distribution. A detailed study 
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b\- has showTi that such interaction is extremely small in the case of 
oxygen dissolved interstitially in tantalum. 

B. INHOMOGENEOUS REL.\X.\TION 

I. CON'CEPT OF TWO-COMPON'ENT SYSTEMS 

In the previous section we considered that in every' element of the speci¬ 
men stress can relax through the relaxation of the fluctuations of some ther¬ 
modynamic potential, such as temperature. It has long been recognized®^ 
and has been repeatedly discussed®® that anelastic phenomena may also 
arise through a mixture of two phases, one of a truly viscous nature, in the 
sense that therein 

shear strain) =6'“' {shear stress) -f 7 “' shear stressJ, (209) 

and one of a truly elastic nature, in the sense that 

Strain-^Stress. (210) 

The coexistence of both types of regions is necessary' in order to interpret 
the observed mechanical behavior of metals. If only the first type of media 
were present, a specimen would behave like a supercooled fluid, creeping 
but manifesting no creep recovery. 

The qualitative behavior of a two-component system may readily be 
understood by reference to Figure 43. The specimen with no load and in a 
completely relaxed state is illustrated in a. The region within the rectangle 
is considered as viscous and therefore obeys equation (209). The region 
outside the rectangle is perfectly elastic and therefore obeys equation 
(210). Let us now consider that an external load is suddenly applied. With 
no essential loss of generality, we may regard the instantaneous (i.e., the 
unrelaxed) modulus of the viscous region as identical with that of the sur¬ 
rounding elastic medium. Just after the load is applied, the deformation is 
then everywhere uniform, as illustrated in b. The shear stress within the 
viscous region at once starts to relax. If the elastic matrix remained rigid 
during this relaxation, the shear stress in the viscous region would relax 
with the time of relaxation, r, in equation (209). The surrounding elastic 
matrix will not, however, remain rigid during this relaxation. The equa- 

66a. T. S. K6, “Stress Relaxation by Interstitial Atomic Diffusion in Tantalum,” Phys. 

Rev. (1948). 

67. C. Maxwell, Scientific Papers (Cambridge, 1890), II, 623; J. H. Poynting and J. J. 
Thomson, A Text Book of Physics—Properties of Matter (London: Griffin, 1907), p. 51 
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tions of equilibrium require that the shear stress in the surrounding elastic 
matrix suffer at least a partial relaxation, accompanied by shear strains, as 
illustrated in c. This relaxation in the elastic matrix results in an over-all 
slow deformation, or creep. Let us now consider that the load is removed. 
The instantaneous strain which accompanies this removal is uniform 
throughout. Far away from the viscous region the elastic matrix will be 
under essentially zero stress. In the immediate vicinity of the viscous re¬ 
gion the elastic matrix will have, however, residual stresses due to the 
strain which has remained in the viscous region, as illustrated in d. These 
residual stresses, in turn, cause shear stresses in the viscous region. The 
gradual relaxation of these latter shear stresses ultimately brings about 
complete relaxation of all stresses, and the specimen finds itself in its 
initial state as illustrated by a. It is this relaxation of all residual stresses 
which gives rise to the elastic aftereffect. 

One remarkable feature of the two-component system is the large mag¬ 
nitude of the anelastic effects which may be produced by an almost negli¬ 
gible amount of viscous material. Thus suppose that the viscous phase is 
distributed in the form of a thin disk. As already explained, shear stress is 
partially relaxed in the elastic matrix for some distance on either side of 
the disk. The distance over which this relaxation extends is essentially in¬ 
dependent of the thickness of the viscous disk and will be comparable with 
the radius of the disk. If the disk is several atomic diameters wide and has 
a radius of several hundred atomic diameters, shear stress will be relieved 
m a volume of the order of magnitude of a hundred times as large as the 
viscous region itself. The concept of a two-component system therefore 
requires only a negligible amount of the viscous phase. 

A second remarkable feature of the two-component system is the build¬ 
ing-up of high stress concentrations within the elastic matrix as the stress 
rela.xes within the viscous region. The formation of such stress concentra¬ 
tions may readily be understood by reference to Figure 43. While the 
shear stress relaxes in the elastic matrix on either side of the viscous re¬ 
gion, high stress concentrations are built up near the ends of the viscous 
region. The necessity for such stress concentrations may be seen also from 
general energy considerations. If L is the applied load, the total strain 
energy stored in the specimen before and after stress relaxation is propor¬ 
tional to L-, Mu and L- ‘Mr, respectively, where Mu and Mr are the unre¬ 
laxed and rela.xed moduli. Since Mr is less than Mu, the strain energy 
stored in the specimen is greater after rela.xation than before. Since, how¬ 
ever, some regions have suffered a relaxation of stress, other regions must 
be under a correspondingly greater stress. An estimate of the order of 
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magnitude of the stress concentration may be obtained from the computa¬ 
tions of Inglis®^ upon the shear stress computations due to cracks. Thus, 
suppose that the viscous region is in the shape of a thin disk and that the 
applied forces result in shear stresses in the plane of the disk, with no nor¬ 
mal tensile stresses. The strains in the surrounding elastic region in the 
completely rela.xed condition will then be precisely the same as if the vis¬ 
cous region had been a cracklike cavity, for which Inglis made his com¬ 
putations. His results may be e.vpressed most conveniently in terms of the 

radius of the disk, a. and in terms of the radius of curvature, r, of the edge 
of the crack. His results are 


Maximum tensile stress at edge of crack 



X (over-all shear stress). 


Thus, if the radius of curvature at the edge of the disk is only one-hun¬ 
dredth the radius of the disk, stress relaxation within the viscous region 

develops a tensile stress within the elastic matrix ten times as large as the 
over-all shear stress. 

The stress concentrations in the elastic matrix developed by stress re¬ 
laxation in the viscous regions may play an important role in the fracture 
of metals. It has long been recognized that the observed over-all fracture 
strength of metals is at least an order of magnitude less than the value 
computed for a perfect lattice. The theoretical value of the tensile stress at 
fracture, T, may be expressed as 

T= e^, 


where Co is the elastic strain at fracture and E is Young’s modulus. The 
theoretical value of eo for crystals is larger than 0.1; thus for NaCl'" it is 
0.15. On the other hand, the elastic strain at fracture is never so large as 
0.01. One is tempted to attribute the discrepancy between the obse*rved 
and the theoretical fracture stress to the presence of small microcracks, 
which, through their stress concentrations, locally raise the tensile stress 
to the theoretical value for fracture. In fact, many fracture phenomena 
may be understood in terms of the presence of such microcracks.'’ As 


69. C. E. Inglis, “Stresses in a Plate Due to the Presence of Cracks and Sharp Corners ” 
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first pointed out by Griffith,'- a microcrack requires a critical over-all 
tensile stress before it will spread, the magnitude of this critical stress in¬ 
creasing as the radius of the crack decreases. This relation may be under¬ 
stood from two viewpoints, from surface energy considerations and from 
stress concentration considerations. According to the first viewpoint, we 
inquire as to whether an increase in the radius of the crack will result in a 
lowering or a raising of the energ>' of the system. Two types of energy must 
be considered, elastic strain energ}- and surface energ>^ If the tensile stress 
is normal to the plane of the crack and the radius of the crack is a, stress 
will be relieved in a volume of the order of magnitude of (Ltt 3)a'L If the 
over-all strain is maintained constant, the strain energy associated with 
the crack is therefore 

Strain energy = 

where C is a numerical constant of the order of magnitude of unity and T 
is the over-all tensile stress. The surface energy associated with the crack 
is 

Surface energy = 2wa^S , 

where S is the energ}' per unit surface, i.e., the surface tension. The crack 
will spread under the influence of the tensile stress only if the decrease of 
strain energy resulting therefrom is greater than the increase in surface 
energy. The critical tensile stress is obtained by equating to zero the deriv¬ 
ative of the total energ}’ with respect to the radius a. One thereby obtains 

7 - = Vg. (211) 


According to the second viewpoint, a crack will propagate when its stress 
concentration locally raises the over-all tensile stress to the theoretical 
fracture stress. Since the stress concentration factor'^ is 2Va r, where r is 
the radius of curvature of the edge of the crack, this second viewpoint 
gives 



( 212 ) 


According to both vie^vpoints, T varies inversely as the square root of a. 
In order that equations (211) and (212) may be compatible, it is nece3sar>' 
that 
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The left-hand side is the work done in opening up a crack of unit area. 
When we take r as one interatomic distance, it is clear that this equation 
will be satisfied by a numerical constant, C, of the order of magnitude of 
unity. Griffith's concept of a critical tensile stress for the projiagation of 
a crack poses the question as to the origin of the original microcracks. 
Their origin very probably lies in the stress concentrations developed by 
stress relaxation within amorphous regions. The role of stress relaxation 
in opening up microcracks will be further discussed in subsequent sections, 
where the different types of viscous regions are treated. The development 
of microcracks by localized stress relaxation has been considered by Mur- 
gatroyd"’ as the origin of stress fatigue in glass, i.e., fracture by a sustained 
load. 

A third remarkable feature of the two-component system is the wide 
variety of types of relaxation spectra which it may have. If all the local¬ 
ized viscous regions had essentially the same size and shape, the rela.xation 
spectrum could be described essentially by a single time of relaxation, as is 
frequently the case for the relaxation of thermodynamic potentials. A re¬ 
view of past e.xperimental work upon anelasticity presumably arising from 
viscous regions is given by Fromm.In no cases did the observations indi¬ 
cate a relaxation spectrum even approaching that for a single relaxation 
time. In fact, all observations indicate a relaxation spectrum whose intens¬ 
ity continually rises or remains stationary with increasing times of relaxa¬ 
tion. No cases have been observed in which the intensity of the rela.xation 
spectrum decreases with increasing times of relaxation. These observa¬ 
tions may be interpreted in several ways. 

As we have already mentioned, a concentrated relaxation spectrum 
would be observed if all the localized viscous regions were of the same 
shape and size. One might interpret the observed anelastic behavior in 
terms of a distribution of size and shape of the viscous regions. Thus sup¬ 
pose that the regions were disk-shaped and that all had the same thickness 
but various radii. Then the total shear strain for complete stress relaxa¬ 
tion would be larger in the disks of greater radii, and hence the larger 
viscous regions would have longer relaxation times. 

On the other hand, the viscous phase might be regarded as forming a 
continuous network which incloses the elastic matrix. If this network 
were in the form of regular polyhedra, then the edges and corners of the 
polyhedra would block viscous deformation, the viscous deformation be- 

J* Murgatroyd, Jour. Soc. Glass Tcchnol., XXVIII (1944), 406: Nature, CLIV 
(1944), 51; J* B. Murgatroyd and R. F. Sykes, Nature, CLVI (1945), 717. 
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ing thereby confined to the faces of the polyhedra. If the polyhedra were 
all of the same size, their relaxation times would all be equal, and the re¬ 
laxation spectrum would be concentrated. On the other hand, a wide dis¬ 
tribution in sizes of the polyhedra would give rise to a wide relaxation spec¬ 
trum. Such a viscous network could give rise to a wide relaxation spectrum 
in still another manner. If the network could be represented by regular 
polyhedra, viscous flow might extend over wide areas unimpeded by edges, 
while viscous flow in other areas might be limited by closely spaced edges. 

The form of the distribution function, A(r), and its methods of deter¬ 
mination have recently become a subject of intense concern among those 
interested in the mechanical properties of high polymers. The activity'® in 
this field during 1944-46 is in marked contrast to the almost complete lack 
of work in the corresponding field of metals. 


II. STRESS RELAXATION ALONG PREVIOUSLY FOR.\IED SLIP BANDS 
(I) ROLE OF SLIP BANDS AND DISLOCATIONS IN PLASTIC DEFORMATION 

That plastic deformation in metals occurs within isolated regions, slip 
bands, rather than homogeneously throughout the material, was first 
clearly recognized by Ewing and Rosenhain." Very careful examination 
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of the material between the slip bands failed to reveal any deformation 
whatsoever. Deformation through isolated slip bands rather than homo¬ 
geneous deformation appears, in fact, to be characteristic of all crystalline 
materials. As will be discussed later in some detail, the evidence is very 
strong that the material in a previously formed slip band is, at least tem¬ 
porarily, viscous in the sense of obeying equation (209) rather than equa¬ 
tion (210). In polycrystalline materials each slip band is confined to a 
single grain. It would therefore appear as if each slip band behaved as an 
isolated viscous region and thereby gave rise to anelastic effects. The 
properties of slip bands are hence ver\' pertinent to a study of the anelas- 
ticity of metals. 

It is, in fact, the viscous nature of slip bands which renders the yield 
stress of pure metals so low, two to three orders of magnitude below the 
theoretical yield stress of where cq is of the magnitude of 0.1. As pre¬ 
viously mentioned, an isolated viscous region will result in shear-stress 
concentration. Just in front of the spearhead of an advancing slip band 
the local shear stress is given by'® 

Local shear stress ^ Vy over-all shear stress, 

where L is the length of the slip band and r the radius of curvature of its 
spearhead. According to this equation, the larger the slip band, the lower 
need be the over-all stress for continued propagation. The relation be¬ 
tween critical stress and slip-band size is corroborated by two observa¬ 
tions. First, the formation of a slip band appears to be a cataclysmic proc- 
ess.'^^ Once the stress is sufficiently high to cause propagation, one would 
expect the slip band to grow at an increasing rate. Second, observations 
upon the proportional limits of copper and copper-base alloys®*^ are in 
agreement®* with the theoretical prediction based upon equation (211) 
that the stress at which yielding is first observed is inversely proportional 
to the grain size, i.e., to the maximum value of L. 

The concept of slip bands propagating under their own stress concen¬ 
tration encounters the same difficulty as does the concept of cracks propa¬ 
gating under their own stress concentration. This concept is powerless to 
describe how the slip bands originate; for, until they attain a certain criti¬ 
cal size, which is smaller, the higher the stress, their stress concentration 

78. Inglis, op. cit. 

79. I. W. Obreimow and L. W. Schubnikoff, Zeitschr.f. Phys., XLI (1927), 907. 
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is not sufficiently high to cause self-propagation. Some other mechanism of 
deformation must e.xist which leads eventually to the initiation of slip 
bands. The most striking evidence for pre-slip-band deformation is to be 
found in the e.xperiments upon a-brass*- and a-iron^^ single crystals. In 
these crystals the application of a load may be accompanied by no im¬ 
mediate observable deformation. Only after an “incubation” period of the 
order of minutes or hours does visible deformation start, and it then pro¬ 
ceeds at an increasing rate. Less pronounced incubation periods have been 
found in commercially pure zinc single crystals.*^ It is possible that during 
the incubation period small slip bands are being formed through the move¬ 
ment of lattice defects, such as dislocations. The concept of dislocations 
was originated by Orowan**^ and Polanyi®® and greatly extended by Tay¬ 
lor,®" in order to explain why metals yield at a stress several orders of mag¬ 
nitude less than the theoretical yield stress. .Although dislocations have 
been the subject of considerable theoretical investigation, a review of 
which has been given by Seitz and Read,®® it is not clear just what is their 
role in plastic deformation, nor is it clear how they can initiate slip bands. 
Polycr>'stalline lead, freshly quenched from an anneal at I20°-I30 C., 
shows an initially increasing creep rate,®* a phenomenon which appears to 
be related to the incubation periods in a-brass and a-iron ciy’stals. 

From the above discussion we conclude that deformation proceeds pri¬ 
marily through the propagation of slip but that another mechanism of de¬ 
formation also exists which, although contributing a negligible amount to 
the total deformation, is nevertheless essential to the initiation of slip 
bands. Before we leave this subject, it will be pertinent to review other ex¬ 
perimental evidence which indicates deformation by some mechanism 
other than slip bands. In his experiments upon tin, Chalmers*® has found 
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that at ver)' low stresses the initial rate of creep is proportional to the 
applied stress but that the total extension produced by this creep is limited 
and is independent of the applied stress. It looks as though his specimens 
contained dislocations which moved at a rate proportional to the stress but 
that they could move only a definite distance, after which they were pre¬ 
sumably stopped by some sort of barrier. 

Plastic deformation usually increases anelastic effects. Examples have 
been reported, however, in which the prior application of a stress decreases 
the anelastic effects. In his experiments upon single and multiple crystals. 
Found*-*’ discovered that the internal friction at low-stress amplitudes 
could frequently be reduced by the prior application of a stress below the 
elastic limit and that the decrease was proportional to the applied stress. 
These e.xperiments might be interpreted in a manner similar to that used 
in Chalmer’s experiments. Found’s specimens may have initially con¬ 
tained dislocations which contributed to the internal friction. A steady 
stress presumably moved these dislocations as far as a barrier, to which 
they appear to become attached and hence no longer contribute to the in¬ 
ternal friction. 

Independent evidence that deformation by slip bands is preceded by 
another type of deformation may be found in experiments upon rock salt. 
Here the yield stress - the stress necessary for the initiation of slip bands 
—is very- dependent upon purity, the addition of only 0.001 mole per cent 
of PbCb increasing the yield by a factor of On the other hand, the initi¬ 
ation of deformation in ionic crystals may also be detected^* by the intro¬ 
duction of “color centers.” The yield stress defined in terms of these color 
centers is, however, found to be quite insensitive to the degree of purity.*^ 
It looks as though the color centers introduced by stressing are the fore¬ 
runners of slip bands and that the hardening influence of impurities is re¬ 
lated to the inhibition which they exert in the transition of these color 
centers into slip bands. Although color centers have been extensively stud¬ 
ied,their relation to plastic deformation is at present unknown, aside 
from the fact that they are introduced by application of stresses. 
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b) RELATION OF PREMOUSLY FORMED SLIP BANDS TO 
MECHANICAL PROPERTIES OF METALS 

i. Hysteresis .—During the last centurj' it was observed that plastic de¬ 
formation introduces hysteresis loops in stress-strain curves.®® Thus the 
stress-strain curve in annealed mild steel is straight and reversible, pro¬ 
vided that the yield stress is not exceeded. Once the yield stress is exceed¬ 
ed. with accompanying plastic deformation, the stress-strain curve is no 
longer reversible, even at stress levels below the original yield stress, as 
illustrated in Figure 44. It was soon discovered that this hysteresis dis¬ 
appears with aging,®' slowly at room temperature, within a few minutes 
at 100® C. Rosenhain®® early recognized that this hysteresis and recovery 
could be interpreted in terms of the then recently discovered slip bands. 
He pointed out that the slip bands behave as if the material therein had 
an initially amorphous structure with corresponding viscous mechanical 
properties and as if this material gradually acquired the ciy'stalline struc¬ 
ture of the surrounding undeformed matrix. 

Unfortunately, considerable confusion has arisen through the associa¬ 
tion of Rosenhain’s concept of viscous-like slip bands with Beilby's®® con¬ 
cept of amorphous layers. Beilby also considered the slip band material to 
be amorphous, but he assumed this amorphous structure to strengthen 
rather than to weaken the slip bands, and he did not believe that the mate¬ 
rial gradually acquired the crystalline structure of the surrounding matrix. 
This confusion has led to a gradual discrediting of Rosenhain’s concept of 
viscous-like slip bands. Alternative theories have therefore been advanced 
to interpret the hysteresis introduced by plastic reformation. One current 
view is that the hysteresis is due to residual stresses introduced by plastic 
deformation. While such residual stresses could give rise to hysteresis at 
high stress levels, such hysteresis would vanish as the stress level was low¬ 
ered. The hysteresis, or internal friction, approaches, how'ever, a constant 
value as the stress level is lowered. Another argument contrary to the view 
point of residual stresses is that there is no evidence that the residual mi¬ 
croscopic stresses are eliminated by low-temperature anneals at 100 C. 
Another current view is that the hysteresis introduced by plastic deforma- 
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tion arises through the partial mobility of dislocations introduced by de¬ 
formation. As previously mentioned, it looks as though all deformation is 
confined to localized slip bands. The dislocations, if present, must then 
be within the slip bands themselves. However, a high concentration of 
dislocations would impart to the slip bands a viscous-like behavior. The 
viewpoint of a high concentration of dislocations within slip bands is 
therefore identical with Rosenhain’s viewpoint of viscous-like behavior. 

The development of electronic methods of measuring internal friction at 
low stress amplitudes has led to a considerable number of experiments^®” 
on the influence of plastic deformation. .All these experiments confirm the 



Fic. 4-1.—Exani[>le of hysteresis introduced by [)lastic deformation 


early observations that plastic deformation introduces hysteresis and that 
this hysteresis is removed by low temperature anneals of from 100° to 
200° C. In his work upon the effect of stress amplitude, Read^®' has found 
that, although the internal friction introduced by elastic deformation is 
independent of stress amplitude at very low stress levels, it increases rap¬ 
idly with stress at comparatively low stress levels. Vibration at the higher 
stress levels does not affect the internal friction measured at the lower 
stress levels, i.e., internal friction is a single-valued function of stress 


100. F. Forster and W. Koster, “Elaslizltats-Modul und Dampfung in .\bhanigkeit von 

Werkstoffzusland,” Zeilschr. f. MtUallk., XXIX (1937), 116; W. Koster and K. Rosenthal, 
“Die Anderung von Elastizitats-Modul und Dampfung bei der Verformung und Rekristallisa- 
tion von Messing,” Zeilschr.f. Metallk., XXX (1938), 343; \V. Koster. “Elastizitatsmodul und 
Dampfung von Aluminum und Aluminiumlegierung,”/«/rc/i/-./. XXXII (1940), 282; 

T. A. Read, “The Internal Friction of Single Metal Crystals,” Phys. Rev., LVIII (1940), 371, 
and “Internal Friction of Single Crystals of Copper and of Zinc,” Tram. A.I.M.E., CXLIII 
(1941), 30; W. A. Lawson, “The Effect of Stresses on Internal Friction in Polycrystalline Cop¬ 
per,”/*Ay5. Rev., LX (1941),330; C. Zener, H. Clarke, C. S. Smith,“Effect of Cold Work and 
Annealing upon Internal Friction of Alpha Brass,” Trans. A.I.M.E., CXLVII (1942), 90. 

101. Op. eit. 



138 


ELASTICITY .AND .ANELASTICITY OF MET.ALS 


amplitude. Similar observations were also made by Norton’®- at much 
higher stress levels. 

Considerable uncertainty e.xists at present regarding the relaxation 
spectrum associated with plastic deformation. On the one hand, many 
measurements’®^ made at and in the vicinity of room temperature indi¬ 
cate that the internal friction introduced by plastic deformation is inde¬ 
pendent of the frequency of measurement and therefore that the relaxation 
spectrum is constant. On the other hand, one single piece of research’®'* in¬ 
dicates that the internal friction associated with prior plastic deformation 
rises rapidly as the period of vibration is increased. By observations upon 
the elastic aftereffect of previously deformed iron at elevated temperatures, 
West'®* concluded that the internal friction increased with an increase in 
temperature and therefore, owing to the close relation of the effect of 
changes in temperature to that of changes in frequency of vibration, that 
it also increased with period of vibration. These conclusions are corrobo¬ 
rated by current experiments by T. S. Ke.’®^ It therefore appears that most 
information may be obtained regarding the mechanical properties of slip 
bands by experiments conducted at long relaxation times. Creep or stress 
relaxation experiments are particularly suited for obtaining this infor¬ 
mation. 

ii. Creep .—As mentioned above, it appears that most information re¬ 
garding the mechanical properties of the material within slip bands may 
be gathered from creep experiments at low stress levels and or from the 
associated elastic aftereffect e.xpcriments. Numerous experiments have 
demonstrated that the elastic aftereffect is increased by cold working. 
These hav? been reviewed by Fromm.’®® The first attempt to gain an in¬ 
sight into the mechanical properties of the slip band material through such 
studies has been made by West,*®* who found that his observations upon 
plastically deformed iron at a series of temperatures could be correlated 
by ascribing to the viscosity of the slip-band material a heat of activation 
of 21,000 cal/mole. West's interpretation of his observations m terms of 
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slip-band viscosity has been questioned by Ke,*°" who has suggested that 
the effects observed by West were probably associated with the high- 
temperature (150°-200° C.) relaxation in carbon- and nitrogen-bearing 
iron induced by plastic deformation,the origin of which is not at 
present understood. 

It is also commonly supposed that deformation necessarily raises the 
resistance to further deformation. This is the case only at lower tempera¬ 
tures. Extensive experiments, a review of which has been given by Bur¬ 
gers,^®® have shown that, at elevated temperatures, previously deformed 
specimens creep at a greater rate than do virgin specimens. Such acceler¬ 
ated creep presumably arises from the viscous behavior of the previously 
formed slip bands. 

Once an insight has been gained, through creep experiments, into the 
viscous behavior of slip bands, it is anticipated, conversely, that this in¬ 
sight will give a better understanding of creep phenomena. In corrobora¬ 
tion of this viewpoint, evidence will first be presented to show that the 
viscous behavior of slip bands is an important factor in the variation with 
temperature and with rate of deformation of the resistance to deformation. 
The dominating influence upon creep will then be discussed. 

It is commonly supposed that the yield strength—the stress necessary 
to produce the first plastic deformation—decreases with a rise in tempera¬ 
ture. When it is realized that the yield stress is commonly taken as the 
stress at which the permanent set is of the order of magnitude of 0.001 and 
that a considerable amount of strain hardening can occur during the first 
0.001 strain, the commonly accepted viewpoint is open to considerable 
doubt. It might, on the other hand, be that the true yield stress, if such 
exists, is relatively insensitive to temperature and that the apparent rise of 
yield stress with a lowering of temperature is due merely to the increased 
rate of strain hardening at lower temperatures. Several experiments indi¬ 
cate that this viewpoint is indeed correct. In their careful observations 
upon the stress-strain curves of very pure single crystals of aluminum and 
silver. Miller and Milligan^’® have found that, as the temperature is raised 
from room temperature up to the recrystallization temperature, the stress 
actually increases, provided that the strains are less than 0.001. Only for 
larger strains does the stress decrease with an increase in temperature. At 
least for these crystals, the usual observation of decrease of yield stress 

107. T. S. K6, personal communication. 

108. J. L. Snoek, Pkysica, VIII (1941), 711. 

109. W. G. Burgers, Uandbuch dtr Melallphysik {Leipzig, 1941), pp. 504-7. 
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with increasing temperature is due merely to a decreased rate of strain 
hardening. In their work upon the initiation of yield in a-brass crystals, 
Burghoff and Mathewson'" found that the critical resolved shear stress at 
which creep finally started, perhaps after a long incubation period, was, 
within their limits of error, independent of temperature from room tem¬ 
perature to 500° F. From the above experiments we are led to suspect that 
the commonly observed decrease in strength with increasing temperature 
is due primarily to viscous flow in previously formed slip bands, such flow 
thereby decreasing the measured strain hardening. The higher the tem¬ 
perature, the more rapid this viscous flow, and therefore the more is the 
strain hardening reduced. Again, the lower the rate of strain hardening the 
more dominant is this viscous flow, and therefore the more is strain hard¬ 
ening reduced. Only at the lowest temperatures is viscous flow within pre¬ 
viously formed slip bands entirely suppressed and hence is strain harden¬ 
ing unaffected by such flow. 

We shall now examine in detail the interrelation between strain harden¬ 
ing and creep, i.e., tensile deformation at constant load. A typical creep 
curve is reproduced in Figure 45. The rapidly decreasing creep rate in the 
first portion, in the primary creep range, is customarily attributed to 
strain hardening. We shall see below that equally important in the primary 


creep range is the variation of stress with strain rate, which, according to 
our present view, depends upon viscous flow in previously formed slip 
bands. The stationary creep rate of the secondary creep range is custom¬ 
arily believed to arise from an interaction of strain hardening and strain 
recovery. A criticism of this viewpoint is presented below, and an alterna¬ 
tive viewpoint"' is discussed. According to the latter viewpoint, the 
stationary secondary creep rate is dependent upon precisely the same 
metallurgical factors as those which govern the primary creep range. Final¬ 


ly, the rapidly increasing creep rate of the tertiary range is due either to 
the initiation of necking or to the formation of a multitude of small cracks 
within the specimen. The origin of these cracks is discussed under the sec¬ 


tion on stress relaxation across grain boundaries. 


We shall now analyze the factors upon which primary creep depends. 
As discussed in detail by Hollomon,"^ the initial strain hardening of 
metals at constant strain rate is represented very well by the expression 


e", e - Constant. 


(213) 
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The strain hardening exponent, m, is equal to the strain at which necking 
begins"^ and has been found to range from 0.06 for low-carbon steel to 0.7 
for large grain size brass. On the other hand, at least in steel, stress varies 
with strain rate at constant strain, as"^ 

S'^ e’', Constant. (214) 

It is anticipated that the exponent n will be dependent upon temperature 
and somewhat upon the stress level. In steel at room temperature and 
below, n is of the order of magnitude of 0.01. Upon combining equations 
(213) and (214), we obtain 

S'^ e”' e" . (215) 



Fig. 45.—Typical creep curve (elongation versus time at constant load) 


If we regard m and n as constants, the solution of equation (215) is 

c = CrA”+"*) , (216) 

where C is a constant. The creep rate, which is therefore given by 

e= /,/-">/("+•") , (217) 

where b is another constant, approaches zero asymptotically. This conclu¬ 
sion is in agreement with the careful work of Andrade,”® who showed that 
at constant stress the creep rate approached zero as The rate at 
which the strain rate approaches zero is dependent both upon the strain 
hardening exponent m and upon the dependence of stress upon strain 
rate, a measure of which is given by the exponent n. 

The customary creep experiments, to which Figure 45 refers, are not 
performed at constant stress but rather at constant load. As the specimen 

U4. Ibid. 
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extends, the cross-section becomes smaller, and therefore the load in¬ 
creases. Although the increase of stress is negligible during the early part 
of the primary creep range, it is not necessarily negligible throughout the 
entire stress range. In order to find its influence, we must replace equation 
(215) by an equation in which load L appears rather than stress. If A refers 
to the cross-section area, then 

L = .15 . 

If we neglect volume changes during creep, 

.1 , 

where € is the base of the natural logarithm. The desired equation is 
therefore 

. ( 218 ) 

The creep rate is seen to be a minimum at the strain m and to increase 
thereafter. An example of the integral of equation (218) is given in Figure 
46 for the particular case of m n - 10. The creep rate is seen to have an 
essentially constant minimum value for a very tong time. This result is 
a general consequence of equation (218) or of any other equation of the 
same general type; for, according to this or any similar equation, the creep 
rate is a minimum at some value of strain, say e,n. The creep rate will 
hence be essentially constant for strains some distance on either side of 
Cm- Since this essentially constant strain rate is also a minimum strain 
rate, a relatively long time will be spent in the region of essentially con¬ 
stant strain rate. 

According to the above viewpoint, which was originally formulated by 
Ludwik,"' the phenomenon of a secondary creep range is merely a conse¬ 
quence of the variation of stress upon strain and upon strain rate, and it is 
not necessary to invoke the additional concept of softening by annealing. 
Although it is known that strain hardening is removed at elevated tem¬ 
peratures, it is also known that the smaller the prior strain, the higher one 
must go in temperature to obtain strain recovery. No experiments have 
been performed to date to determine whether in the range of strain and 
temperature which are of practical importance to creep, appreciable 
strain recovery takes place. 

iii. Stress relaxation .—As pointed out in chapter v, section B, precisely 
the same information may be obtained by stress relaxation experiments at 
low stress levels as by creep experiments at low stress levels. The viscous 
properties of slip band material may therefore be investigated by stress 

117. P. Ludwik, FJemenIe der technologischen Meclhinik (Berlin: Springer, 1909). 
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relaxation experiments. Conversely, any information so gained may be ap¬ 
plied to reach an understanding of stress relaxation under practical condi¬ 
tions of operation. 

Stress relaxation at constant strain is intrinsically a simpler physical 
phenomenon than is creep under ordinary conditions. The major fraction 
of an applied stress may be relieved through viscous flow in previously 
formed slip bands or along grain boundaries, and no generation of new 
slip bands is required. However, under ordinary creep conditions an in¬ 
crease in strain is presumably attended by the generation of new slip 



Tim# 


Fig. 46.—Theoretical creep curve according to equation (180) (drawn for case m,n = 10) 


The results of stress relaxation experiments are usually plotted in a 
semilogarithmic fashion, i.e., stress versus logarithm of time. In their work 
upon lead, Trouton and Rankine,”* presumably the first to study stress 
relaxation extensively, found that in such a plot the major part of their 
data lay upon a straight line. Thus in this region 

S = a-b\ogl, (219) 

where a and b are constants. Later investigators have found essentially 
the same result. 

Equation (219) does not behave properly for small times, since S{1) 
must approach a constant value 5(0) as t approaches zero. In order to 
remedy this defect, another constant is necessary. The new equation may 
be written as 

5(/) =5(0)1(220) 
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This new equation, as well as the original equation (219), behaves im¬ 
properly as t becomes very large. In order to avoid this difficulty, Trouton 
and Rankine have suggested still another modification, with, of course, the 
introduction of still another constant. This new equation is 

5a)=5(0)l-C/«j-|^[ (221) 


where 
and where 


r2»Ti 



( 222 ) 


The equality sign in equation (222) is to hold in case of complete stress 
relaxation at infinite time. 

No satisfactory theoretical justification has been presented for an equa¬ 
tion of the form of (219) for intermediate time ranges. A theoretical deri¬ 
vation has been given by Tobolsky and Eyring.'*^ An assumption neces¬ 
sary to their argument is, however, that the creep rate is not linear in 
stress, but rather 

sinh A 5 ; 


and their conclusion is valid only when A5 is so large that the sinh func¬ 
tion cannot be replaced by its argument. Equation (219) has been found 
valid, however, in the low stress range where creep rate is proportional to 

stress. 

iv. Moduli of elasticity 1907 Beilby*^® observed thatplastic 
deformation decreased the moduli of elasticity and that this decrease was 
removed by a low-temperature anneal. Some confusion arose in later years 
from observations which indicated that, at least in certain cases, the mod¬ 
uli were increased. A clarification was finaUy given by Kawai''* in 1930. 
Kawai found that a slight deformation always decreases the moduli and 
that in certain metals continued deformation leads finally to an increase 
in the moduli. A review of the literature on the subject has been given by 
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Elam.'" Since this review, a thorough investigation by Koster'-^ shows 
that the restoration of the original moduli on annealing closely parallels 
the removal of the internal friction introduced by plastic deformation. 

A ready interpretation of the effect of plastic deformation in lowering 
the elastic moduli is to be found in the viscosity of the slip bands produced 
by this deformation. In fact, the decrease in elastic moduli is just one man¬ 
ifestation of this viscosity, of which hysteresis, creep, and stress relaxation 
are other manifestations. The simplest relation of the decrease of elastic 
modulus is with the stress relaxation function. Thus let f{i) be the force 
necessary to cause a constant unit deformation to be suddenly applied at 
i — 0, and let M{P) be the corresponding dynamic modulus measured at 
a period of vibration of P. Then from equation (67) we see that 

M{P) =/(0,^P/8- 

This interpretation of the origin of the decrease in elastic moduli with 
plastic deformation leads to the prediction that no decrease in modulus 
would be observed if the measurements were carried out at very low tem¬ 
peratures. Such measurements have not been made. 

The final rise in elastic moduli in certain metals due to continued plastic 
deformation presumably arises from reorientation effects which are as¬ 
sociated with extensive deformation. 

V. Bauschinger effect .—Annealed metals usually have the same yield 
stress in compression as in tension. Bauschinger'*' discovered in 1881 that 
this symmetry was destroyed by plastic deformation. If a metal is 
stretched in tension, its strength is raised with respect to a second applica¬ 
tion of a tensile stress but is lowered with respect to the application of a 
compressive stress. In fact, the metal yields plastically almost as soon as 
even a small compressive stress is applied. 

The customary interpretation of the Bauschinger effect is in terms of 
the anisotropy of individual crystals with respect to their yield strength 
and their at least partially random orientation. When a tensile stress is 
slowly applied to a specimen, those grains yield first which are most favor¬ 
ably oriented for plastic deformation. As the stress is further increased, 
the stress increases less rapidly in the grains undergoing plastic deforma¬ 
tion than in those grains which have remained elastic. Now, if the applied 
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stress is removed, the grains which had undergone plastic deformation and 
hence had the smallest peak tensile stress will now be in compression, 
while the grains which had remained elastic will be in tension. 

This customary interpretation of the Bauschinger effect can at most be 
a very rough approximation to the actual conditions in a plastically de¬ 
formed specimen, since the deformation is not uniform throughout an in¬ 
dividual grain. It does not, therefore, have much meaning to speak of the 
“residual stress” within an individual grain. It would appear more appro¬ 
priate to speak in terms of individual slip bands. Just after a slip band has 
formed, the shear stress within it and in the elastic matrix on either side is 
to a great extent entirely relieved. If. now. the applied stress is removed, 
these regions will have shear stresses of the same magnitude but of a sign 
opposite to the original over-all stress. These regions will therefore yield 
under a very slight shear stress applied in the opposite direction. An inter¬ 
pretation of fatigue in metals may be based upon just these residual 
stresses.*-^ 

The observations of Sachs and Shoji'-* upon single crystals corroborate 
the above view'point. They found that single crystals show the Bauschinger 
effect to just as marked a degree as do polycrystalline specimens. It 
therefore appears that it is the residual stresses associated with individual 
slip bands rather than with individual grains which are responsible for the 
Bauschinger effect. 

The foregoing interpretation leads to certain interesting predictions. 
The residual stresses across individual slip bands should be to a large ex¬ 
tent removed by low-temperature annealing. The Bauschinger effect 
should likewise be so removed. On the other hand, it should be possible to 
reintroduce the Bauschinger effect without further plastic deformation. 
Application of a high stress at a slightly elevated temperature should re¬ 
introduce the residual stresses. 

vi. Fracture —delayed fracture phenomena are most readily 
interpreted in terms of the gradual relaxation of shear stress across pre¬ 
viously formed slip bands, with a resulting high accumulation of stress 
concentration at the edges of the slip bands. Most striking examples have 
been observed*” by the author in armor-piercing projectiles. When an un¬ 
capped armor-piercing projectile strikes an armor plate under sufficiently 
severe conditions at normal incidence, the projectile penetrates in an in- 
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tact condition but with the front portion having sulTcred a slight comi^res- 
sion. After the projectile is recovered, its tip frc(iucnt!>' dies olT. An exam¬ 
ple is given in Figure 47. The delay between recovery and fracture varies 
from several minutes to several days. The force system acting u]>on the 
projectile is such that, after penetration, the central part of the front por¬ 
tion is under tension, the outer parts under compression. After recover)- 
the residual tensile stress is presumably partially relaxed by shear stress 
relaxation across previously formed slip bands, this relaxation resulting 
in such high stress concentrations that small microcracks are eventualU- 
formed, which, in turn, are further of^ened by the residual tensile stress. 
The flying-otT of the tip finally results from the acceleration by the com¬ 



pressive stresses along the outer portion after the spreading of the crack 
throughout the central regions. 

The rather frequent cases of delayed cracking in steel following cjuench- 
ing‘-* may also be the result of the gradual building-up of high stress con¬ 
centrations through stress relaxations across previously formed slip bands. 
A spectacular example has been described by Howe. “In the early days of 
making armor-piercing shells, spontaneous and violent aging rupture was 
so common that shells, after hardening, used to be stored for a considerable 
time in a room to which no one was admitted. 


III. STRESS RELAXATIO>J ACROSS GRAIN BOU.VDARIES 
a) VISCOSITY OF GRAIN BOUNDARIF.S 

Since about 1912 the existence of amorphous grain boundaries has been 
the subject of almost continual controversy. On the one hand, the evi- 
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dence, as reviewed first by Rosenhain’^® and later by Jeffries and Arch¬ 
er,'^* is very convincing that grain boundaries behave as if they were 
amorphous. On the other hand, metallurgists have been reluctant to ac¬ 
cept the concept of an amorphous phase. No controversy need arise, once 
it is realized that it is not necessary for any portion of the metal to be 
amorphous in order that the grain boundaries may behave in a viscous 
manner. It is necessary to assume only that the resistance to slipping of 
one grain over an adjacent grain obeys the laws commonly associated with 
amorphous materials rather than the laws associated with crystalline ma¬ 
terials. Since the surface atoms of one grain cannot fit into the lattice posi¬ 
tions of an adjacent grain, the binding across the interface of two grains 
may reasonably be expected to have the characteristics associated with 
amorphous materials. 

The primary difference between the laws governing the mechanical be¬ 
havior of amorphous and crystalline materials lies in the rapid increase in 
resistance to deformation of amorphous materials with respect to an in¬ 
crease in rate of deformation and to a decrease in temperature, compared 
with the relative insensitivity, with respect to these same variables, of the 
resistance to deformation of crystalline materials. Thus in amorphous ma- 

4 

terials of simple structure 

Shear stress = ?? X shear strain rate , 


where the coefficient of viscosity, t;, varies with temperature as 

e»/nT ^ 

where R is the gas constant and Hy the heat of activation, has a magnitude 
of several tens of thousands cal/mole. On the other hand, in the typical 
case of steel in the vicinity of room temperature, the shear stress varies as 
only about the one-hundredth power of the strain rate,*^- and the heat 
of activation, which describes the dependence upon temperature, is of 
the order of magnitude of only 100 cal/mole. It is therefore to be expected 
that, as the rate of deformation is lowered or as the temperature is raised 
the resistance to slipping across the grain boundaries wiU be lowered with 
respect to the resistance to deformation within the interior of the grains. 
At sufficiently low rates of deformation and/or at sufficiently high tem¬ 
peratures, effects should be observed which are attributable only to a slip- 
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ping across the grain boundaries. Such effectswere first observed by Rosen- 
hain and Humphrey When specimens of copper were polished and then 
extended at a slightly elevated temperature, those authors found that the 
grain boundaries became delineated by a relative displacement normal to 
the surface of adjacent grains. Such delineation does not occur at room 
temperature. Hanson and Wheeler*^^ made a thorough study of grain 
boundary movement in aluminum. Here considerable relative movement 
is observed at the grain boundaries when the specimens are extended slow¬ 
ly at temperatures of 250° C. and above, but not when the specimens are 
extended rapidly an equivalent amount at room temperature. Andrade 
and Chalmers^^^ found that a slight extension of cadmium at room 
temperature results in a decrease in electrical resistivity, while a like ex¬ 
tension at liquid-air temperature raises the resistivity. It looks as though 
at room temperature the grains rotate in such a manner as to bring the 
axis of low resistivity closer to the tensile axis, without appreciable plastic 
deformation within the grains. Such a rotation could occur only by slip¬ 
ping along the grain boundaries. Moore, Betty, and Dollins’^® have ob¬ 
tained most striking evidence of viscous flow at grain boundaries. When 
specimens of lead are polished, scratched with parallel straight lines, and 
then pulled several per cent at room temperature, Moore el al. found that 
the scratches remained parallel in each grain but assumed different orien¬ 
tations in adjacent grains when the extension was slow and maintained the 
same orientation when the extension was rapid. Schumacher’®^ has pointed 
out that discontinuities in the creep curves of lead are favored by low 
stresses and by large grains, and he has interpreted such discontinuities in 
terms of viscous grain-boundary slip. A review of the literature comparing 
the low creep resistance of fine-grained specimens in contrast to that of 
large-grained specimens has been given by Burgers.’®* The conclusive 
proof that grain boundaries do, in fact, behave in a truly viscous manner 
has recently been given by Ke, whose extensive work is referred to in the 
following sections. 
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b) RELATION OF THE NTSCOSITY OF GRAIN BOUNDARIES TO 
THE MECHANICAL PROPERTIES OF METALS 

i. Hysteresis .—Viscous flow is inevitably associated with the dissipa¬ 
tion of mechanical energy into heat energy. Viscous flow at the grain 
boundaries should therefore contribute to the internal friction of metals 
subject to cyclic vibration. As long as the frequency of vibration is suffi¬ 
ciently high so that the shear stress across the grain boundaries rela.xes 
only slightly during a half-cycle of vibration, the internal friction is ex¬ 
pected to be higher, the higher the temperature and the greater the grain 
boundary area per unit volume, i.e., the smaller the grain size. These ex¬ 
pectations have been confirmed by observations upon polycrystalline 
zinc'-'**-’ and brass^^® specimens at a variety of grain sizes and over a wide 
range of temperature. In all cases the internal friction varied inversely as 
the grain size and increased with rising temperature according to a heat 
of activation type of law. 

If the individual grains are essentially equiaxed and the grain-size dis¬ 
tribution is uniform, it is anticipated that only a fraction of the over-all 
stress can be relaxed by grain-boundary slip. The situation is analogous 
to the case of a jigsaw puzzle, in which the over-all configuration possesses 
rigidity in spite of the fact that no shearing stresses exist between adja¬ 
cent pieces. In this case a continued rise in temperature will not lead to a 
continued rise in the internal friction. Thus the energy dissipated across a 
single boundary may be written as 

Energy dissipated'^ (relative displacement) X (shear stress). 

At low temperatures the first factor is negligible, since here the viscosity 
is so high that no appreciable relative displacement of adjacent grains oc¬ 
curs during a half-cycle of vibration. At high temperatures the second 
factor is negligible, since here the viscosity is so low that the shear stress 
across all grain boundaries is essentiaUy completely relaxed at all times. 
Only in an intermediate temperature range, where there is some relative 
displacement and some shear stress across the grain boundaries, is the dis¬ 
sipation of energy appreciable. An example of such an optimum tempera¬ 
ture range has been demonstrated by Ke.^^^ His results are shown in 

Figure 48. 

The magnitude of the internal friction in the optimum temperature 

139. A. Barnes and C. Zener, "Internal Friction at High Temperatures,” Phys. Rev., LV III 

C. Zener, D. Van Winkle, and H. Nielsen. "High Temperature Internal Friction of 
.Mpha Brass,” .<4./.i/.£., CXLVTI (1942), 98. j • in 

141 T S. K6, "Experimental Evidence of the Viscous Behavior of Gram Boundaries m 
Metals,” Phys. Rev., LXXI (1947), 533. 
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range should be independent of grain size. This conclusion may be ar¬ 
rived at from two vie^^"{)oints. In the optimum temperature range the in¬ 
ternal friction is proportional to the grain-boundary area per unit volume 
times the average relative displacement across a grain boundar)’. Now the 
grain-boundary area is inversely proportional and the relative displace¬ 
ment is directly proportional to the grain size. Their product is therefore 
independent of grain size. On the other hand, we may inquire into the 
total over-all stress rela.\ed by the grain-boundary flow, the over-all strain 
being maintained constant. The fraction of stress rela.xed is a dimension- 



Fic. 48.—Variation of internal friction with temperature in polycrystalline and “single- 
crystal” aluminum (frequency of vibration = 0.8 cycles per second at room temperature). 
(After K6.) 

less quantity and cannot therefore depend upon grain size, since there is 
no other quantity having dimensions of length upon which the stress re¬ 
laxation can depend, as long as the grain size is small compared with the 
smallest dimensions of the specimen. This nondependence of maximum in¬ 
ternal friction upon grain size is illustrated in Figure 49 from further work 
of Ke.^^* It is seen that, as long as the grain size is small compared with the 
diameter of the specimen, the effect of an increase of grain size is merely to 
raise the temperature for the maximum internal friction. 

An increase in the frequency of vibration has essentially the same effect 
as an increase in grain size, namely, a higher temperature is required for 

142. T. S. K6. “Stress Relaxations across Grain Boundaries in Metals,” ibid., LXXII 
(1947), 41. 
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Fig. 49.—Effect of grain size on internal friction and rigidity relaxation in aluminum plotted 
against l/F (frequency of vibration = 0.69 cycles per second at room temperature). Average 
grain diameter: G, 0.02 cm.; X, 0.04 cm. (After K6.) 
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the same degree of relaxation. The interrelation of frequency, grain size, 
and temperature is beautifully demonstrated in Figure 50, taken from 
Ke. 

Ke*'*^ has shown that the heat of activation associated with grain bound¬ 
ary slip is essentially identical to the heat of activation for self-diffusion 
and for creep, as is demonstrated in Table 13 taken from his work. Be¬ 
cause of this correlation he has suggested that grain boundary slip involves 
the same mechanism as does volume diffusion. 

ii. Creep and creep recovery .—Shear stress relaxation across grain 
boundaries under constant stress conditions must necessarily be accom¬ 
panied by an extension, i.e., by creep. The creep rate will be smaller, the 
larger the grain size. This viewpoint is in accord with the general recog- 

T.VBLE 13 

Different T\tes of Activ.\tio\ Energy 
FOR Different Metals 

(In Calories/Mole) 


Metal 

Volume 

Diffusion 

Grain 

Boundary 

Slip 

Creep 

a-brass. 

41,700 

41,000 

42,000 

a iron .. 

78,0(K) 

85,000 

90,000 

.■\luminum. .. 

37,500 

34,500 

1 

37,000 


nition*^"* that increase in grain size favors an increased resistance to creep. 
As one example, in his work upon high-purity lead, McKeown^'^^ found 
that the minimum creep rate for a given load was essentially inversely 
proportional to the grain size and hence to the total grain area per unit 
volume. As another e.xample, Burghoff^"*® and collaborators found that the 
creep resistance of brass wire increased with an increase in grain size. 
If creep tests are not made sufficiently slowly, the usual correlation be¬ 
tween grain size and creep rate will not be observed, as was, for example, 
the case in the experiments of Parker and Riisness.'^^ 

143. T. S. K6, “The Structure of Grain Boundaries in Metals,” Phys. Rev., LXXJII (1948), 
267. 

144. H. Carpenter and T. M. Robertson, Metals (London: Oxford University Press, 1939), 
p. 1415. 

145. J. McKeown, “Creep of Lead and Lead Alloys. I. Creep of Virgin Lead,” Jour. Inst. 
Metals, LX (1937), 201. 

146. H. L. Burghoff, A. I. Blank, and S. E. Maddigan.“The Creep Characteristics of Some 
Copper Alloys at Elevated Temperatures,” Trans. A.S.T.M., XLII (1942), 668. 

147. E. R. Parker and C. F. Riisness, “Effect of Grain Siae and Bar Diameter on Creep 
Rate of Copper at 200“ C,” Trans. A.I.M.E., CLVT (1944), 117. 
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An increase in creep rate with a decrease in grain size will be obtained 
only if the temperature is sufficiently high, as well as the rate of creep suf¬ 
ficiently low, that the grain boundaries manifest a smaller resistance to 
deformation than does the grain interior. Thus, if creep tests are made at 
a set of temperatures, a critical temperature will be found below which 
small-grained specimens are more creep resistant and above which large¬ 
grained specimens are more creep resistant. E.xamples have been cited by 
Clark and White."'^ 

The relaxation of shear stress across grain boundaries must be accom¬ 
panied by a readjustment of stresses within the interior of the grains. Such 
stress readjustment will result in residual stresses of such a sign as to pro¬ 
duce slip across the grain boundaries in an opposite direction to that 
caused by the original load. In other words, these residual stresses will 
give rise to an at least partial creep recovery. Shear stress relaxation across 
grain boundaries as one possible source of creep recovery was recognized 
by Jeffries and Archer"^ in their discussion of the creep and marked creep 
recovery of a very fine-grained ternaiy' eutectic. A demonstration that the 
creep occasioned by viscous slip across grain boundaries is essentially 
completely recoverable has been given by Ke,***® whose observations are 
reproduced in Figure 51. 

iii. Stress relaxation .—We have seen (p. 50) that, as long as the stress 
level is sufficiently low that the superposition principle is valid, the 
stress relaxation function is essentially the reciprocal of the creep func¬ 
tion. One therefore anticipates stress relaxation at constant deformation 
under the same conditions as for creep and expects that viscous slipping at 
grain boundaries will be responsible for stress relaxation when, under the 
same conditions, it is responsible for creep. Numerous stress relaxation 
experiments have been conducted in recent years.The tests have been 
made on commercial materials, with no attempt to understand the physi¬ 


cal origin of the stress relaxation. 
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The anticipated reciprocity of stress relaxation and of creep associated 
with grain-boundary slip has, in fact, been found by Previously, 

West^®" had conducted stress relaxation experiments in iron and had in¬ 
terpreted his results in terms of grain boundary viscosity. These interpre¬ 
tations have been questioned*^^ by Ke, who has observed that in iron a 
carbon content such as that contained in West’s specimens is sufficient to 
block grain-boundary slip. 

iv. Elastic moduli .—It has long been recognized that the elastic moduli 
decrease with an increase of temperature. In 1891 Sutherland‘S^ found that 
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Fig. 51.—Creep under constant stress and creep recovery at 175® C. in polycrystalline 
aluminum. (After K6.) 

the data then existing for the rigidity modulus, G, of nonferrous metals 
were in agreement with the formula 



where Go is the rigidity modulus at the absolute zero and is the melting 
temperature. Although his measurements extended to a T/T^ of only 
0.75, Sutherland was sufficiently confident of the extrapolation to r« to 
believe that the rigidity approaches zero as the melting temperature is 

152. W. A. West, “Elastic After-effects in Iron Wires from 20® to 550® C.,” Trans. A J M E 
CLXVII (1946), 192. 

153. T. S. K6, personal communication. 
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Phil. J/oj-., XXXII (1891), 31. 
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approached. Later observations by Jasper*®^ were found to be in agree¬ 
ment with Sutherland’s formula. The data for iron and steel, which have 
been reviewed by Sisco,show the same general trend as Sutherland’s 
formula. The concept that the rigidity modulus becomes essentially zero 
as the melting temperature is approached has received rather widespread 
acceptance, as may be seen from the discussion of Brillouin.^^" 

Whether the observed rapid decrease of the rigidity modulus at high 
temperatures is a characteristic of the interior of the grains or is merely a 
manifestation of stress relaxation across grain boundaries can be decided 
only by a comparison of the temperature variation of the rigidity modulus 
of single crystals and of polycry-stalline specimens. The work of Hunter 
and Siegel'^^ upon a single crystal of NaCl, and of Siegel and Cummerow*^® 
upon a lead crystal, indicates that the latter is the case. Thus in the former 
example the rigidity modulus, Cu, decreased only 28 per cent as the speci¬ 
men was raised from room temperature to just below melting temperature. 
The recent e.xperiments of Ke‘®” upon single crystalline and polycrystal¬ 
line aluminum, shown in Figure 52, demonstrate conclusively that the 
usually observed rapid decrease of the rigidity modulus at high tempera¬ 
tures is, in fact, due to viscous slipping along the grain boundaries. The 
dynamic rigidity modulus of the single crystal decreased slowly in a linear 
manner as the temperature was increased to 400° C., while that of the 
polycrystalline specimen departed rapidly from this linear relation at 
about 200° C. Prior measurements by Birch and Bancroft*®* on the dynani- 
ic shear modulus of polyciy'stalline aluminum agree precisely with Ke’s 
measurements on a single crystal up to 400° C., departing from a linear 
relation only above this temperature. Such a correlation is understood 
when cognizance is taken of the high frequency used by Birch and Ban¬ 
croft, 8,000-5,000 cycles; sec in contrast to a frequency of about 1 cyde/sec 
used by Ke. The temperature at which grain boundary relaxation be¬ 
comes appreciable is raised from about 250° to 450° C. by this increase in 
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frequency. Similar observations have been reported by Ke for grain 
boundary^ relaxations in the case of a-brass’®'* and in the case of a-iron.^®^^ 
A computation has been made*®- of the maximum effect that stress re¬ 
laxation across grain boundaries can have upon Young’s modulus in the 
particular case in which the grains are regular polyhedra of identical size. 
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Fig. 52.—Relaxation of rigidity modulus in polycryslalUne aluminum. (/Vfter K6) 

The ratio of the completely relaxed modulus, Er, to the unrelaxed modu¬ 
lus, Euj was found to vary from 0.62 to 0.68 over the usual range of Pois¬ 
son’s ratio of 0.3-0.4, From the analysis on pages 56-57 we saw that the 
relaxation in the shear modulus is about 15 per cent greater than the re¬ 
laxation in the Young’s modulus, so that the ratio Gr/Gu should lie within 
the range 0.56-0.63 for the usual range of values of Poisson’s ratio. If the 

161a. T. S. K6, “Viscous Slip along Grain Boundaries and Diffusion of Zinc in Alpha 
Brass/’ J. Appl. Phys., XIX (1948), 285. 
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grains are not regular polyhedra or of uniform grain size, the relaxation 
will be greater, and Gr/Gu may approach zero. The experimental value of 
the ratio Gr/Gu in equiaxed aluminum may be seen from Figure 52 to 
be in complete agreement with the theoretical value. 

V. Fracture .—As discussed on pages 126-30, the relaxation of shear 
stress in an isolated region will lead to high stress concentration, which 
may result in the formation of small microcracks. An example of how 
shear stress relaxation across a grain boundary may lead to a microcrack 
is illustrated in Figure 53. Here tensile stress is taken as applied along an 
axis normal to the grain boundary, BB\ This boundary is intersected by 
another grain bounda-y, AA\ across which the tensile stress has a shear 



Fio. 53.—Illusiration of how microcracks can f)e formed by shear stress relaxation along 
the grain boundaries. 


stress component. Relaxation of shear stress across the boundary A A ‘ will 
result in an extremely high tensile stress in the corner AA^BK When this 
tensile stress reaches the intrinsic fracture strength, a small microcrack 
will form. Further slip along the plane AA^ will result in a spread of the 
microcrack until the crack has such a size that it may propagate under 
its own stress concentration. Further propagation need not then wait for 
further slip along the plane AA^ but proceeds at a great rate. Fracture of 

the specimen ensues. 

The role played by the viscosity of grain boundaries at elevated tem^- 
peratures in the initiation of creep was first recognized by Rosenhain'*'" 
and his colleagues in their work upon the deformation of 7 -iron at high 
temperatures. When pulled under the usual test conditions, this iron mani- 
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fested fair ductility and fractured in a fibrous manner. When extended 
very slowly at the same temperature, thereby allowing time for viscous 
slip, the iron fractured in an intercrystalline manner, with hardly any 
prior plastic deformation. The occurrence of this type of intercrystalUnc 
failure in a wide range of metals was later discussed by Rosenhain and 
Archbutt,who demonstrated that such failure is facilitated by small 
plane-grain boundaries, grain boundaries that can slip without mechanical 
hindrance. Season cracking in brass was ascribed by them to grain bound¬ 
ary slip. More recently Hanson and Wheeler'®^ have made a study of the 
role that grain boundaries play in the fracture of aluminum under creep 
conditions. They found that at elevated temperatures cracks invariably 
started at the grain boundaries under slow extension and propagated into 
the grains only just before fracture, when the rate of extension was large 
and hence when the resistance to deformation of the viscous grain bound¬ 
aries became greater than that of the grain interior. As spectacular evi¬ 
dence of crack formation at grain boundaiies, they showed that under 
creep conditions the density of polycrystalline specimens decreased ap¬ 
preciably, while that of single crystals remained constant. A beautiful ex¬ 
ample of intercrystalline cracks in copper found by slow extension at ele¬ 
vated temperatures has been given by McAdam‘®® and is reproduced 
as Figure 54. 

It is well recognized that the conditions which favor viscous flow at 
grain boundaries—high temperature and low rates of extension—also 
favor low ductility and that the resulting essentially brittle failure is even 
more deleterious than creep under service conditions. As pointed out by 
Carpenter and Robertson,’®’ those factors which increase the resistance of 
metals to creep do not necessarily increase their resistance to brittle failure 
by the same amount. 

IV. STRESS RELAXATION ACROSS TWIN INTERFACES 

Twins induced by stress usually form by a discontinuous process, their 
appearance being accompanied by a characteristic clicking sound. Once a 
twin band has formed, however, its boundaries may move in a continuous 
manner. A well-known example is furnished by recrystallized a-brass. In 
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this metal each crystal is traversed by several twin bands, the number of 
twin bands, and the relation of their width to the diameter of the individ¬ 
ual grains being essentially independent of grain size. Apparently, as the 
grains grow, the width of the twin bands likewise grows by a continuous 
displacement of the twin interface boundaries. 

As may be seen by reference to Figure 42, a shear stress across a twin 
interface tends to induce a movement of this interface. Thus in the case 
illustrated, a mass movement of the right-hand twin member in the direc¬ 
tion of the applied traction is accompanied by the movement of the twin 



Unetched (X250) Etched with a mixture of equal parts of concen¬ 

trated nitric and glacial acetic acids (X250) 


Fig. 54.—Longitudinal sections after creep tests of cold-drawn monel metal 

interface from AA' to BB\ Conversely, twin interfaces will tend to adjust 
their positions in such a manner as to minimize the shear stress acting 
across them. This continual striving of the twin interfaces to minimize the 
shear stresses will give rise to characteristic anelastic effects. 

The anelasticity of copper manganese alloys (~90 per cent Mn) has 
been studied by F. Worrell,^®* and is thought to arise from the stress-in¬ 
duced movement of the twin interfaces of this alloy. The temperature 
variation of the internal friction of a typical specimen is given in Fipire 
55. In the low-temperature range the motion of the twin interfaces is so 
slow that no appreciable displacement occurs during a half-cycle of vibra¬ 
tion, and so the internal friction is very low. On the other hand, in the 
high-temperature range the twin interfaces are so mobile that the shear 

168. F. Worrell, Pkys. Rev., LXXII (1947), 533. 
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stress is completely relieved at all times, so the internal friction is again 
very low. Only in the intermediate-temperature range, where both the dis¬ 
placement and the net shear stress are appreciable, is the internal fric¬ 
tion large. 


Temperoture (*C) 

50 25 0 -25 .50 



Fig. 55. —Typical variation of internal friction of CuMn alloy with temperature. This 
internal friction is supposedly due to movement of twin interfaces. (After Worrell.) 


The alloys of manganese with f.c.c. metals are especially suited to a 
study of the laws governing the movement of twin interfaces. It is com¬ 
monly thought that the high-temperature 7 -phase of manganese is 
face-centered tetragonal (f.c.t.) and that the tetragonality ratio, c/a, in¬ 
creases gradually from 0.934 for pure Mn to unity as the concentration of 
a f.c.c. metal is increased. As yet no high-temperature X-ray structures 
have been determined, the f.c.t. structure being deduced from measure- 
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ments on quenched specimens. The highly twinned nature of the CuMn 
alloys, first discovered by Worrell and illustrated in Figure 56, suggest 
that the high-temperature structure is, in fact, cubic and that the tetrag- 
onality is induced by a diffusionless transition during quenching. Accord- 



Fig. 56.—Example of finely spaced (101) twins of face-centered tetragonal CuMn alloy 
superposed upon twins characteristic of a face-centered cubic lattice. (After Worrell.) 

ing to this viewpoint, heavy twinning must arise in order that high resid¬ 
ual stresses may be avoided. Whatever the origin of the twins in the man¬ 
ganese alloys, they are of unusual theoretical interest because of the small¬ 
ness of the atomic movements associated with the motion of the twin in¬ 
terfaces. Thus, in the interface movement illustrated in Figure 40, the 
shear strain between the planes AA' and BB' varies from 0.066 in pure 
"y-manganese to less than 0.015 in 7 -manganese containing about 18 per 
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cent copper. In the latter case the atomic displacements associated with 
the movement of a twin interface by one interatomic distance are of the 
same order of magnitude as the amplitude of vibration at the absolute 
zero temperature. 

A second example of stress-induced movement of twin interfaces occurs 
in the case of magnetic domains. In a stress free specimen of iron, the local 
magnetization is along one of the principal axes, resulting in local tetrag- 
onality through magnetostriction. The boundary of two adjacent do¬ 
mains is along a (110) plane'®® and thus may be regarded as a twin inter¬ 
face, the two domains joining with no elastic distortion whatsoever. The 
application of a stress visibly moves these interfaces, and the time lag 
in this movement gives rise to anelastic effects. Since the tetragonality in 
the magnetic domains is very small, the c/a ratio being 1.000032,^"® a very 
small amount of stress relaxation results in a large movement of the twin 
interface. 
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